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2.5.1 

*  2. 5. 1,1 

Soils  Survey  and  Productivity  Assessment  Studies 

Objectives 

Objectives  of  the  soils  program  remain  imchanged.  Consult 
Progress  Report  2,  Section  2.  5. 1. 1  for  specifics. 

2.  5. 1.2 

Methods 

Methods  employed  in  the  soil  studies  are  described  in 
Section  2.  5. 1.  2  of  Progress  Report  2. 

2. 5. 1.3 

Results 

Soil  surveying  and  mapping  is  being  carried  out  by  the 

Soil  Conservation  Service  (SCS)  according  to  standard 
techniques.  The  SCS  has  completed  trace  element  and 
soil  horizon  analyses  of  two  pits  in  the  area.  A  soil 
map  has  been  prepared  and  will  be  finalized  and  become 
available  in  June. 

A  soils  contractor  will  be  selected  in  a  few  weeks  and 
this  phase  of  the  program  will  be  initiated  in  July.  The 
program  will  include  the  collection  and  analyses  of 
soils  samples  associated  with  major  vegetation  types. 

Trace  metal  concentrations  will  be  determined  and 
soil-plant  relationships  will  be  interpreted. 
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2.  5.  2. 
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2. 5. 2.1. 

Archaeological  Survey 

Objectives 

The  archaeological  program  is  designed  to  prevent  the 
inadvertent  loss  of  valuable  historical  or  archaeological 
data  as  development  of  the  area  progresses. 

2. 5. 2. 2. 

Methods 

Methods  by  which  archaeological  surveys  are  performed 
are  described  in  Progress  Report  2,  Section  2.  5. 2. 2. 

2.  5.  2.  3. 

Results 

Several  potential  development  areas  were  examined  for 
the  presence  of  archaeological  or  historical  artifacts 
during  the  study  period.  They  included  proposed  plant 
sites  2  and  3,  shallow  alluvial  water  monitoring  test 
holes  S-5,  S-16,  S-17,  S-18,  S-19,  S-22  and  S-23,  an 
unnumbered  location  immediately  south  of  84  Ranch 
monitoring  hole  M-5.  Clearances  were  granted  for  each 

I 
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2.  5.  3  Re  vegetation 

2.  5.  3. 1  Objectives 

Objectives  of  the  revegetation  program  are  discussed  in 
Section  2.  5.  3. 1  of  Progress  Report  2.  The  revegetation 
research  is  designed  to  provide  sufficient  information  to 
specify  rehabilitation  plans. 


2.  5.  3.  2  Methods 

Native  vegetation  on  the  selected  sample  plots  will  be 
completely  removed  before  seeding.  Seed  mixtures  will 
be  planted  at  each  plot  over  a  3-year  period  starting  in 
the  fall  of  1975.  Discing  on  the  various  sites  will  be 
carried  out  during  early  July,  if  possible,  so  that  the 
seedbed  can  settle  before  fall  seeding. 

2.  5.  3.  2. 1  Selection  of  Species 

Species  have  been  selected  on  the  following  criteria: 
drought  resistance,  palatability,  growth  response,  growth 
form,  ability  to  stabilize  soils  and  origin.  Extremely 
palatable  species  (for  either  livestock  or  big  game)  will  be 
avoided  to  prevent  excessive  grazing  pressure  on  the 
revegetated  areas. 

The  species  listed  in  Table  2.  5.  3  -  1  are  suggested 
because  of  their  wide  suitability  to  various  site  conditions 
found  on  Tract  C-a.  Plant  responses  will  indicate  optimal 
locations  for  revegetation  on  Tract  C-a.  The  seeds 
selected  provide  a  diversity  of  plant  forms  including  grasses, 
forbs  and  shrubs. 


2.  5.  3.  2.  2  Sample  Plots 

Each  treated  plot  will  be  90  x  40  m,  with  at  least  a  1-m 
buffer  zone  around  the  entire  plot. 

Two  conventional  means  of  planting  (drilling  and  broad¬ 
casting)  will  be  used  to  place  the  seed  in  the  surface  soil. 
A  small  row  seed  planter  or  a  small  single  or  double  disc 
grain  drill  equipped  with  both  small  and  large  seed  boxes 
will  be  used  for  drilling.  Rows  will  be  35.  5  cm  apart  and 
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Table  2.  5.  3  -  1.  Suggested  list  of  species  to  be  seeded  on  experimental 

revegetation  plots  on  and  near  Tract  C-a  for  the  RBOSP. 


Common  Name 

(Scientific  Name) 

Suitability 

Code'^ 

Amount  of  Seed 
(kg/ha) 

Drilled  Broadcast 

Pubescent  wheatgrass 

(Agropyron  trichophorum) 

2,  4,  6,  8,9 

2.  80 

5.  62 

Western  wheatgrass 

(Agropyron  smithii) 

4,  5,8,9 

2.80 

5.  62 

Bluebunch  wheatgrass 

(Agropyron  spicatum) 

1,5,7 

2.80 

5.62 

Streambank  wheatgrass 
(Agropyron  riparium) 

3,  5,  8,9 

2.  80 

5.  62 

Sand  dropseed 

(Sporobolus  cryptandrus) 

2,  4,  5,  7,  9 

0.  56 

1.12 

Indian  ricegrass 

(Oryzopsis  hymenoides) 

2,4,  5,  7,9 

1.  12 

2.25 

American  vetch 

(Vicia  americana) 

5 

0.  56 

1.  12 

Cicer  milkvetch 

(Astragalus  cicer) 

6,  11 

0.  56 

1.  12 

Yellow  sweetclover 

(Melilotus  officinalis) 

2,4,5,9,11 

0.56 

1. 12 

Big  sagebrush 

(Artemisia  tridentata) 

2,5,  9,10,11 

0.22 

0.45 

Little  rabbitbrush 

(Chrysothamnus  vicidiflorus) 

2,4,5,9,11 

0.22 

0.45 

Bitterbrush 

(Purshia  tridentata) 

1,  5,  9, 10,  11 

1.  12 

2.25 

Mountain  mahogany 

(Cercocarpus  montanus) 

1,  5,  10 

0.  56 

1.  12 

Fourwing  saltbush 

(Atriplex  canescens) 

3,  4,  5,9 

1.  12 

2.25 

Totals 

17.  80 

35.  73 

1  -  Highly  palatable;  2  -  Moderately  palatable;  3  -  Unpalatable;  4  -  Soil 
stabilizer;  5  -  Native  species;  6  -  Introduced  species;  7  -  Bunch  grass; 

8  -  Sod  grass;  9  -  Drought  resistent;  10  -  Browse;  11  -  Cover  for  game 
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40  m  in  length.  A  whirlwind  seeder  which  broadcasts 

#  larger  seeds  in  one  operation  and  smaller  seeds  in  a 

second  operation  will  be  used  to  broadcast  seeds.  All 
broadcasted  seeds  will  be  covered  lightly  with  soil. 

Three  mulches,  hay,  wood  fiber  (silvifiber)  and  crushed 
rock,  will  be  applied  as  treatments  to  all  seeded  plots  in 
each  of  the  locations  to  measure  their  effectiveness  in 
aiding  seedling  emergence  and  survival.  All  mulches 
will  be  applied  to  cover  approximately  70%  of  the  soil 
surface  as  determined  by  a  point  hit  plot  method.  A 
balanced  fertilizer  will  be  applied  on  all  seeded  plots 
after  the  first  growing  season  at  the  rate  of  80  pounds  per 
acre. 

Each  subplot  sample  will  be  10  x  10  m.  This  will  provide 
two  samples  1  m  long  with  rows  .  5  m  wide  apart  and 
equivalent  space  in  broadcasted  areas.  Two  replicates 
will  be  established  and  permanently  marked  for  subsequent 
yearly  samplings. 

2.  5.  3.  2.  3  Parameters  to  be  Measured 

f  During  the  first  spring  following  fall  planting,  the  number 

of  emerged  seedlings  will  be  counted  per  plot  for  each 
species. 

At  the  end  of  the  second  growing  season,  the  number  of 
seedlings  will  be  counted  on  the  original  plots  to  determine 
seedling  survival  for  each  species. 

At  the  end  of  the  third  growing  season,  production  of  total 
above-ground  biomass  will  be  determined.  Total  herbage 
of  seeded  and  alien  species  will  be  determined  separately. 
Production  of  the  major  seeded  species  making  up  the 
majority  of  the  total  herbage  yield  will  be  recorded  by 
species.  Along  with  production  of  biomass  of  the  seeded 
species,  closed  community  ratings  will  be  ranked  for  each 
treatment.  These  rankings  will  range  from  1  to  10,  with 
1  being  the  best  soil  cover  and  10  being  the  poorest.  The 
criterion  used  for  rankings  will  be  the  size  of  interspaces 
among  seeded  perennial  plants.  If  spaces  are  1  m2  or  less, 
the  rating  would  be  1;  if  the  interspaces  are  as  large  as 
5  m2,  the  rating  would  be  5;  and,  if  the  interspaces  are  as 
large  as  10  m2,  the  rating  would  be  10. 
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2. 5. 3. 2. 4 


Treatments 


§ 


Nine  individual  treatments  are  suggested: 

1.  Control.  No  seeding,  no  mulches  and  no 
fertilizer  (natural  succession). 

2.  Drilled  seeded. 

3.  Drilled  seeded  with  hay  mulch. 

4.  Drilled  seeded  with  wood  fiber  mulch. 

5.  Drilled  seeded  with  crushed  rock  miilch 

(2  to  5  cm  diameter). 

6.  Broadcast  seeded  with  soil  coverage. 

7.  Broadcast  seeded  with  soil  coverage  and 
hay  mulched. 

8.  Broadcast  seeded  with  soil  coverage  and 
wood  fiber  mulched. 

9.  Broadcast  seeded  with  soil  coverage  and 
crushed  rock  mulched  (2  to  5  cm  diameter). 

2.  5.  3.  2.  5  Experimental  Matrix 

The  sample  plots  will  be  designed  as  a  matrix  with  nine 
treatments  and  four  replicates  of  each  treatment.  Three 
40  X  90  m  sample  plots  shall  be  established.  Each 
treatment  within  every  large  plot  shall  constitute  an 
experimental  unit.  An  additional  comparison  shall  be 
made  between  drilled  and  broadcasted  planting  techniques. 
Treatment  placement  shall  be  deterministic  rather  than 
random;  therefore,  Model  I  factorial  analysis  of  variance 
shall  be  used.  Three  parameters  will  be  compared 
(treatment,  species  and  planting  techniques).  A  three-way 
analysis  of  variance  will  be  used.  Units  of  interest  shall 
be  for  a  given  species  within  a  given  treatment  and  a  given 
planting  method.  Analysis  of  variance  tables  shall  have: 
a  treatment,  a  species,  a  treatment  x  species  interaction, 
a  planting  method,  a  treatment  x  planting  method  inter¬ 
action  and  species  x  planting  interaction,  as  well  as  the 
error  or  residual  term.  All  analyses  shall  be  tested  at  the 
95  confidence  level.  Parameters  shall  be  compared  using 
the  Newman-Keuls  multiple -range  technique  to  minimize 
Type  II  error. 

2.  5.  3.  2.  6  Schedule 

Each  year  for  3  years,  seeding  will  be  accomplished 
during  the  fall.  Treatment  procedure  will  be  the  same 
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each  year  unless  results  warrant  a  change.  Therefore, 

^  the  first  yearns  seeding  would  yield  data  during  the 

3 -year  study  on  seedling  emergence,  seedling  survival, 
herbage  yield  and  closed  community  ratings.  The 
second  year’s  seeding  would  provide  data  on  seedling 
emergence  and  seedling  survival  and  the  third  year’s 
seeding  would  yield  data  only  on  seedling  emergence. 

2.  5.  3.  2.  7  Associated  Studies 

Percent  soil  moisture  will  be  determined  by  means  of  a 
soil  probe  at  15,  30  and  60  cm  depths.  Three  samples 
will  be  collected  at  each  depth.  These  determinations 
will  be  made  during  the  fall  at  the  time  of  seeding,  in  the 
spring  when  moisture  is  most  available,  and  again  in 
mid-July  when  soil  moisture  is  least  available.  Soil 
moisture  will  be  determined  only  for  the  level  plots  that 
serve  as  controls  and  will  be  made  for  each  of  the  3  years 
of  the  study. 

Sufficient  area  at  each  location  will  be  prepared  at  one 
time  to  satisfy  the  needs  for  plantings  to  be  made  during 
the  entire  study,  if  possible.  The  area  not  seeded  each 
#  year,  except  for  the  control  plots,  will  be  fallowed  in 

the  summer  to  control  plant  growth. 

Initially,  the  native  vegetation  will  probably  be  removed 
from  the  area  by  light  blading  with  a  minimum  disturbance 
to  the  topsoil,  followed  by  several  deep  plowings  or 
discings.  Deep  plowing  or  discing,  followed  by  rotary 
sweep  raking  to  remove  the  debris  from  the  area  prior  to 
planting,  may  be  used  as  an  alternative.  The  soil  will  be 
thoroughly  disturbed  and  mixed  to  a  depth  of  45  cm. 

The  soils  at  each  site  will  be  analyzed  for  texture,  pH, 
exchangeable  sodium  percentage,  electrical  conductivity 
and  organic  matter.  Two  samples  will  be  taken  from 
each  site  at  0  to  15  cm,  18  to  30  cm  and  33  to  61  cm  for 
analysis. 
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2.5.5 


Trace  Metals 


Trace  metal  analyses  were  performed  on  soils  from  two 
SCS  pits  during  the  quarter.  Results  of  these  analyses 
will  be  available  within  a  few  weeks.  Analyses  of  plant 
and  animal  tissues  will  be  performed  if  interpretations 
of  the  soils  data  indicate  that  trace  metal  concentrations 
in  the  soils  are  sufficiently  high  enough  to  warrant 
further  investigation. 
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Seismicity 
Objectives 

The  purpose  for  performing  a  seismicity  investigation  of 
the  Piceance  Creek  Basin  in  which  Tract  C-a  lies  was  to 
determine  the  extent  of  seismic  activity  near  the  site  due 
to  earthquakes. 

Methods 

To  develop  the  information  necessary  for  this  study,  a 
literature  search  was  made  to  establish  both  recorded 
and  reported  events  in  the  area.  Contact  was  made  with 
various  governmental  and  private  organizations  to 
gather  the  relevant  data. 

Results 

The  following  information  was  submitted  by  Dr.  Alcock 
at  the  completion  of  his  literature  review. 

"The  Piceance  Creek  Basin  in  which  oil  shale  Tract  C-a 
lies  appears  to  be  an  area  of  very  low  seismicity. 

a.  Earthquake  Distribution  -  The  seismic  risk  map 

(Ref.  1)  shows  that  all  of  the  Piceance  Creek  Basin 
is  a  quiet  seismic  area  and  is  classified  as  zone  1. 
Structure  in  this  zone  with  fundamental  periods  greater 
than  1.  0  sec.  can  be  expected  to  suffer  no  more  than 
minor  damage  from  large  distant  earthquakes.  This 
opinion  of  the  seismic  activity  of  the  area  is  also 
held  by  Mrs.  Ruth  Simon  (Ref.  2)  and  Dr.  Frank  Stead 
(  personal  communication  )  both  of  the  U.  S.  Geological 
Survey.  Examination  of  the  Seismic  Activity  of 
Colorado  maps  issued  annually  between  1966  and  1971, 
(Fig.  1,  Ref.  3)  shows  that  the  only  earthquakes  in 
Colorado  within  40  miles  of  the  Tract  C-a  (magnitude 
2  or  greater)  were  at  the  Harvey  Gap  Reservoir  and 
at  Rangely.  The  last  available  map,  Fig.  2,  shows 
the  same  seismicity  for  the  Basin.  Those  events  at 
the  Harvey  Gap  Reservoir  are  believed  to  have  been 
due  to  the  effect  of  the  loading  by  the  stored  water 
(Ref.  2)  on  the  underlying  geologic  structure.  Those 
events  at  Rangely  have  been  postulated  to  result  from 
%  the  water  flooding  of  the  Rangely  oil  field  (Ref.  4). 
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A  search  was  made  of  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  historical 
seismic  data  in  an  area  within  50  and  100  miles 
of  the  C-a  tract.  This  search  covered  the  period 
from  1941  to  1975  for  earthquakes  greater  than 
magnitude  3.  The  results  of  this  search  are  shown 
in  Figure  3.  (A  computer  printout  and  legend  of 
the  search  of  the  earthquake  data  file  appears  as 
Reference  9. ) 

CER  Geonuclear  Corp.  operated  a  seismic  station 
in  connection  with  Project  Rio  Blanco  from 
November  20,  1971  to  December  17,  1973.  This 
station  was  located  in  Sec.  28,  T2S,  R98W,  only 
about  10  km  from  Tract  C-a  (Ref.  5).  In  this  time 
interval,  aside  from  the  earthquake  and  aftershocks 
resulting  from  the  detonation  of  Project  Rio  Blanco, 
only  2  microearthquakes  (local  magnitude  less  than 
1)  were  recorded  with  40  km  of  the  station. 

In  the  preparation  for  the  seismic  investigation  of 
the  Rangely  oil  field  the  U.  S.  Geological  Survey 
made  a  search  of  the  seismic  data  recorded  for  a 
seven  year  period  at  the  Uinta  Basin  Observatory 
for  earthquakes  originating  near  the  oil  field,  65  km 
ESE  of  the  observatory  (Ref.  6).  This  search  covers 
the  area  in  which  Tract  C-a  lies  and  shows  no  more 
than  one  earthquake  of  local  magnitude  greater  than  1 
in  the  seven  year  period.  This  seems  to  be  in  close 
agreement  with  the  results  from  the  CER  seismic 
station. 

b.  Crustal  Stress  -  The  strain  release  map  of  the  area 
(Ref.  1)  indicates  very  little  probability  of  any  strain 
stored  in  the  subsurface  to  cause  earthquakes.  In 
preparation  for  the  Rio  Blanco  Project,  I  made 
several  calculations  of  subsurface  stress  existing  at 
various  points  in  the  Rio  Blanco  unit  area,  using 
Kehle’s  method  of  determining  the  stress  from 
hydraulic  fracturing  data  (Ref.  7).  In  all  cases  the 
principal  stress  was  compressive  and  vertical  and 
within  the  limits  of  the  method  equal  to  the  lithostatic 
pressure.  There  was  no  evidence  of  any  significant 
horizontal  stress. 
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c.  Recommendations  -  In  view  of  the  available  data 
g  quoted  above,  I  believe  that  the  seismicity  of  the 

oil  shale  Tract  C-a  is  very  low  and  is  amply 
documented.  Sanford  &  Singh  (Ref.  8)  have  shown 
than  in  general  about  six  months  of  recording  will 
give  a  reliable  estimate  of  the  earthquake-frequency 
relation  from  which  short  term  (10  yrs. )  seismicity 
is  calculated. 

I  believe  that  additional  seismic  monitoring  at  the 
present  time  is  redundant  and  could  be  postponed 
until  6  or  8  months  in  advance  of  actual  mine 
operations.  At  that  time  a  review  of  the  situation 
could  change  the  seismic  program  due  to  the  mining 
techniques  to  be  used  and  the  specific  objectives  of 
the  seismic  program,  i.  e.,  environmental  safety, 
mine  safety,  etc.  " 


Respectfully  submitted, 

(Original  signed  by 
Dr.  Alcock) 

Dr.  E,  D.  Alcock 
Registered  Progessional 
Engineer 
Colorado 

Copies  of  the  references  referred  to  in  this  letter  appear 
in  the  appendix. 
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FIGURES 


Figur©  1  Fig.  IC-1  of  CER  EIE  PNER-B4 

Figure  2  Seismicity  Map  for  1972 

Figure  3  Fig,  IC-2  of  CER  EIE  PNER-B4 
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Dr.  E.  D.  Alcock 


Geophysicist 


Dr.  Alcock  holds  a  B.  S.  in  Mathematics  (1930)  from  the  University  of 
California  at  Los  Angeles  and  a  M.  S.  (1933)  and  Ph.  D.  (1935)  in  Physics 
from  the  California  Institute  of  Technology. 

Dr.  Alcock  was  employed  by  the  Bell  Telephone  Laboratories  as  a 
research  engineer  working  on  accoustical  problems  from  1930  to  1932. 

From  1935  to  1965,  he  was  employed  by  the  oil  and  gas  industries, 
gaining  experience  in  seismic  exploration,  drilling  and  production  of  oil 
and  gas  wells,  property  evaluation  and  oil  company  management.  From 
1965  to  1968,  Dr.  Alcock  conducted  basic  research  for  the  Atomic  Energy 
Commission  on  the  effect  of  geology  on  seismic  response.  He  joined  CER 
Geonuclear  in  1970  and  has  been  actively  engaged  in  measuring  and  evaluat¬ 
ing  the  effects  of  underground  nuclear  explosives  on  the  environment. 
Specifically,  he  was  responsible  for  the  evaluation  of  the  seismic  effects  of 
Project  Rulison.  He  also  designed  and  fielded  the  seismic  network  at 
Project  Rio  Blanco  and  analyzed  the  results.  As  an  outgrowth  of  the 
information  received.  Dr.  Alcock  has  devised  a  microtremor  survey 
system  for  evaluating  existing  structures,  building  sites  and  building 
designs  for  possible  earthquake  damage. 

He  is  a  member  of  the  Society  of  Exploration  Geophysicists,  Canadian 
Society  of  Exploration  Geophysicists  and  the  Seismological  Society  of 
America.  He  is  a  Registered  Professional  Engineer  and  a  Registered 
Professional  Geophysicist. 

Dr.  Alcock  has  published  several  papers  in  the  area  of  geophysics  and 
seismicity. 
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Seismic  activity  of  uencraL  area. 
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SEISMIC  RISK  STUDIES  IN  THE  UNITED  STATES 

•  1  * 

S.  T.  Algermissen 

*  i 

Synopsis,  •  i 

A  new  Seismic  Risk  Map  of  the  United  States  is  presented, 
together  with  strain  release  and  maximum  Modified  Mercal.li 

intensity  maps  of  the  country.  Frequency  of  occurrence  of  ; 

damaging  earthquakes  was  not  considered  in  assigning  ratings  to  i 

the  various  zones  on  the. risk  map,  but  studies  of  earthquake 
frequency  are. included  as  an  aid  in  using  the  risk  map.  The 
Seismic  Risk  Map  is  suggested  as  a  revision  of  the  Seismic^ 

Probability  Map  prepared  by  the  Coast  and  Geodetic  Survey  in 
1947  and  withdrawn  in  1952. 


# 


I 


I  ntroductlori 

The  zoning  of  the  United  States  for  seismic  risk  l^s  not 
received  the  attention  in  the  earthquake  engineering  literature  ^ 

that  the  subject  has  enjoyed  in  other  parts  of  the  world,  notably  . 

in  countries  such  as  Japan,  the.  U.S.S.R.  and  New  Zealand.  Risk 
studies  in  the  United  States  have  been  limited  by  the  relatively 
short  historical  record  of  seismicity  available  for  the  country 
as  compared  with  such  countries  as  Japan,  the  difficulty  in 
correlating  known  seismicity  with  geological  evidence  of  tectonic  ^ 

activity  in  many  areas  and  the  general  scarcity  of  pertinent-  .  ' 

geodetic  data  that  might  bear  on  the  problem. 


Materials  Used  t 

The  historical  record  of  earthquakes  in  the  conterminous  ; 

United  States  consists  of  data  on  approximately  28,000  shocks.  ^ 

Within  the  past  two  years  the  Coast  and  Geodetic  Survey  has 
catalogued  all  of  the  pertinent  seismicity  information  on 

punched  cards  and  magnetic  tape.  All  earthquakes  have  been  i 

assigned  a  reference  number  and  origin  times  and  hypocenters 

have  been  recorded  if  available  from  instrumental  data.  Hypocenters 


% 


^Chief,  Geophysics  Research  Group,  Coast  and  Geodetic  Survey,  ESSA, 
U.  S.  Department  of  Commerce,  Rockville,  Maryland,  U.  S.  A. 
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are  estimated  from  accounts  of  the  shock  and  isoseismal  data 
if  no  instrumental  data  are  available.  Magnitudes  and  reported 
intensities  are  also  keypunched  into  each  card.  Sources  of  the 
information,  availability  of  strong  motion  records,  and  other 
related  data  are  also  catalogued.  A  more  complete  description 
of  the  forrrat  of  the  data  stored  on  punched  cards  and  magnetic 
tape  has  been  given  by  Algermissen  (l).  It  is  possible  to 
prepare  quite  a  number  of  different  statistical  presentations 
and  naps  of  the  various  parameters  of  seismicity  using  the 
basic  data  retained  on  magnetic  tape  and  punched  cards.  Two 
types  of  maps  were  prepared  which,  it  is  believed,  are  particu¬ 
larly  useful  in  estimating  seismic  fisk.  These  two.  maps  show 
the  distribution  of  maximum  Modified  Mercalli  (M.M.)  intensities 
and  strain  release  in  the  United  States. 


Maximum  Modified  MercalLi _ 1  ,nten?,it.y. 

The  maximum  M.  M.  intensities  reported  throughout  the  United 
States  from  the  first  recorded  earthquake  in  1534  through  1965 
are  shown  in  Figure  1.  Isoseismal  maps  were  prepared  for  sig¬ 
nificant  earthquakes  not  previously  mapped  by  using  whatever 
data  were  available.  Figure  1  shows  clearly  that  high  M.M. 
intensities  have  been  reported  throughout  a  large  port  ion  of  the 
country,  even  though  the  frequency  of  occurrence  of  high  inten¬ 
sities  in  rrany  areas  may  be  quite  low  . 


Strain  Release  '  ■ 

, 

The  cumulative  strain  release  resulting  from  earthquakes  .*• 

in  the  time  span  1900-1965,  is  plotted  in  Figure  2.  Using  the  ’ 

concept  developed  by  Benioff  (2),  the  strain  release  of  an  ^  , 

earthquake  is  taken  to  be  proportional  to  the  square  root  of  '  i-rH 

its  energy  release.  The  energy,  E,  was  related  to  earthquake'  ,  * 

magnitude,  M,  using  the  formula  log  E  =  11.8  +  1.5  M  (3).  The 
magnitudes  of  earthquakes  below  six  are  known  only  for  shocks  ^ 

occurring  during  approximately  the  past  five  years.  Earthquakes 
whose  magnitudes  are  six  and  above  are  known  from  about  1900. 

When  riB-gnitude  values  were  not  available,  they  were  estimated 
from  maximum  M.M.  intensity  using. the  relation  M  =  1  +  2/3  T, 
where  I  is  M.M.  intensity  (4).  The  above  empirical  relationship 
between  magnitude  and  intensity  is  only  an  approximation  and 
possibly  subject  to  considerable  error.  The  error,  however,  is 
not  large  in  this  calculation  of  strain  release  since  it  was 
necessary  to  use  this  conversion  from  maximum  intensity  to  magnitude 
only  for  small  shocks  for  which  the  strain  release  is  small.  The 
entire  United  States  was  divided  into  squares  of  10,000  square 
kilometers  and  the  strain  release  associated  with  earthquakes  within 
each  square  was  summed.  The  strain  was  represented  as  the  equivalent 
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number  of  magnitude  four  earthquakes  in  each  block.  The  „ 
equivalent  number  of  nmgnitude  four  earthquakes  is  then  obtained 
by  dividing  the  sum  of  the  strain  release  associated  with  earth¬ 
quakes  in  a  particular  block  by  the  strain  release  associated 
with  a  magnitude  four  earthquake. 

The  strain  release  trap  is  useful  in  that  it  shows  the 
current  relative  rate  of  tectonic  activity  in  various  parts 
of  the  United  States.  Areas  of  the  country  with  a  continuing, 
but  a  low  level  of  seismic  activity  are  revealed.  These  areas 
do  not  normally  show  up  on  a  maximum  intensity  mapl  The  exact 
shape  and  pattern  of  the  areas  shown  in  Figure  2  are  not  signifi¬ 
cant  since  the  size  of  the  block  over  which  the  strain  energy  is 
summed  is  somewhat  arbitrary.  A  comparison  of  the  relative  strain 
release  throughout  the  United  States  is,  however,  interesting  and 
i nforrrat i ve,  as  it  is  an  index  of  current  tectonic  activity. 
Inspection  of  Figure  2  shows  that  the  United  States  may  be  divided 
roughly  into  four  areas  that  have  different  rates  of  strain  release 
during  the  66  year  period  considered.  The  four  areas  are? 


Pacific  West 
Rocky  Mountains 
Central  Plains 
Eastern  U.  S. 


-west  of  114  oW.  longitude 

-  106°  -  114  °W.  longitude 

-  92°  -  106°  W.  longitude 

-  east  of  92°  'W.  longitude 


Dividing  the  total  equivalent  number  of  ringnitude  four  earth¬ 
quakes  in  each  area  by  the  total  area  v/e  obtain  the  following 
rates! 


Area 


Pacific  West 
Rocky  Mountains 
Central  Plains 
Eastern  U.  S. 


No.  of  Equivalent  Magnitude _4 

earthquakes/1000  km  ^/66  years 

12 

2.0 

.14 

.74 


Another  calculation  of  the  relative  strain  release  in  the  four 
areas  was  made  considering  only  those  portions  of  each  area  in 
which  more  than  0.25  equivalent  magnitude  four  earthquakes 
occurred.  The  strain  release  per  unit  area  was  computed 
considering  only  the  shaded  portions  of  figure  2.  The  results  are! 


Area 


No.  of  Equivalent  Mapnitude__4 

earthquakes/lOOO  km  ^/66  years 


Pacific  V/est  13 

Rocky  Mountains  2.6 

Central  Plains  *38 

Eastern  U.  S.  1  •  1 
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Using  the  first  set  of  data  on  the  preceding  page,  the  strain 
release  during  the  past  66  years  in  the  Pacific  West  ha's  been 
approximately'six  times  greater  than  in  the  Rocky  Mountain  area, 

86  times  greater  than  in  the  Central  Plains  and  16  times  greater 
than  in  the  Eastern  United  States.  Using  the  secondset  of  ^ta 
on  the  preceding  page,  the  strain  release  in  the  Pacific  West  has 
been,  respectively,-  five,  34  and  12  times  greater  than  the  other 

areas. 

Earthquake  Frequency, 

The  difficulties  encountered  in  any  analysis  of the  recurrence 
rates  of  damaging  earthquakes  are  very  great.  It  has  been  stated 
repeatedly  in  the  literature  that  the  historical  record  of 
seismicity  in  the  United  States  is  inadequate  to  estirrate  recur¬ 
rence  rates.  It  is  believed,  however,  that  the  historical  record 
of  seismicity  in  the  United  States  does  provide  some  guidelines 
to  the  relative  seismicity  of  various  areas  of  the  country. 

One  approach  to  the  problem  to  determining  the  relative 
seismicity  of  various  regions  is  to  plot  earthquake  magnitude 
versus  frequency  of  occurrence  for  the  areas  considered.  The^ 
magnitudes  of  earthquakes  are  known  for  shocks  of  magnitude  six 
and  greater  since  about  1900  but  for  small  shocks  the  records 
are  far  less  complete.  In  contrast,  M.M.  intensities  have  been 
assigned  by  the  Coast  and  Geodetic  Survey  and  other  investigators 
to  nearly  all  significant  earthquakes  known  to  have  occurred 
within  the  United  States.  In  order  to  rrake  use  of  the  complete 
historical  earthquake  record,  maximum  M.M.  intensity  was  plotted 
versus  earthquake  frequency.  A  number  of  assumptions  are  necessary 
when  the  maximum  intensity  is  used  in  a  study  of  earthquake 
frequency.  The  principal  assumptions  are  that  ail  of  the  earth¬ 
quakes  in  the  area  considered  occurred  at  approximately  the  same 
focal  depth  and  that  the  maximum  intensity  for  each  earthquake 
has  been,  in  fact,  reported.  Earthquakes  in  the  United  States- 
are  nearly  all  of  normal  focal  depth,  although  California  and 
Nevada  shocks  are  somewhat  shallower  (averaging  about  TO-15  km/ 
than  in  the  remainder  of  the  country  (about  30-60  km).  It  is 
possible  to  select  geographic  areas  for  study  such  that  the 
earthquakes  within  each  area  are  of  approximately  the  same  depth. 

The  adequacy  of  earthquake  reporting  is  also  a  significant  variable 
and  this  point  will  be  discussed  in  greater  detail. 

In  1860,  all  states  in  the  western  portion  of  the  country, 
with  the  exception  of  California,  had  an  average  population 
density  of  less  than  one  person  per  square  mile.  states  had 

population  densities  greater  than  one  per  square  mile  by  I9l0  with 
the  exception  of  Nevada.  it  is  clear  that  the  low  population 
density  in  the  western  United  States -prior  to  1900  adversely  affects 
the  completeness  of  intensity  data  in  that  area,  particularly  for 
intensities  less  than  VI  or  VII. 
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Reasonably  complete  records  of  great  earthquakes  wi 1 1 
obviously  be  available  f or  a  longer  span  of  years  than  for 
earthquakes  of  small  or  moderate  maximum  intensity.  Accord¬ 
ingly,  the  intensity  data  were  examined  and  an  estimate  was  made 
of  the  total  span  of  years  for  which  maximum  intensity  data  at 
each  level  of  the  intensity  scale  are  nearly  complete.  To 
illustrate  the  technique  further,  Figure  3  shows  the  number  of 
maximum  intensity  VI  earthquakes  reported  for  a  north-south 
trending  zone  of  seismicity  through  western  Montana,  eastern 
Idaho,  central  Utah  and  Arizona  (region  4  of  Figure  5).  The 
data  seem  relatively  complete  from  1900  to  present.'  To  obtain 
the  frequency  of  intensity  VI  earthquakes  in  this  particular 
zone,  the  number  of  intensity  VI  shocks  was  simply  divided  by  the 
time  span,  1900-1967,  or  68  years.  Similarly,  for  all  areas 
considered,  the  average  number  of  earthquakes  at  each  intensity 
level  was  determined  and  plotted  against  intensity. 

,•  The  recurrence  relationship  for  all  intensity  data  in  the 
United  States  is  shown  in  Figure  4.  It  was  assumed,  in  plotting 
the  data,  that  a  relationship  of  the  form  log  N  (per  year)  =  a  +  b  1 
exists,  where  N  is  the  number  of  earthquakes  of  intensity  X* 

Two  conclusions  can  be  drawn  for  this  plot:  1.  maximum  in¬ 
tensities  of  V  and  below  in  general  have  not  been  completely  ' 
reported  for  United  States  earthquakes  and  2.  the  occurrence 
of  great  earthquakes  (M.M.  intensity  XI  and  XII)  is  erratic 
over  the  time  span  considered.  The  recurrence  curve  for  all 
United  States  earthquakes  obviously  combines  data  from  areas 
of  different  geologic  structure  and  tectonic  activity.  Accord- 
ingly,  the  country  was. divided  into  nine  areas  (Figure  5)  in 
an  attempt  to  obtain  recurrence  curves  for  areas  that  seem  to 
be  related  to  more  or  less  definite  trends  or  regions  of  seismic 
activity  and  are  somewhat  geologically  similar  throughout..  . 

Figure  6  shows  the  recurrence  intensities  for  California, 
based  on  approximately  16,000  earthquakes.  The  slope  of  the  least 
squares  line,  log  N  =  a  -f  b  I,  is  0.54  and  corresponds  to  a  value 
of  b  of  0.81  in  the  equation  log  N  =  a  f  b  M  where  M  is  the 
magnitude  and  the  conversion  M  =  1  ■}•  2/3  X  is  used.  For  the 
equation  log  N  =  a  +  b  (8-M),  Gutenberg  and  Richter  (5)  obtained 
a  value  of  0.88  for  Southern  California  and  Allen,  ^  §JL»  (3)  ■  - 

obtained  a  value  of  0.86.  Both  of  these  slopes  are  in  close 
agreement  with  the  value  obtained  from  the  data  in  Figure  6. 

The  curve  for  Nevada  shows  a  similar  slope. 

Area  8  (Figure  5)  includes  the  southeast  Missouri  area 
where  great  earthquakes  occurred  during  1811-12  together  with 
the  St.  Lawrence  Valley  region.  The  recurrence  relationship 
is  shown  in  Figure  7.  The  1811-12  New  Madrid  earthquakes  (maximum 
intensity  =  XI  I )  were  not  used  in  computing  the  least  squares  line. 
There  is  considerable  question  as  to  the  validity  of  considering 
these  two  active  seismic  regions  as  one  area  for  statistical  purposes 
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Recurrence  curves  were  constructed  for  the  regions  separately, 
but  there  were  only  slight  differences  in  the  slopes  of  the 
lines. 

A  summary  of  recurrence  relationships  are  shown  in  Table  1. 

The  areas  are  arranged  in  order  of  decreasing  V1 1  M.M.  activity. 

The  seismic  risk  zones  associated  approximately  with  the  areas 
used  in  the  frequency  studies  are  also  indicated.  Table  11 
shows  the  frequency  of  occurrence  of  earthquakes  of  various 
maximum  intensity  levels  per  lOO  years  per  100,000  km  ^  The 
frequencies  in  Table  1  have  been  divided  by  the  area  of  each 
region.  Assuming  that  the  recurrence  relat ionshi p ' has  some^ 
validity,  these  tables  are  measures  of  the  relative  seismicity 
of  each  area  and  could  be  used  on  a  broad  scale  for  planning 
purposes.  The  results  shov^n  must  be  used  with  many  qualifications 
and  reservations.  The  numbers  in  parentheses  are  recurrence 
rates  based  on  the  least  squares  lines  only,  that  is,  earthquakes 
with  the  indicated  maximum  intensity  are  not  known  to  have  occurred 
/in  that  particular  area.  It  should  be  stated  again  tfmt  the 
above  results  depend  upon  the  assumption  that  a  relationship  of 
the  type  log  N  ~  a+b  I  is  val id  for  these  data  and  that  the 
decisions  made  regarding  the  completeness  of  the  data  sample 
were  reasonably  correct. 
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Saismic  Risk  Map 

A  Seismic  Probability  Map  of  the  United  States  was  prepared 
in  1948  by  F.  P.  Ulrich  with  the  advice  of  seismologists  through¬ 
out  the  United  States  (6)  and  was  issued  by  the  Coast  and  Geodetic 
Survey  in  1948.  In  1949,  the  Seismic  Probability  Map  was  revised 
such  that  the  Charleston,  South  Carolina  area  was  changed  from 
Zone  3  to  Zone  2  and  a  Zone  3  was  set  up  for  the  Puget  Sound 
region  of  Washington  which  had  formerly  been  included  in  Zone  2 
-  (7).  The  revised  map  was  adopted  by  the  Pacific  Coast  Buildi'ng 
Officials  Conference  for  inclusion  in  the  1952  edition  of  the 
Uniform  Building  Code.  Subsequent  editions  of  the  Uniform 
Building  Code  have  included  this  nap  with  no  changes  (8). 

The  map  was  withdrawn  by  the  Coast  and  Geodetic  Survey  in 
1952  as  "subject  to  misinterpretation  and  too  general  to  satisfy 
the  requirements  of  many  users."  The  annual  publication  Unitsd 
States  Earthquakes  has  included  a  map  showing  locations  of  ^ 

■  destructive  and  near  destructive  earthquakes  throughout  the 
United  States,  but  the  Coast  and  Geodetic  Survey  has  offered 
no  other  probability  nap. 

The  Seismic  Probability  Map  has  been  most  useful  as  part 
of  the  criteria  in  the  establishment  of  lateral  force  requirements 
for  buildings  but  it  does,  however,  have  several  obvious  limitations* 
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As  pointed  out  by  Richter  (9),  the  circular  spot  numbered  3  in 
Montana  corresponds  to  damagine  earthquakes  in  1935,  but  the  map 
i snores  other  earthquakes  and  known  active  structures  in  the  same 
seolopical  province.  Along  many  zone  boundaries,  Zone  3  borders 
Zone  1  and  Lne  2  borders  Zone  0.  If  a  major  earthquake  is 
considered  possible  in  Zone  3,  the  nor^l  distribution  of  in- 
tensities  would  suggest  that  an  area  of  Zone  2  should  surround  it, 
the  same  reasoning  being  true  for  the  peripheries  of  Zone  2 
areas.  The  original  Seismic  Probabi 1 i ty  Map  describes  U 

through  3  only  in  terms  of  ”no  damge,  minor  (^mage,  moderate 
damage  and  major  damage”  respectively.  No  attempt  was  made  to 
correlate  the  zones  with  intensity. 

The  rrap  presented  in  Figure  8,  which  will  be  termed  a 
"Seismic  Risk  Map,”  is  meant  to  be  an  interim  revision  of  ' 
original  Coast  and  Geodetic  Survey  Seismic  Probabi 1 ity  Map. 

The  map  in  Figure  8  is  actually  more  than  a  revision  of  the 
old  map  in  that  the  factors  used  in  constructing  the  new  risk 
rap  are  somewhat  different  from  those  considered  in  the  original 

map. 

The  Seismic  Risk  Map  presented  here  is  based  on  the  following 
factors! 

1.  the  distribution  of  M.M.  intensities  associated 
with  the  known  seismic  history  of  the  United  States, 

2.  strain  release  in  the  United  States  since  1900,  and 

3.  the  association  of  strain  release  patterns  with 
large  scale  geologic  features  bel ieved  to  be 
related  to  recent  seismic  activity. 

Where  seismic  activity  has  occurred  intermittently  along  a 
recognizable  geologic  trend,  it  has  been,  in  general,  assumed 
that^seismic  activity  could  occur  with  equal  1 ikel ihood  anywhere 
along  that  trend  or  structure.  In  areas  where  the  relations^ 
between  seismicity  and  geologic  structure  are  not  clear  or  where 
only  limited  geologic  information  is  available,  the  risk  zones 
are^based  on  the  distribution  of  M.M.  intensitip  and  strain 
release.  Inall  cases,  the  size  and  shape  of  the  zones  have  been 
heavily  unfluenced  by  the  historical  distribution  of  intensities 
(Figure  l).  It  is  realized  that  it  would  be  desirable  to  relate 
risk  to  a  measurable  quantity  less  subjective  than  M.M.  intensay, 
but  the  available  strong  motion  and  other  data  appear  to  be  in¬ 
sufficient  to  attempt  this  at  present.  No  special  corrections 
are  presented  for  types  of  surficial  geology  and  soils  ttet 
may  increase  or  decrease  the  intensity  of  shaking.  It  should  be 
rerrerked,  however,  that  the  M.M.  intensities  reported  for  an 
earthquake  are  usually  the  maximum  effects  observed  on  the 
worst  ground. 
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The  risk  zones  in  the  northeast  portion  of  the  United  States 
are  not  greatly  changed  from  the  original  risk  -map.  The  St. 
Lawrence  Valley  continues  to  be  Zone  3  because  of  the  great 
earthquakes  in  1663  and  1925  and  lesser  shocks  near  Attica,  N.  Y. 
in  1929  and  near  Massena,  N.  Y.  in  1944.  The  Boston,  Massachusetts 
area  was  zoned  3  in  consequence  of  the  earthquake  east  of  Cape  Ann 
in  1755  and  many  other  lesser  earthquakes  in  the  area.  The  Zone 
3  rating  in  the  Charleston,  South  Carolina  area  is  a  consequence 
of  the  1886  earthquake.  Zone  2  has  been  extended  to  include  the 
main  structures  of  the  Appalachian  area  from  Pennsylvania  to 
Alabama  which  are  associated  with  numerous  moderate  earthquakes. 
Zone  3  and  Zone  2  Tn  the  Mississippi  Valley  have  been  very  much 
enlarged  over  those  zones  on  the  old  Seismic  Probability  Map. 

The  new  zone  more  nearly  reflects  the  probable  distribution  of 
intensities  in  any  repetition  of  the  I8ll~l2  series  near  New 
Madrid,  Missouri.  Zone  2  has  been  extended  northeast  across 
a  zone  of  relatively  minor  seismicity  in  Indiana  to  southwestern 
Ohio.  Although  earthquakes  with  maximum  M.M.  intensities  of 
V11  are  known  in  western  Ohio  (lO),  connection  of  Zone  2  in 
Ohio  with  the  Zone  2  surrounding  the  New  Madrid,  Missouri 
area  is  a  matter  of  interpretation.  This  zone  is  based  on 
the  apparent  alignment  of  epicenters  from  the  St.  Lawrence 
Valley  to  southeast  Missouri.  There  is  no  conclusive  geologic 
evidence  at  present  for  the  alignment,  and  other  interpretations 
are  possible  (lO). 

A  Zone  2  area  has  been  drawn  along  the  Nenraha  Ridge  structure 
in  Nebraska,  Kansas  and  Oklahoma  which  includes  the  1952  earth¬ 
quake  near  El  Reno,  Oklahoma  and  numerous  other  shocks  along  the 
Nemaha  Ridge.  The  north-central  plains  states  have  been,  for 
the  most  part,  rated  in  Zone  1  because  of  the  lack  of 'data  relating 
seismicity  with  geologic  structure  and  because  of  the  rather 
widespread  and  erratic  distribution  of  moderate  earthquake  activity. 
North  Dakota  is  a  good  example.  No  earthquakes  (11)  are  known  , 
to  have  occurred  in  the  state  until  July  8.  1968  when  a  magnitude 
4.4  earthquake  shook  the  Bismarck  area  (12;.  The  Zone  3  area 
restricted  to  southeastern  Montana  on  the  old  map  has  been  extended 
southward  through  Idaho  and  Utah.  The  large  earthquake  near  Hebgen 
Lake,  Montana  in  1959,  the  Kosmo,  Utah  shock  of  1934,  the  series  of 
shocks  in  1920  and  1921  near  Elsinore,  in  southern  Utah  and  the 
related  north-south  trending  faults  amply  justify  the  zoning. 

Only  minor  changes  have  been  made  from  the  original  map  in 
the  states  bordering  the  Pacific.  The  original  Zone  1  rating  in 
Oregon  has  been  slightly  changed.  Other  minor  changes  have 
been  made  in  the  shape  of  the  Zone  3  area  in  Washington  at  the 
south  end  of  Zone  2  in  the  Great  Valley  of  California. 
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n|o:rM^.ciion  and  Summai:/ 

The  Seismic  Risk  Map  presented  is  offered  a  revision 
of  the  Seismic  Probability  Map  originally  prepared  by  the  Coast 
and  GeodeUc  Survey  in  1947.  The  new  ^p  is  only  an  interim  one 
and  doL  not  repreLnt  the  final  form  of  a  risk  map  for  the 
United  States.  It  does  not  consider  the  frequency  of  occurrence 
of  the  seismic  events  in  each  zone.  An  estimate  of  the 
rates  of  occurrence  of  earthquakes  of  various  maximum  M.M. 
intensity  in  different  zones  has  been  made  in  this  paper 
using  strain  release  and  earthquake  recurrence  data. 

It  is  suggested  that  eventually  it  may  be  possible  to 
s.;+imate  seismic  risk  in  the  United  States  using  three  separate 

naps:  one  map  would  be  geological 

citructure  thought  to  be  related  to  seismic  activity,  with  assign 
menrof  relative  risk  to  each  of  the  structures;  a  second  map 
would  estimate  the  frequency  of  occurrence  of  ^^['^'"Jijakes  of 
various  magnitudes  or  maximum  intensities  throughout  . 

States;  and  a  third  map  would  show  the  distribution  o 
M  M  intensity  in  the  United  States  such  as  is  shown  in  Figure  1 
oAhis  rper!  If  proper  coefficients  are  assiped  to  each  area 
or  featuK  of  each  of  the  three  maps,  it  would  be  possible  to 
obtain  both  the  long  term  (risk  over  a  very  long  time!  and  short 
term  (rUk  over  25-50  years)  risks  in  the  fol lop ng  manner. 

The  short  term  risk  at  a  particular  site  could  be  obtaped  by 
a  weighted  consideration  of  coefficients  from  maps  1,  2  and  3. 
Long  term  risk  could  be  obtained  from  a  similar  integrati 
of  Lta  from  maps  1  and  3.  Much  additional  wpk  remains  to 
be  tone,  particularly  with  regard  to  the  relationship  betpen 
oeolopic  Sructure  and  earthquake  occurrence.  It  is  hoped  that 
the  data  and  maps  presented  will  be  useful  in  developing  u  ure 
seismic  risk  maps  for  the  United  States. 
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SIMMARV  nr  FARTHOIIAKE  RECURRENCE  FORMULAS 

LOGN=a+bl  ■  EARTHQUAKES  PER  100  YEA FB 

AREA  -RISK  ZONE  (per  year) 
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SUMMARY  OF  EARTHQUAKE  RECURRENCE 

FORMULAS 
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AREA 

RISK  ZONE 
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100  YEARS  PER 

V11  • 
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Figure  l.^  Maximum  Modified  Mercall<i'  intensities  throughoui  the  United 
'•  V  States  -  ail  historical  data  through  1965 
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Figure  5  Location  map  showing  the  areas  for  which  recurrence  formulas  were 
computed 
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SEISMIC  RISK  MAP  OF  THE  UNITED  STATES 

ZONE  0- No  damage. 


ZONE  1  -  Minor  damage;  distant  earthquakes  may  cause  damage 
to  structures  with  fundamental  periods  greater  than  ' 
1.0  seconds:  corresponds  to  intensities  V  and  VI 
of  the  M.M.*  Scale. 


ZONE  2  -  Moderate  damage;  corresponds  to  intensity  VII  cf  the  M.M.*  Scale. 

I  ZONE  3  -  Major  damage;  corresponds  to  intensity  VIII  and  higher  of  the  M.M.*  Scale. 


This  map  is  based  on  the  known  distribution  of  damaging  earthquakes  and  the 
M.M.*  intensities  associated  with  these  earthquakes;  evidence  of  strain  release; 
and  ponsidcation  of  major  geologic  structures  and  provinces  believed  to  be 
associated  with  earthquake  activity.  The  probable  frequency  of  occurrence  cf 
damagirig  earthquakes  in  each  zone  was  rot  considered  in  assigning  ratings  to 
the  various  zones.  See  accompanying  text  for  discussion  of  trequency  of 
earthquake  occurrence. 


•Modified  Mercal'i  Intensity  Scale  of  1931. 
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Seismicity  of  Colorado:  Consistency  of  Recent 

Earthquakes  with  Those  of  Historical  Record 

« 

Abstract  Earthquakes  instrumentally  recorded  from  1966  to  1968  have  oc¬ 
curred  in  the  same  regions  of  western  Colorado,  the  Arkansas  and  Platte  river 
valleys,  as  those  felt  back  to  1870  (from  newspaper  reports),  despite  the 
increasing  cultural  effects  of  mining,  highway  construction,  reservoir  building, 
and  loading.  Thus  it  appears  unnecessary  to  explain'  the  Denver  earthquakes  in 
terms  of  pressure  induced  by  the  introduction  of  waste  fluid.  However,  the 
assumption  of  preexistent  tectonic  strains  in  the  area  of  the  Rocky  Mountain 


Arsenal  seems  to  be  confirmed. 

Colorado  is  considered  a  region  of 
minor  seismicity.  However,  the  advent 
of  hundreds  of  earthquakes  in  the  area 
of  northeastern  Denver  since  1962, 
when  the  Rocky  Mountain  Arsenal  be¬ 
gan  pumping  fluids  under  pressure  into 
a  3600-m  disposal  well,  caused  concern 
that  a  more  destructive  earthquake 
might  possibly  occur  there  in  the  fu¬ 
ture.  This  study  does  not  consider  a 
sufficiently  long  interval  of  time  to  yield 


seismic  information  on  the  distant  past 
before  the  Rocky  Mountain  region  be¬ 
came  populated.  The  data  collected  do 
show  a  picture  of  seismicity  for  1966 
through  1968  consistent  with  that  de¬ 
rived  from  reports  of  felt  shocks  in 
newspaj>ers  published  100  years  ago. 
Figure  1  shows  the  known  locations 
of  earthquakes  in  Colorado  from  1870 
to  1967.  This  map,  adapted  from  Had- 
sell’s  study  (/),  indicates  that  the  ma¬ 


jority  of  Colorado  earthquakes  in  this 
period  have  occurred  west  of  the  Rocky 
Mountain  Front  Range  and  within  30 
km  of  post-Oligoccnc  extrusive  rocks. 
Reports, of  felt  shocks  from  25  news¬ 
papers,  published  in  1882  and  later, 
were  used  by  Hadsell  to  estimate  the 
intensities  listed  in  Table  1. 

A  detailed  study  of  seismograms 
from  the  Colorado  School  of  Mines 
station  (GOL)  (2)  with  additional  data 
from  the  Uinta  Basin  Observatory 
(UBO)  at  Vernal,  Utah  (3),  produced 
sufficiently  accurate  information  on  the 
epicenter  locations  to  make  possible 
reports  of  events  recorded  in  the  years 
1966,  1967,  and  1968.  Figure  2  shows 
locations  of  seismic  events  for  these 
years.  The  trend  of  epicenters  for  all 
3  years  is  north-south,  parallel  with 
the  Continental  Divide,  with  most  epi¬ 
centers  placed  west  of  the  divide.  East¬ 
ern  Colorado  is  nearly  ascismic,  with  a 
few  epicenters  in  the  Arkansas  and 
Platte  river  valleys.  There  are  con¬ 
spicuous  areas  of  the  state  where  no 
earthquakes  occurred,  either  historical¬ 
ly  or  in  the  3-year  study.  These  inactive 
regions,  in  addition  to  the  eastern 
plain,  are  in  the  northwestern  and 
southwestern  corners  of  Colorado.  The 
Sangre  de  Cristo  Mountains  in  the 
south-central  part  of  the  state  also  re¬ 
main  inactive.  Table  1  lists  the  annual 
number  of  events  in  each  epicentral 
region. 

Nine  of  these  epicentral  regions  have 
produced  earthquakes  of  Richter  mag¬ 
nitude  greater  than  4.0.  Only  one  of 
these  locations  had  shocks  of  large 
magnitude  in  both  1966  and  1967.  This 
was  at  Dulce,  New  Mexico,  on  the 
southern  border  of  central  Colorado; 
these  shocks  were  widely  felt  in  south¬ 
ern  Colorado.  Other  epicenters  in  1966 
having  earthquakes  larger  than  magni¬ 
tude  4.0  were  Grand  Valley,  Leadville, 
Rangely,  and  Trinidad.  In  1967  north¬ 
eastern  Denver  (Derby)  earthquakes 
were  the  largest  in  the  state,  with  three 
shocks  above  magnitude  5,  apd  the 
greatest  number,  314  events,  recorded 
•  at  GOL  (4).  Other  epicentral  locations 
in  1967  experiencing  earthquakes  great¬ 
er  than  magnitude  4.0  were  Dulce, 
Montrose,  Silverton,  and  Pueblo.  In 
1968  there  were  no  earthquakes  with 
magnitude  >  4.0  recorded  at  GOL. 
Seismicity  of  Colorado  for  1968  was 
less  than  for  the  previous  2  years. 

Earthquakes  reported  in  this  study 
were  of  magnitude  ^  1.0;  most  were 
greater  than  magnitude  2.0,  with  the 
exception  of  the  Denver  events.  Den- 
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Table  1.  Summary  of  earthquakes  in  Colorado  for  the  years  1966,  1967,  and  1968  com¬ 


pared  to  historical  record. 


Nos. 

(Fig. 

2> 

Location 

1966 

1967 

1968 

Historical  records! 

Nos. 

(Fig. 

1) 

1 

Aspen 

10 

12 

13 

1880  VI,  1889  V  (Glenwood 

Springs),  1941  IV,  1960  V 

1 

2 

Baxter  Mountain 

4 

20 

6 

3 

Cabin  Creek 

19 

73 

17 

1871  IV  (Georgetown),  1894  V 

2 

4 

Castle  Rock 

1 

0 

0 

5 

Climax 

12 

31 

0 

6 

Colorado  Springs 

0 

0 

3 

7 

Denver  (Derby) 

189 

314* 

99 

1-882  VII,  1916  IV  (Boulder), 

1962-65  VII 

3 

8 

Dulce, 

New  Mexico 

161* 

31* 

14 

1952  V  (Conejos  County) 

4 

9 

Eikhead  Mountains 

152 

177 

85 

10 

Four  Corners 

0 

3 

5 

1941  V  (Durango) 

5 

11 

Golden  Fault 

0 

31 

0 

12 

Grand  Valley 

2* 

6 

2 

1954  IV  (Grand  Junction) 

6 

13 

Harvey  Gap 

Reservoir 

0 

1 

0 

14 

Lake  Hattie 

Reservoir 

16 

4 

1 

15 

Leadville 

28* 

13 

4 

1901  VI  (Buena  Vista), 

1961  IV 

7 

16 

Montrose 

12 

22* 

11 

1960  VI,  1962  IV  (Cimarron) 

8 

17 

Mount  Evans 

0 

5 

0 

18 

Mount  Gunnison 

47 

67 

29 

1882  IV,  1944  VI 

9 

19 

Mount  Wilson 

8 

12 

5 

1913  V  (Ouray) 

10 

20 

North  of  Eikhead 

11 

7 

17 

21 

Pueblo 

5 

12* 

10 

1870  VI,  1888  IV  (San 

Isabel  National  Forest), 

11 

1955  IV  (Rocky  Ford), 

1956  V,  1963  IV 

12 

22 

Rangely  . 

6* 

25 

4 

1891  VI  (Axial  Basin) 

13 

23 

Ruth  Mountain 

16 

51 

13 

24 

Sheep  Mountain 

1 

0 

0 

1955  VI  (Lake  City),  1928  V 

(Creede) 

14 

25 

Silverton 

2 

1* 

2 

1882  IV 

15 

26 

Snowmass  Mountain  1 

5 

8 

27 

Steamboat  Springs 

5 

20 

59 

1895  V,  1955  V 

16 

28 

Tenderfoot 

Mountain 

0 

16 

3 

29 

Trinidad 

5* 

2 

1 

30 

Urad 

15 

15 

71 

31 

Wi  ite  River 

1 

6 

4 

32 

Yampa 

0 

7 

16 

Totals 

729 

989 

502 

•Magnitude  >4.0.  t  Roman  numerals  refer  to  estimated  Intensities  on  the  modified  Mercalli 
scale  UO). 


AS-26 


Fig.  t.  Historical  seismicity  of  Colorado 
adapted  from  studies  of  F.  A.  Hadsell  (7). 
Circles  represent  locations  of  epicenters. 
Numbers  on  map  refer  to  numbers  in  the 
last  column  of  Table  1.  Large,  irregular,  . 
solid-black  regions  represent  post-Oligo- 
cene  (less  than  26  x  10®  years)  extrusive 
volcanic  rocks, 

ver  is  45  km  from  the  Bergen  Park 
Observatory;  thus  events  of  magnitude 
1,0  could  be  easily  interpreted,  even  in 
winter  months  when  the  magnification 
of  the  instruments  is  lowered  to  reduce 
background  noise. 

Possibly  20  percent  of  the  total 
number  of  events  recorded  are  due  to 
mining  operations.  At  Urad,  for  in¬ 
stance,  reports  were  regularly  received 
of  charge  sizes  and  approximate  times 
of  blasting.  In  1968  this  operation 
reached  its  peak,  with  the  planned 
cave-in  of  a  mountain.  Fifty  natural 
events  at  Urad  occurred  within  a  month 
of  this  collapse,  half  of  this  number 
within  a  few  days.  At  other  locations 
it  is  not  always  possible  to  distinguish 
between  natural  earthquakes  and  ex¬ 
plosive  sources.  Several  epicentral  loca¬ 
tions  where  blasting  operations  are 
known  to  occur  are  Elkhead  Moun¬ 
tains,  Climax,  and  Mount  Gunnison. 
Usually  mining  is  carried  on  only  dur¬ 
ing  daylight  hours,  so  that  events  oc¬ 
curring  between  the  hours  of  2400  and 
1200  G.M.T.  are  probably  natural 


Fig.  2.  Map  represents  locations  of  earth¬ 
quakes  during  1966,  1967,  and  1968  in¬ 
strumental  ly  recorded  in  Colorado  at 
GOL  and  in  Utah  at  UBO.  Circles  rep¬ 
resent  areas  approximately  15  km  in 
diameter.  Numbers  within  each  circle  re¬ 
fer  to  numbers  in  the  first  column  of 
Table  1. 

events.  Blasting  at  Climax  has  been 
carried  on  here  around  the  clock.  No 
large  blasts  were  recorded  from  Climax 
in  1968.  Highway  blasting  has  been 
omitted  from  the  data  when  it  is  known 
to  have  occurred.  Construction  com¬ 
panies  notified  us  in  advance  of  large 
blasts. 

Seismic  activity  has  been  observed 
at  some  reservoir  sites,  namely  Cabin 
Creek,  Lake  Hattie,  Harvey  Gap,  and 
Mount  Gunnison  (Blue  Mesa  Reser¬ 
voir).  Carder  described  earthquakes 
occurring  along  the  southern  edge  of 
Lake  Mead,  Nevada,  with  water-load¬ 
ing  of  the  reservoir  in  1932  (5).  On  the 
other  hand,  no  earthquakes  have  oc¬ 
curred  at  sites  in  Colorado  where  the 
U.S.  Bureau  of  Reclamation  has  been 
constructing  reservoirs,  several  of  which 
are  already  filled.  Nine  of  these  sites 
are  within  a  radius  of  150  km  of  the 
Bergen  Park  Observatory,  a  distance 
where  even  small  events  would  be  well 
recorded.  No  events  were  observed 
from  these  nine  sites  in  the  3  years 
studied.  Earthquakes  did  occur  in  the 


Table  2.  Number  of  Denver  earthquakes  per  month,  magnitude  >  1.0,  recorded  at  Bergen 
Park  Observatory  (GOL)  for  the  years  1962  through  1968,  ~ 


Month 

1962 

1963 

1964 

1965 

1966 

1967 

1968 

January 

0 

24 

5 

4 

25 

4 

10 

February 

0 

35 

2 

32 

11* 

19 

8 

March 

0 

38 

8 

12 

20 

2 

3 

April 

2 

72 

9 

27 

11 

72 

9 

May 

30 

26 

3 

25 

18 

31 

11 

June 

45 

14 

4 

68 

12 

14 

10 

July 

24 

21 

5 

132 

18 

11 

14 

August 

30 

11 

9 

65 

11 

74 

9 

September 

23 

14 

4 

105 

22 

11 

10 

October 

9 

20 

13 

23 

9 

16 

6 

November 

8 

4 

3 

43 

23 

44 

7 

December 

20 

5 

8 

14 

9 

9 

2 

Totals 

191 

284 

73 

550 

189 

314 

99 

*  Pumping  stopped  at  the  Rocky  Mountain  Arsenal  disposal  well. 
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historical  record  at  locations  where 
reservoirs  have  since  been  constructed, 
especially  at  Mount  Gunnison  and 
Harvey  Gap,  in  western  Colorado,  and 
Cabin  Creek  (Georgetown)  in  central 
Colorado.  The  occurrence  of  earth¬ 
quakes  at  reservoir  sites  would  seem  to 
depend  not  on  the  size  of  the  body  of 
water  but  rather  on  the  geologic  struc¬ 
ture  beneath  it. 

Recent  earthquakes  began  at  Derby 
in  April  1962.  This  was  shortly  after 
the  U.S.  Army  began  pumping  waste 
Iluids  under  pressure  into  the  deep 
disposal  well  at  the  Rocky  Mountain 
Arsenal  in  northeastern  Denver.  Evans, 
a  Denver-based  consulting  geologist,  hy¬ 
pothesized  a  connection  between  the 
Denver  earthquakes  and  this  well  (6). 
Healy  et  a/,  theorized  that  the  Denver 
earthquakes  are  man-made  (7). 

The  arsenal  well  was  constructed 
through  3650  m  of  sedimentary  layers 
into  fractured  Precambrian  basement 
rock.  Studies  were  carried  out  at  the 
Colorado  School  of.  Mines  with  a  dense 
network  of  instruments  to  determine 
whether  there  was  any  relation  between 
the  Denver  earthquakes  and  the  arsenal 
well  (4).  Pumping  of  wastes  into  the 
well  was  stopped  hi  February  1966. 
The  earthquakes  continued.  The  num¬ 
bers  of  earthquakes  recorded  at  Bergen 
Park  Observatory  for  the  years  1962 
through  1968,  with  monthly  and  yearly 
totals,  are  listed  in  Table  2, 

The  largest  earthquakes  in  the  Den¬ 
ver  area  occurred  in  1967,  a  year  after 
the  cessation  of  pumping  of  waste 
fluids  into  the  Rocky  Mountain  Arsenal 
well.  The  relation  of  frequency  of  oc¬ 
currence  to  magnitude  of  the  Denver 
earthquakes  is  similar  to  that  observed 
in  other  geographic  regions  where 
earthquake  activity  is  not  influenced 
by .  man.  However,  in  1 967,  the  ex¬ 
pected  great  number  of  small  shocks 
that  usually  accompany  a  few  large 
ones  did  not  occur. 

As  debate  about  the  causal  relation 
between  earthquakes  in  Denver  and  the 
arsenal  well  continued.  Army  Engineers 
began  a  series  of  pumping  tests  in  Sep¬ 
tember  and  October  1968  to  determine 
whether  removal  of  small  amounts  of 
fluid  from  the  well  would  change  the 
earthquake  occurrence.  A  report  of 
their  results  indicates  that  many  small 
shocks  (as  low  as  magnitude  —1.0) 
did  occur  near  the  well  during  the  pe¬ 
riod  studied  (S).  Accof;^i:^ly,  the  tests 
were  discontinued.  Tti^riWniber  of 
earthquakes  in  the  Dea.ver  area,  re¬ 
corded  at  GOL  in  1968,  was  99.  The 
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largest  magnitude  in  1968  was  3.2,  re¬ 
corded  in  November.  The  total  number 
of  Denver  earthquakes  recorded  since 
the  Bergen  Park  Observatory  began 
operations  in  1962  is  1698  of  magni¬ 
tude  55  1.0. 

Rangcly,  Colorado,  a  seismically  ac¬ 
tive  region  during  the  years  studied,  is 
near  the  southeastern  edge  of  a  large 
oil  field.  An  earthquake  of  magnitude 
4.5,  located  by  the  U.S.  Coast  and  Geo¬ 
detic  Survey,  occurred  there  iT\  1966. 
Munson  tried  to  correlate  numbers  of 
earthquakes  near  Rangely  with  millions 
of  gallons  of  water  pumped  into  the 
field  as  oil  was  withdrawn  (9).  He 
found  no  change  in  numbers  of  earth¬ 
quakes  whether  pumping  rate  increased 
or  decreased  in  the  period  of  time  con¬ 
sidered. 

Colorado  seismicity  has  been  virtu¬ 
ally  unchanged  in  the  last  100  years. 
The  same  areas  that  were  seismically 
active  in  the  historical  record  were 
seismically  active  in  the  years  1966, 
1967,  and  1968;  those  areas  that  were 
seismically  inactive  100  years  ago  were 
inactive  in  1966,  1967,  and  1968.  The 
intervening  years  have  seen  a  large 
population  growth  in  the  Rocky  Moun¬ 
tain  region,  along  with  increased  inter¬ 
est  in  seismic  activity  and  the  earth¬ 
quake  hazard  problem  for  construc¬ 
tion  projects. 

The  possibility  of  a  devastating  earth¬ 


quake  of  greater  than  magnitude  6, 
which  might  cause  loss  of  life  and  ex¬ 
tensive  damage  to  property,  always  ex¬ 
ists.  On  the  basis  of  this  seismicity 
study,  however,  the  possibility  that  Col¬ 
orado  will  experience  a  larger  magni¬ 
tude  shock  than  has  yet  occurred,  and 
its  location  and  time,  cannot  be  pre¬ 
dicted. 

Ruth  B.  Simon 

Department  of  Geophysics,  Colorado 
School  of  Mines,  Golden  80401 

References  and  Notes 

1.  F.  A.  Hadsell,  Colo.  Sch.  Mines  Quart,  6J 
(1),  57  (1968). 

2.  Tlic  Cecil  H.  Green  Geophysical  Observatory 
of  the  Colorado  School  of  Mines  (GOL), 
located  at  Bergen  Park,  Colorado,  is  one  of 
the  World  Wide  Standard  Seismograph  Sta- 

•  tions  network  supported  by  the  Environmental 
Science  Services  Administration  and  the  U.S. 
Coast  and  Geodetic  Survey. 

3.  The  Uinta  Basin  Seismological  Observatory 
(UBO).  located  at  Vernal,  Utah,  is  operated 
by  Geotech  Corporation  for  Advanced  Re¬ 
search  Projects  Agency.  I  thank  these  agendea 
for  data  provided  for  this  study. 

4.  M.  W.  Major  and  R.  B.  Simon,  Colo.  Sch. 
Mines  Quart.  63(1),  9  (1968). 

5.  D.  S.  Carder,  paper  given  at  the  meeting  of 
the  Geological  Society  of  America,  Rocky 
Mountain  Section,  10-14  May  1967. 

6.  D.  M.  Evans,  Mountain  Ceol.  3,  No.  1  (1965). 

7.  J.  H.  Hcaly,  W.  W.  Rubey,  D.  T.  Griggs, 
C.  B.  Raleigh,  Science  161.  1301  (1968). 

8.  D.  B.  Hoover  and  J.  A.  Dietrich,  U.S.  Geol. 
Surv.  Circ.  613  (1969). 

9.  R.  C.  Munson,  thesis,  Colorado  School  of 
Mines  (1968). 

10.  C.  F.  Richter,  Elementary  Seismology  (Free¬ 
man,  San  Francisco,  1958),  p.  137, 

11.  Supported  by  the  State  of  Colorado,  the 
Colorado  School  of  Minos,  and  the  U.S. 
Army  Corps  of  Engineers. 

22  April  1969;  revised  17  June  1969 


bbosp 


M 


EARTHQUAKE  NOTES 

« 

Vol.  XLIII  June  1972  Number  2 

SEISMICITY  OF  COLORADO 
1969  -  1970  -  1971 

* 

Ruth  B.  Simon 

Colorado  School  of  Mines 
Golden,  Colorado  ' 


The  last  three  years  of  Colorado  seismic  events  closely  resemble  ^ 
those  of  the  preceding  three  years,  from  1966,  when  these  studies  began. 
Epicentral  locations  are  approximately  the  same,  with  one  or  tv/o  new  sites 
each  year,  and  an  absence  of  activity  at  one  or  two  previously  active 
sites,  so  that  the  number  of  locations  remains  about  30  in  the  state. 

In  1971  the  new  Colorado  network  of  short  period  seismograph  stations 
was  set  up  under  joint  auspices  of  NOAA  and  CSM,  These  stations:  Alamosa, 
Craig,  Durango,  Eagle,  Grand  Junction,  Montrose,  and  Trinidad,  have  oper¬ 
ated  intermittently  since  June  1971.  Using  arrival  times  from  at  least 
three  of  these  stations  recording  simultaneously,  48  events  of  Mag  2-4.5 
were  located  with  an  accuracy  of  1-5  km,  in  the  second  half  of  1971.  (See 
map  attached).  The  accuracy  of  locations  using  the  two  standard  stations, 
GOL  seismograms,  and  UBO  bulletins,  has  been  15  km. 


Telemetered  recordings  are  received  at  CSM  and  three  of  the  stations 
are  telemetered  also  at  NOAA  Earthquake  Research  Laboratory  in  Boulder. 
Stringent  financial  problems,  however,  make  future  operation  of  the  net¬ 
work  for  the  rest  of  1972  doubtful.  This  is  regrettable  as  seismic  veloc¬ 
ity  studies  in  the  complex  Rocky  Mountain  region  had  barely  begun.  The 
need  for  carefully  timed,  v/ell-recorded  data  over  a  long  period  of  time 
is  essential. 

Table  1  lists  the  epicentral  locations,  their  distances  and  azimuths 
from  GOL  (the  Cecil  H.  Green  Geophysical  Observatory  at  Bergen  Park)  and 
UBO  (the  Uinta  Basin  Observatory  at  Vernal,  Utah).  This  table  lists  also 
the  numbers  of  events  at  each  location  for  the  years  1969,  1970,  and  1971. 
In  1969  there  was  one  event,  at  Montrose,  with  a  magnitude  of  4.0,  but  no 
events  in  1970  of  this  magnitude.  The  AEC  Plowshare  explosion  at  Rulison 
on  September  10,  1969,  had  a  Mag  5.0,  but  its  aftershocks  were  all  much 
smaller  than  4.0.  In  1971  there  were  seven  locations  with  events  greater 
than  Mag  4.0.  Maps  for  these  years  are  attached. 

Table  2  lists  all  seismic  events  for  three  years  according  to  hour  of 
day  of  occurrence.  It  is  not  possible  to  distinguish  blasting  activities 
from  natural  earthquakes  on  seismograms.  Assuming  that  most  mining  oper¬ 
ations  are  conducted  during  daylight  hours,  there  should  be  an  increase  in 
numbers  of  events  from  8  a.m.,  1400  GCT  to  5  p.m. ,  2300  GCT  ±  1  hour. 
Plotting  Table  2,  Figure  1  shows  that  the  increase  in  activity  begins 
about  three  hours  later  the  working  day,  11  a.m.  MST,  1700  GCT,  and 

decreases  abruptly  about  8  p.m.  MST,  0200  GCT,  for  all  three  years.  This 


Simon,  Ruth  B. ,  Seismicity  of  Colorado;  Consistency  of  Recent  Earth 
quakes  with  Those  of  Historical  Record,  Science,  Aug.  29,  1969,  Vol.  165, 
pp  897-699. 
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SEISMICITY  OF  COLORADO 


Location 

Aspen 
Baxter  Mt. 

Bellyache  Mt. 

Cabin  Creek 
Climax 

Colo.  Sprgs.  Vic. 

Del  Norte 
Derby  (Denver) 

Dulce 

Eaglets  Roost 
Elkhead  Mts. 

Four  Comers 
Glenv/ood  Sprgs.  Vic.** 
Grand  Valley 
Harvey  Gap  Res. 

Lake  Hattie  Res. 

Leadville 

Louisville 

Montrose 

Mt.  Gunnison 

Mt .  Wilson 

Pueblo 

Rangely 

Ruth  Mt. 

San  Antonio  Peak 
Sheep  Mt. 

Silverton 
Snowmass  Mt. 

Steamboat  Sprgs. 
Tenderfoot  Mt. 
Trinidad 
Urad 

Walsenberg 
White  River 
Yampa 


Totals 


TABLE  1 

’&  AZ  GOL 

130  Km  SW 
175  Km  W 
110  Km  W 
35  Km  W 
80  Km  SW 
100  Km  SSE 
235  Km-  SSW 

45  Km  NE 
320  Km  SSW 

55  Km  WIW 
200  Km 
420  Km  SW 
160  Km  W 
325  Km  WSW 
175  Km  W 
180  Km  mVif 
90  Km  SW 

46  Km  NNE 
250  Km  SW 
200  Km  SW 
300  Km  SW 
175  Km  SSE 
320  Km  WI'JW 
200  Km  SW 
390  Km  SSW 
260  Km  SW 
280  Km  SW 
160  Km  SW 
150  Km  NW 
180  Km  SW 
300  Km  SSE 

50  Km  W 
225  Km  SSE 
175  Km  WNW 
135  Km  NW 


No.  of  Events 


& 

AZ 

UBO 

1969 

1970  1971 

266 

Km 

SE 

33 

11 

16 

188 

Km 

E 

17 

25 

21 

250 

Km 

E 

0 

1 

3 

550 

Km 

E 

47 

44 

21 

525 

Km 

ESE 

3 

5 

4 

500 

Km 

ESE 

5^ 

2 

2 

455 

Km 

SE 

0 

0 

1 

400 

Km 

ESE 

74 

56 

35* 

520 

SSE 

68 

56 

33 

300 

Km 

BSE 

15 

11 

2 

235 

Km 

NE 

79 

59 

55 

400 

Km 

S 

16 

9 

2 

220 

Km 

ESE  ' 

0 

0 

1* 

180 

Km 

E 

5 

1 

2* 

190 

Km 

E 

0 

3 

3 

570 

Km 

E 

0 

2 

1 

350 

Km 

ESE 

7 

2 

0 

330 

Km 

E 

0 

0 

1  ' 

250 

Km 

SE 

21* 

9 

10* 

250 

Km 

SE 

15 

13 

20 

300 

Km 

SE 

8 

10 

6 

445 

Km 

SSE 

5 

11 

2 

65 

Km 

ESE 

3 

10 

7 

200 

Km 

SE 

10 

17 

23 

495 

Km 

SSE 

0 

0 

3* 

330 

Km 

SE 

2 

2 

5 

300 

Km 

SE 

1 

9 

9 

235 

Km 

SE 

27 

14 

10 

225 

Km 

NE 

99 

95 

62* 

320 

Km 

SE 

8 

6 

6 

520 

Km 

SE 

1 

0 

3* 

310 

Km 

E 

46 

26 

29 

550 

Km 

SE 

0 

1 

0 

175 

Km 

E 

5 

9 

3 

235 

Km 

NE 

38 

62 

32 

690 

561 

433 

(at  28 

loca- 

(at  30 

(at  33 

tions) 

! 

loca-  loca¬ 
tions)  tions) 

*  Magnitude  4.0 
**  Dotsero 

AEC  Plov/share  Program  Explosion  -  September  10,  1969,  Mag  5*0 
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can  possibly  be  accounted  for  by  the  fact  that  blastiniz  mav  ba  rnnHnr+oH 

shifts,  or  at  the  conclusion  of  the  working  day.  The 
dashed  line  through  23  events  across  the  plot  of  Figure  Improbably  re- 

level  of  naturally  occurring  earthquakes  throughout  the^day 

^  small  percentage  of  these  could  be  blasts,  but 

as  natural  working-day  hours  must  be  included 

as  natural  seismicity.  Frequently,  many  small  blasts  not  recorded  will 

etc.,  which  occur  naturally,  are  larger  than 
tnerefore  are  recorded  on  seismograms.  No  statistical  estimate 
of  the  percentage  of  blasts,  as  compared  with  natiorally  occurring  earth- 

pretation  accurate  because  of  the  imcertainty  of  inter- 

TABLE  2 


Events 

per  hour/day 

GCT 

1969 

1970 

1971 

00 

54 

28 

15 

01  ■ 

43 

14 

22 

02 

18 

15 

8 

.  03 

13 

5 

17 

04 

20 

9 

.  23  . 

05 

8 

15 

10 

06 

7 

li 

7 

5 

7 

4 

03 

4 

7 

8 

09 

18 

7 

7 

10 

6 

17 

7 

11 

9 

20 

■-  v>C 

'  •  9 

12 

12 

11 

8 

13 

.  5 

9 

5 

14 

i 

15 

10  . 

-  13 

15 

17 

8 

8  ‘  ' 

;  16 

12 

19 

11 

17 

26 

'  19 

4 

18 

44 

30 

■27 

19 

45 

32 

26  ^  \ 

20 

55 

40 

'  42  > 

*  ■*  * 

21 

.  72 

■  81 

*  57 

22 

102 

85 

59  ■ 

23 

80 

62 

36 

■  ?r 

/ 

*■  *  p.  i.. 

■ 
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Faulting  and  Crustal  Stress  at  Rangely,  Colorado 

C  B.  RaLeigh  and  J.  ri.  Healy 

U^,  Geological  Survey,  Menlo  Park,  California  94025 


J.  D.  Bredehoeft 

X/-S.  Geological  Survey,  Washington,  D.C.  20242 

In  Rangely,  Colorado,  the  maunitudes  and  directions  of  the  principal  stresses  have  been 
determined  from  hydraulic  fracturing  pressure  data  in  \\  tber  sandstone  at  the  depth  of  the 
eartliquake  foci.  The  total  principal  srros.=es  are  compressive;  Si  r=  590  bars  (SooO  psi), 
Ss  =:  430  bars  6200  i<i).  and  Sa  =  315  bars  (4550  psi).  The  direction  of  the  maximt 
compressive  stress,  N  TO^E,  is  'onsi.'tent  with  the  orientations  of  the  ma.ximuni  horizontal 
compressive  stress  measured  bj  ovorcoring  surface  exposures  of  sandstone  oi  the  Mesa 
Verde  formation.  The  eaivtiquakes  at  Rangely  occur  on  a  pre-existing  fault.  Tno  orientation 
of  the  fault  plane  and  tiie  slip  direction  are  known  from  focal  plane  solutions.  Resolving 
the  sUesscS  calo  the  'ujH  plane  in  the  direction  of  slip  gives  a  shear  stress  of  SC  bars 
(1120  psi)  and  a  total  normal  stress  of  350  bars  (5030  psi).  I.aboratorv'  data  on  the  initiation 
of  frictional  slidiug  on  saw  cuts  in  samples  of  V  eber  sand.stoue  and  the  r<''r.'t>utcd  she.ar 
and  normal  stresses  on  the  fault  at  Rang('ly  indicate  that  a  pore  pre.ssurti  of  260  bars  (.3130 
psi)  is  required  to  reduce  the  etioctivc  normal  stre.<5.s  sufucientiy  to  induce  slip.  Pore  pressures 
in  the  seismically  active  part  of  tlie  rc.scrvoir  at  a  time  of  fre<iucnt  earthquakes  were  275 
bars  (40C0  psi)  accord  with  the  calculated  pressure. 


David  Grig~s  first  suggested  during  the  de- 
uberat'orc  oi  the  ad  hoc  Panel  on  Earthquake 
Prediction  [Press,  1905]  that  earthquakes  might 
be  controlled  by  fluid  injection  into  a  fault 
zone.  This  radical  hypothesis  was  not  included 
in  the  panel's  final  report.  However,  without 
the  knowledge  of  the  panel,  the  U.S.  Army  had 
unintentionally  been  conducting  the  first  ex¬ 
periment  in  earthquake  control  with  ver>'  prom¬ 
ising  results. 

Earthquakes  were  triggered  by  injection  of 
waste  fluid  into  the  basement  rock  beneath  the 
Rocky  Mountain  Arsenal  of  the  Army  near 
Denver,  Colorado.  Evans  [19G0]  sliowed  con¬ 
vincingly  that  tlie  fluid  injection  wa.s  linked  to 
the  earthquakes  and  suggested  that  fluid  injec¬ 
tion  might  be  u.sed  to  control  earthquakes.  The 
ba-sis  for  the  hypothesis  of  Griggs  and  Evans 
regarding  earthq’'.ake  control  came  from  the 
work  of  Ilubbcrt  and  Rxihey  [1959],  who  showed 
that  the  fracture  strength  of  rock  is  propor¬ 
tional  to  the  difference  between  the  total  normal 
stress  across  the  fracture  and  the  pressure  of 


the  pore  fluids  within  the  rock.  Hcalu  et  al. 
[19G$]  .presented  evidence  that  the  Hubbert- 
Rubey  effect  provided  a  satisfactory  explan.a- 
tion  for  the  triggering  of  tiu?  Denver  earth¬ 
quakes  by  the  injection  of  fluid  at  high  p’'Cs  ar« 
into  rock  already  stressed  to  a  point  n^ar  its 
breaking  strength. 

Laboratory  experiments  on  the  fracture  of 
rocks  have  shown  that  failure  occurs  when  the 
shear  stress  on  some  potential  shear  surface 
exceeds  the  product  of  tlio  normal  stress  across 
it  (less  the  fluid  pressure)  and  a  friction.al  co¬ 
efficient  plus  the  cohesive  strength.  This  simple 
law  has  been  used  successfully  by  Dieterich 
[1971]  in  modeling  two-dimensional  fault.«  by 
finite  element  codes  in  digital  comjuircr.s.  The 
coefficient  of  friction  and  cohesive  strength  can 
be  determined  adt'quately  in  the  InboraroiA'.  but 
the  naturally  occurring  shear  and  normal  stresses 
can  only  be  determined  by  in  situ  measurements. 
Of  the  several  questions  left  unresolved  by  the 
incidence  of  the  Denver  earthquakes,  one  of  the 
most  important  and  experimentally  most  difll- 
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FAULTING  AND  STRESS  AT  RANGELY  > 


cult  was  that  of  the  state  of  stress  in  the  rocks. 
The  Denver  experiment  has  had  two  impor¬ 
tant  implications  aside  from  confirmation  of 
the  Hubbert-Rubey  theor\'.  Naturally  occurring 
earthquakes  might  be  modified  beneficially  by 
fluid  injection,  but  earthqn.ikcs  might  bo  trig¬ 
gered  inadvertently  by  fluid  injection  in  places 
where  no  hazardous  seismic  activity  would  nor¬ 
mally  be  anticipated.  If  earthquakes  arc  ever 
to  be  controlled  or  the  seismic  effects  associated 
with  subsurface  fluid  injection  predicted,  the 
stresses  within  the  rock  mass  must  be  known. 
This  paper  discusses  measurements  of  oticss  at 
the  Rangely  oil  field  and  their  bearing  on  the 
seismicity  associated  with  fluid  pressures  within 
the  reservoir  rock. 

EARTHQUAKES  AT  RANGELY 

Rangely,  Color.'>do,  is  the  site  of  a  detailed 
experiment  by  the  L'.S.  Geological  Survey  on 
the  relationship  between  fluid  pres.«’"e  and 
earthquakes.  Water  has  been  injecter’  by  the 
field’s  producers  into  the  resenmir  rock,  the 
Weber  sandstone,  for  secondary  reeovciw’  of  oil 
since  195S.  The  fluid  pressure?  in  the  pcripher>' 
of  the  field  have  risen  to  240-27,’)  bars,  about 
50%  in  excess  «.f  the  original  field  prcssiuo  of 
170  bars,  normal  hydrostatic  prt.'.-ure  for  the 
reserv'oir  depth  of  1.8  km.  Earthqu.akos  have 
been  occurring  in  the  vicinity  or  Rangely  at 
least  since  19b2.  when  the  Uinta  Basin  Scisrao- 
logical  Observatory  was  in.stnllcd  and  began  re¬ 
cording  [MiWfton,  1970],  Since  detailed  record¬ 
ing  of  the  seismic  activity  was  begun  by  the 
U.S.  Geological  Survey  in  19G9,  over  1000  earth¬ 
quakes  of  magnittule  Ml  greater  tlian  —0,5  have 
been  located  (Figure  1)  (C.  B.  Raleigh,  J.  H. 
Healy,  J.  D.  Brcdchoeft,  and  J.  Bohn,  unpub¬ 
lished  data,  1971), 

The  earthquakes  occur  along  the  southwest¬ 
ern  extension  of  a  fault  mapped  through  the 
field  by  Chevron  Oil  Company  geologists  (T. 
Larson,  personal  comintini cation.  1909)  from 
well-log  data.  Displacements  of  th.*  upper 
boundarj'  of  the  Weber  sandstone  along  the 
fault  arc  small,  about  15  metars.  As  is  shown 
in  Figure  1,  the  earthquakes  occur  only  along 
that  part  of  tlic  fault  at  which  pore  pressures 
are  greater  than  about  200  l)nrs.  Most  occur  at 
or  somewhat  below  the  depth  of  the  injection 
horizon,  the  Wchcr  sand.stone.  The  experiment 
is  designed  to  lower  and  raise  fluid  prc-ssures  in 


the  seismically  active  zone  to  determine  to  what 
extent  the  pore  pressure  controls  the  occur-  c, 
rence  of  the  earthquakes.  The  experiment  sh''uld 
provide  the  first  adequate  field  test  of  the  _ 
Hubbert-Rubey  hypothesis  if  the  fluid  pressure 
distribution  away  from  the  pumping-injcction 
wells  is  known  and  the  rock  stresses  are  meas-  , 
ured,  '*> 

When  the  orientation  and  the  magnitudes  of 
the  principal  stre.sscs  in  the  Weber  sandstone  ,  , 

are  given,  the  orientation  of  the  fault  plane  and  . 
the  slip  direction  still  must  be  determined  to 
c.alcul.ate  the  shear  and  normal  stresses  on  the 
active  fault  surface.  The  direction  of  the  prin-  ^ 
cipal  stresses  in  the  rock  should  be  consistent 
with  the  knov\Ti  sense  of  shear  (right  lateral)  > 

along  the  fault.  Fault  plane  solutions  have  been 
determined  from  the  sen.se  of  the  fir.st  motion  of 
P  waves  for  Rangely  event.s  recorded  on  10  or 
more  stations.  Tlic  fir>t  motions  are  distributed 
in  a  quadrar.tal  pattern  giving,  for  most  earth¬ 
quakes  (Figure  2),  strike  slip  on  vertical,  north¬ 
west-trending,  or  northe.ast-trcnding  nodal 
X)lanc3.  The  nodal  planes  showing  right  lateral 
strike  slip  are  nearly  parallel  to  the  mean  trend 
of  th#»  epicenters  and  thus  are  likdy  to  repre¬ 
sent  the  active  fault  planes.  This  trend  of  fault¬ 
ing  is  consistent  with  the  orientation  of  frac¬ 
turing  in  the  stirface  outcrops  of  ^I.ancos  sh.nle 
at  Rangely  [Thomas,  194.5]  and  the  surface 
faults  in  exposures  to  the  south  of  the  field. 

STRESS  MEASUREMENTS 

Two  modes  of  stress  meastirement  were  used 
at  Rangely.  Near-surface  stresses  were  meas¬ 
ured  by  strain  relief  techniques  in  the  Mesa 
Verde  sandstone  that  rims  the  Rangely  anti¬ 
cline  [de  la  Cruz  and  Raleigh,  1972].  Overcoring 
of  strain  gages  cemented  to  the  rock  surface 
and  of  the  U.S.  Bureau  of  Mines  borehole 
defonnation  gage  planted  in  sh.allow  holes  gave 
similar  results.  The  magnitudes  of  the  stresses,  j 
all  compressive,  were  small,  10  bars  or  so  [de  la 
Cruz  and  Raleigh,  1972],  bec.ause  tlicy  were 
apparently  relieved  by  surficial  weathering  and  ^ 
jointing.  The  directions  of  the  principal  hori- 
zunt.al  stress  components  wore  reasonably  con¬ 
sistent  from  measurerr  nts  at  three  different 
outcrops.  The  principal  directions  at  the  three 
sites  are  shown  in  Figure  3.  The  maximum  hori- 
zont.al  compressive  stress  is  within  20*  of  c.ist- 
west  from  averages  of  C-10  measurements  at 
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FAULTING  AND  STRESS  AT  RANG  ELY 


Fig.  2.  Well-determinod  focal  plane  solutions  for  Raneely  earth* (uaUes  for  October  1969 
to  October  1970.  Compre.'sional  first  motion.s  are  shown  in  blacE.  and  dilatational  nrst 
motions  in  white.  The  ro.<e  dinerum  givc.s  the  azimuths  of  the  right  lateral  strike  slip 
nodal  planes  from  each  solution.  The  d.ashed  line  is  the  field  boundary. 


each  jsite.  This  direction  is  consistent  with  both 
the  sense  of  shear  on  tiic  fault  determined  from 
the  focal  plane  solutions  and  the  greatest  com¬ 
pressive  stress  detennined  from  the  orientation 
of  an  induced  hydraulic  fracture  within  the 
Weber  sandstone  (N  TU'^E)  (Figure  3). 

The  other  technique  for  measurement  of  in 
situ,  stress  employed  at  Rangeiy  Is  h\uiauiic 
fracturing  of  a  borehole  into  the  Weber  sanu- 
stone.  The  method,  analyzed  by  Hubbert  and 
Willis  [1957],  Kchle  [1904],  and  Ilaivison  and 
Fairhurst  [1970],  permits  a  determination  of 


the  magnitude  of  the  raaNimum  and  least  hori¬ 
zontal  stresses,  provided  that  the  Took  is  net 
already  fractured  and  its  tensile  strength  is 
known.  As  perfonne*!  by  Ilahnson  [1972]  at 
Rai'igcly,  the  mcasuri-ments  involved  pai-ki: 
otf  a  section  of  the  open  hole  and  raising  ti.e 
fluid  p,»'f^sure  in  the  packed-ofif  «<»ction  un^il 
lUe  Tor::  fractured  in  tension.  When  an  unira*:- 
ture<i  cobh.m  of  core  5  nteters  lojig  was  removed 
during  drilling,  a  packer  was  set  5  meters  above 
the  bottom  of  the  hole  so  that  the  fr.'friunn 
could  be  accomplished  in  intact  ro*"k.  The  flui 


Fig.  3.  Map  showing  orientations  of  horizontal  princlj  il  strr.«503  dettrminod  from  ovor- 
-«nDg  (sing!c-slcmmc<l  arrows)  in  surface  expo.^'ures  of  t.;^  ^ai.A.one  ot  the  Me.'.a  Vf-rde 
formation,  from  earthquake  focal  plane  .•'olution.s  (double-stemmed  ojicn  arrows),  an<l  inun 
the  orientation  of  an  imluced  hydraulic  fracture  in  the  Weber  sandstone  (heavy  solid 
arrows). 
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well  rock  would  have  been  nrh^eved  over  the 
interval  of  less  than  hour  dtirin^  which  fluid 
pressure  was  raised  from  normal  hydrostatic  to 
the  breakdown  pressure.  With  regard  to  the 
second  condition  the  fracture  of  the  wall  of 
the  borehole  was  recorded  over  the  full  lencth 
on  both  sides  of  the  soft  inflatable  rubber  packer 
following  the  experiment.  Careful  examination 
of  the  core  from  that  section  of  the  hole  showed 
no  such  fracturing;  in  fact,  the  core  was  intact 
except  for  a  few  horizontal  partings  along 
shaley  Layers.  The  fracture  recorded  by  the 
packer  was,  therefore,  produced  by  pressuriza¬ 
tion  of  the  borehole. 

The  total  principal  stresses  calculated  frony 
(1)  are: 

St  “  590  bars  N  70®E  horizontal 
3%  =  427  bars  vertical  (assuming  the  litho- 
static  pressure  to  be  0.23  b/m 
depth) 

S,  »  314  bars  N  169®E  horizontal 

The  orientations  are  accurate  to  within  ±5®; 
the  fluid  pressure  measurements  during  fractur¬ 
ing,  from  which  the  magnitudes  of  the  stress 
are  derived,  are  believed  accurate  to  within 
dhl5  bars  on  the  basis  of  obs''n’ntions  from  two 
surface  pressure  transducers  and  an  Amerada 
gage  recording  at  the  bottom  of  the  hole. 

The  total  least  principal  stress,  being  equal 
to  the  ISIP,  is  very  well  documented  from 
numerous  other  tests  in  the  field.  During  routine 
hydraulic  fracturing  for  well  stimulation  at 
Rangely  both  the  ISIP  and  the  breakdown 
pressure  are  recorded.  The  initial  breakdown 
is  achieved  through  two  opposing  rows  of  ver¬ 
tical  perforations  in  the  casing,  and  the  near¬ 
in  fracture  orientation  may  be  more  controlled 
by  the  perforations  than  by  .the  far-field  stresses. 
\^atever  the  explan.ation  may  be,  it  is  observed 
that  the  breakdown  pressure  is  quite  variable  in 
such  tests.  However,  provided  that  the  fracture 
attains  an  orientation  normal  to  5»  at  short 
distances  from  the  hole,  the  ISIP  should  reflect 
the  discontinuous  decrease  in  tran.'jmissibility 
expected  when  the  crack  closes  as  the  fluid 
pressure  falls  to  the  magnitude  of  the  normal 
stress  Sj.  For  45  w’clls  in  the  central  and  eastern 
parts  of  the  field  the  ISIP  in  Ihese  more  re¬ 
cently  treated,  cased  holes  averaged  290  bars 
[Baimson,  1972].  Of  these,  35  wells  gave  an 
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ISIP  within  l\i%  of  that  v.alue.  The  ISIP  In  the 
western  end  of  the  firid  is  les.s,  averngine  220 
bars  for  12  mcastircments  [Ilnimson,  1972]. 
Therefore,  althoiigh  aiiditional  measmements  of 
the  breakdown  pressure  in  such  controlled  tests 
arc  desirable,  determinations  of  Sj  from  the 
ISIP  over  tlie  oil  field  arc  consi.-^icui,  and  in  good 
agreement  with  the  value  used  in  the  following 
analysis. 

EFFECTIVE  STRESS  DURING  FAULTING 

If  the  above  values  for  the  principn’  stresses 
and  the  orientation  of  the  fault  plane  and  slip 
direction  are  given,  the  shear  and  normal  stresses 
acting  iu  the  slip  direction  on  the  fault  plane 
in  the  vicinity  of  well  UP  72X31  (Figure  1) 
can  be  calculated.  The  nodal  plane  ami  slip 
direction  orientations  vary*  by  ±20®  for  differ¬ 
ent  earthquakes  along  different  parts  of  the 
fault  zone.  Furthermore,  the  measurements  of 
the  orientations  of  the  principal  horizontal 
stresses  \'ar>’  by  as  much  as  .±20*.  The  follow¬ 
ing  calculations  are  therefore  best  applied  Only 
to  the  near  vicinity  of  well  UP  72X31.  where 
the  detailed  stress  measurements  were  con¬ 
ducted.  The  orientations  of  the  nodal  planes  of 
the  nearest  eartliqtiakes,  O.S  km  to  the  south, 
are  plotted  in  Figure  i>  wl.ii  the  Ciiontations  of 
the  principal  stresses.  The  light  lat'^ral  nodal 
plane  parallel  to  the  trend  of  \he  epicenters 
strikes  N  50°E.  dipping  SO®  NW.  The  slip  direc¬ 
tion  plunges  20*  toward  N  2.34®E.  The  norm.*! 
to  the  fault  and  the  slip  direction  in  t’le  fault 
are  taken  as  the  coordinato  axes,  .‘i'  ami 
respectively,  and  the  normal  to  tlie  slip  direc¬ 
tion  in  the  fault  plane  is  taken  as  a:/.  Rotation 
of  the  principal  stresses  Sti  to  this  nc~*  ?•. - 
ordinate  system  gives  the  stresses  Sk/  in  tensor 
notation: 

5fci'  «  anatiSii  (2) 

The  direction  cosines  for  the  transformation 
are  determined  from  angles  between  the  coordi- 
.  natc  axes  measured  from  the  equal-area  projec¬ 
tion  (Figure  5).  . 
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direction  on  the  fault  plane  is  r. 

Tlie  normal  stress  S„  across  the  fault  plane  is 
=  S33'  Siflu*  +  Sta22^  +  ^3033* 

=*  347  bars 

The  shear  stress  t  along  the  slip  direction  is 


in  the  experimental  x\-ells  during  the  injection 
phase  of  the  experiment  when  earthquakes  were 
most  frequent  and  the  pore  pressure  was  2<5 
bars.  After  1  month  of  backllowing  the  experi¬ 
mental  wells  the  bottom  hole  pressure  dropped 
by  35  bars,  and  the  earthquakes  in  the  \'’cinity 
of  the  wells  ceased. 

It  is  useful  to  examine  the  effect  of  the 
variability  in  the  relative  orientations  of  the 
fault  and  the  stresses.  For  a  veitical  fault  plane 
that  is  slipping  in  the  horiccutoi  direction  and 
inclined  at  angle  a  to  S:,  the  pore  pressure  re¬ 
quired  for  faulting  can  readily  be  calculated. 
Figure  7  shows  a  plotted  against  the  pore  pres¬ 
sure  p  for  5i  =  590  bars,  =  314  bars,  and 
t/(S»  —  p)  =  O.Sl  at  failure.  The  variarion 
of  the  pore  pressure  re<iuired  to  initiate  faulting 
varies  between  about  235  and  200  bars,  where  a 
is  between  45’  and  55''  ;  The  second  cur  »  e  shows 
the  pore  pressure  variation  for  initial  stresses 
5i  —  0*0  biirs  and  Sj  —  ooO  bars.  Pro\  lued  mat 
the  magnitudes  of  the  stresses  moasureii  nre  cor¬ 
rect  to  trithin  about  15  bars,  the  pore  pressures 
existing  along  the  scismically  active  part  of  the 
fault  at  Rangely  are  quite  close  to  pressures 
calculated  by  Ilubbert  and  Rubey's  [1950]  ei- 
ffcvcive  st’-ejs  law  to  initiate  faulting  in  pre¬ 
viously  fracf’Ted  rock. 

CONCLUSIONS 


T  *=  Ti3*  =*  fluflaiSi  "b  012022*^2  4"  <*13^^23*5^3 

T«  77  bars  (4) 

Laboratory  experiments  on  the  strength  of 
Weber  sandstone  cores  have  been  conducted  by 
Byerlee  [1971].  The  data  in  Figure  6  show  the 
shear  stress  and  effective  normal  stresses  at 
failure  on  pre-existing  fracture  or  saw  cut  sur¬ 
faces.  The  reported  stresses  are  the  peak  stress 
encountered  where  sliding  begins  for  each  speci¬ 
men.  The  ratio  of  the  shear  stress  to  the  effec¬ 
tive  normal  stress  is  t/o-*  =  O.Sl  and  defines 
i  the  coefficient  of  friction  for  the  initiation  of 
slip  on  pre-existing  faults  in  Weber  sandstone. 

If  the  laboratory  data  are  applied  to  faulting 
•*  la  the  Weber  sandstone  at  Rangely,  slip  should 
take  place  when 

t/(S.  -  p)  =  0.81  (5) 

With  the  r  and  S,  values  from  hydraulic  fractur¬ 
ing  and  the  known  orientation  of  the  faulting, 
a  pore  pressure  p  of  257  bars  satisfies  (5).  This 
pressure  is  very  near  the  bottom  hole  pressure 


To  control  earthquake:  or  to  avoid  triggering 
earthquakes  inadvertently,  we  need  to  know  the 
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NORMAL  STRESS  IN  KB 

Fig.  6.  Shear  vcr*:3  effective  normal  stresjcs 
for  fracturing  of  intact  rock  and  slidins  on  .-aw 
cuts  in  Weber  sandstone  [Bycrlce,  1971].  Tri¬ 
angles,  fracture  of  intact  rock;  circles,  maximum 
friction  for  sliding  on  saw  cuts. 
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Fig.  7.  Pore  pressure  p  required  to  initiate 
shear  failure  in  Weber  sandstone  reservoir  as  a 
function  of  a,  the  angle  between  the  normal  to 
the 'fault  plane,  and  S,  the  maximum  prin<-ipal 
stress.  Si  and  Sj  arc  measured  as  590  and  314  bars, 
respectively  ;  values  of  670  and  330  bars  are  as¬ 
sumed  tc.snow  effect  of  a  variation  in  the  state 
of  stress. 

state  of  stress  in  the  rocks  and  a  criterion  for 
failure  that  takes  into  arronnt  the  chert  of 
pore  fluids  ou  rock  strength  and  behavior.  At 
Rangoly,  measurement  of  the  state  of  stress  in 
situ  has  becu  carried  out  by  hj'draulic  fractur¬ 
ing  at  the  depth  of  the  earthquake  foci.  In  the 
laboratory  the  -shear  strength  of  previously 
fractured  Weber  sandstone  has  been  determined 
by  Byerlee  [1971]  to  be  the  product  of  a  co- 
elHcient  of  friction  and  the  effective  normal 
stress  (the  total  normal  stress  less  the  pore 
pressure  of  the  ambient  fluid).  At  Rangcly,  if 
the  earthquakes  are  assumed  to  occur  when  slip 
is  initiated  on  a  pre-existing  fault,  the  resolu¬ 
tion  of  the  measured  stresses  onto  the  fault 
plane  in  the  direction  of  slip  indicates  that  a 
pore  pressure  of  at  least  255  bars  is  required 
for  shear  failure  in  tlie  fault.  Earthquakes  occur 
near  the  bottom  of  wells  at  which  the  pore 
.  pressure  is  v/ell  known  to  bo  around  275  bars. 
The  seismicity  ceased  after  the  pressure  in  the 
Wv!!?  was  dropped  by  35  bars.  The  obsor\*od 
.  pore  pressure  and  that  calculated  from  a  simple 
laboratory-derived  criterion  for  failure  corn 
blued  with  stresses  measured  in  situ  are  thus 
in  good  agreement. 

Tho  results  reported  here  suggest  that  the 
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measurements  of  stress  in  situ  at  considerable 
depth  can  be  made  and  tiiat  such  measurements 
can  be  used,  in  conjunction  with  the  h’boratorx'- 
determined  physical  properties,  to  predict  the 
pore  pressures  at  which  earthquakes  may  be 
triggered  or  propagating  fractures  possibly  ai- 
rested.  Although,  at  present,  little  is  known 
about  the  state  of  stress,  the  temperatures,  the 
poro  pressures,  or  tlie  strength  of  the  fault  zone 
material  in  active  faults,  it  appears  that  these 
quantities  are  measurable.  Moreover,  even 
though  adequate  failure  criteria  that  include 
the  effect  of  temperature  are  not  completely 
determined,  conventional  laboratory  methods  of 
obtaining  failure  criteria  by  experimentation 
can  be  applied  to  shear  failure  in  active  faults 
with  reasonable  hope  of  success.  As  a  result  of 
the  Rangely  experiment,  Griggs’s  speculation  on 
earthquake  control  can  now  be  tested  by  -an 
experiment  on  an  active  fault.  r  •  . 

APPENDIX 

Stress  directions  and  focal  plane  solutions. 
One  interesting  result  of  this  study  is  the  ob¬ 
servation  that  the  directions  of  the  prlnclivil 
horizontal  stresses  from  in  situ  measurements 
agree  with  those  derived  from  tl’.c*  '^oc"’  plane 
solutions.  The  radiation  pattern  nf  P  waves 
from  each  of  these  earthquakes  defines  a  double 
couple  solution  for  which  the  nodal  planes 
represent,  in  effect,  mutually  perpendicular  con¬ 
jugate  shear  planes.  The  maximum  end  least 
principal  stresses  gudng  the  maximum  shearing 
stress  on  these  planes  in  the  proper  direction 
bisect  the  dilatational  and  compressionai  quad¬ 
rants,  respectively,  of  the  focal  sphere  [Schei- 
degger,  1964].  For  the  unambiguous  solutions 
for  Rangely  earthquakes  the  maximum  and  least 
compressive  stresses  arc  derived  accordingly  and 
shown  in  Figure  S.  The  maximum  and  least 
horizontal  stresses  from  the  in  situ  measure¬ 
ments  plotted  on  the  same  projection  agree 
.  quite  well  with  those  from  the  earthquakes. 

McKenzie  [1969a]  has  pointed  out  that  pre¬ 
existing  faults  may  slip  even  at  quite  low  re¬ 
solved  shear  stresses  before  new  fractures  in 
intact  rock  will  develop  [also,  cf.  llandin, 
1969].  His  analysis  indic.atcs  that  the  P  direc¬ 
tion,  otSu  the  maximum  principal  stress  in  the 
present  terminology,  may  lie  anywhere  within 
the  dilatational  quadrant  of  tlic  focal  sphere.  If 
this  analysis  b  correct,  stress  directions  derived 
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Fig.  8.  P  (laaximum  principal  .'tre-'S,  compres¬ 
sion)  directions  (circles)  and  T  (least  principal 
stress,  tension)  directions  (plus  signs)  from  fault 
piano  solutions  of  Rangely  eartlu'iuakes.  The 
arrows  on  the  circle  of  projection  are  the  greatest 
and  least  horizontal  principal  stresses  from  in 
situ  stress  measurements. 


from  focal  plane  solutions  would  clearly  be  of 
only  limited  tectonic  significance. 

However,  even  intact  homogeneous  rock  has 
finite  strength,  so  that  new  faults  can  be  gen¬ 
erated  even  though  existing  faults  are  under 
finite  shear  stress. 

The  shear  strength  r,  of  a  fracture  in  rock  is 
given  by  the  product  of  the  effective  normal 
stress  times  a  coefficient  of  friction,  i.e.,  for 
normal  stresses  corresponding  to  shallow  earth¬ 
quakes, 


r,  «  <r,M  (Al) 

Likewise,  the  shear  strength  r/  of  intact  rock 
is  given  by  the  normal  stress  times  the  coeffi¬ 
cient  of  internal  friction  tan  plus  a  cohesive 
strength  r,: 

r/  =*  To  +  <r»  tan  <i>  (A2) 

These  equation,  are  represented  by  straight 
lines  on  the  Mohr  diagnmi,  a  plot  of  t  versus 
<r„  the  shear  and  nonnal  stresses.  With,  for 
example,  the  data  for  Weber  B  sandstone  from 
Byerlee  [1971]  in  kilobars. 


T.  =  0.81<r,  (A3) 

t/  «=  0.7  +  (A4) 

Failure  by  faulting  in  the  intact  Weber  sand¬ 
stone  will  occur  when  i  ^^ohr  circle  (Figiire  9) 
is  tangent  to  the  envelope  for  r/,  or/.  The 
circle  has  its  center  at  (o-i  -f-  a»)/2  and  is  the 
locus  of  points  where  r  and  cr,  are  given  by 


r  =  •*  " — -  sin  2a  (A5) 

O’l'f'O’s  ,  O’]  CTs _  /  K  C\ 

<r.  == - - 1 - - cos  2a  [A 6] 

Thus,  for  stresses  o-i  and  <r,,  jequired  for  fractur¬ 
ing  intact  material,  tlicre  will  be  a  range  of  <x 
over  which  slip  on  a  pre-existing  fracture  will 
take  place  (Figure  OF  The  values  of  2rt  at  the 
intercept?  of  the  Molir  circle  on  tlie  t.’.veicpe 
for  sliding  on  existing  faults  are  given  by 

sin  2a  ^  1  2rrt 

- cos  2a  =  "  —  7 - r~  T 

H  sm  <t>  (cTj  —  (T3)  tan  q> 

(A7) 

Two  values  for  a,  at,  and  a,  satisfy  the  above 
equation  hud 

p  =  —  cot  (ai  aj)  (AS) 

The  angles  at  and  a,  can  also  be  solved  graphi¬ 
cally  ir  the  Mohr  con.«'ruction  (Figure  9). 
Table  1  gives  the  ranine  of  a  over  which  faulting 
will  take  place  on  existing  faults  in  ’^reference 
to  fracture  in  intact  Weber  sandstone  as  a 
function  of  the  stresses  (for  simplicity,  <r,  >  a. 


Fig.  9.  Mohr  diacram  showing  envrlnprs  for 
failure  in  intact  Weber  sandstone  .and  for  slip  | 
on  pre-existing  faults  in  Weber  sandstone  an<l  a 
Mohr  circle  tangent  to  the  envelope  for  failure 
in  the  intact  rock. 
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TABLE  1.  R-inse  of  AorIcs,  a,  ax  <  a  <  between  ^und  the  Normal  to  the  Fault  Phuie  at  Strer=ses  over 
Which  Slip  Will  Occur  on  a  Pre-existine  Fatilt  Ibither  Than  on  a  New  Shear  Fracture  in  Intact  Rock* 


St*’es3,  kb 

Range,  deg 

(<ri  +  <rj)/2 

«ri 

<r» 

on 

at  —  ax 

1.6 

2.95 

0.25 

2.0 

3.60 

0.40 

2.5 

4.38 

0.62 

3.0 

5.12 

0.88 

3.5 

5.88 

1.12 

4.0 

6.61 

1.39 

4.5 

7.38 

1.62 

44 

8.5 

41 

46 

ai 

38 

48 

81 

3.3 

51 

78 

27 

54 

76 

•>2 

56 

.  72 

16 

59 

•70 

11 

% 


*  Calculated  for  Weber  B  sandstone. 


=  a,).  The  range  is  greatest  at  low  stress. 
From  Table  1,  o-i  (or  tlie  P  axis  in  McKenzie's 
[19G9a]  terminology)  may  lie  anywhere  be¬ 
tween  44°  and  S5°  to  tiie  fault  plane  normal 
for  slip  to  occur  preferentially  alon"  the  fault 
at  stresses  ci  =  2.95  kb  and  u,  =  0.25  kb.  The 
mean  stress  is  1.15  kb,  corresponding  approxi¬ 
mately  to  a  depth  of  5  km.  At  greater  depths 
the  higher  stresses  reejuired  for  sliding  on  exist¬ 
ing  fractures  approach  those  requiretl  for  irac- 
turmg  in  intact  rock  because  ,u  >  tan  <}>. 

A  more  extreme  example  is  the  Westerly 
granite,  which  has  a  very  high  fiacturt  c. renuih 
but  a  frictional  sliding  coefficient  verx*  do.ic  to 
that  of  the  Weber  sandstone  {Bijerlce.  1007]. 
When  ffi  lies  less  than  40°  or  greater  than  S2° 
to  a  pre-existing  fracture,  faulting  wilt  occur  in 
the  intact  rock  in  preference  to  frictional  sli<ling. 
These  angles  hold  over  a  wide  range  of  stress. 


The  rna.ximum  principal  .'tress  is  taken  to  bj 
45°  from  the  fault  plane  wh<^ii  derived  from 
focal  pl.ane  solutions.  This  arigle  could  be  in 
error  by  as  much  as  40°  where  sliding  on  an 
existiiu:  fault  in  an  otherwise  in'act  rock  pro¬ 
duces  the  earthquake.  If  the  nou:ii  plane  corre-  - 
spending  to  the'  fault  is  known,  sr.ould  b? 
plotted  at  00°  to  the  normal  to  the  fault  plaiK 
and  30°  to  the  slip  direction.  In  this  case  the 
o'ientation  of  would  not  be  in  error  by  more 
than  20°. 
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leauue^.  J.  D.  Byf'rlee,  K.  cle  la  Cruz,  and  J. 
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criticisms.  Publication  of  this  paper  was  author¬ 
ized  by  the  Director  of  the  U-'i.  GooluW'’ul  ia.M  y. 
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Departm''nt  of  Defense  under  ARP-\  orders  1169 
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PREFACE 


The  objective  of  Project  Rio  Blcirico  'wsis  to  8timvil3.te^  of  rmcle^r 

explosives,  the  production  of  natural  gas  from  a  1, 300-foot  section  of  the 
Fort  Union  and  Mesaverde  formations  in  the  Piceance  Basin  of  Northwest 
Colorado. 

The  three  nuclear  explosives  for  Project  Rio  Blanco  were  detonated  at 
1000:00.  12  hr  -0.  01  second.  Mountain  Daylight  Time,  or  1600:00.  12  hr 
to.  01  second,  Greenwich  Mean  Time,  on  May  17,  1973.  The  three 
explosions  occurred  within  the  Fort  Union  and  Mesaverde  formations  at 
depths  of  5,  838.  5  feet,  6,  229.  7  feet,  and  6,  689.  5  feet.  The  three  explo¬ 
sions  occurred  nearly  simultaneously  as  planned  and  were  completely 
contained.  The  preliminary  indications  are  that  the  yields  of  the  three 
explosives  totalled  approximately  90  kt.  The  elevation  of  the  ground  at 
the  emplacement  well,  RB-E-01,  is  6,629.  9  feet  above  mean  sea  level. 
RB-E-01  is  located  1,  080.  50  feet  south  of  the  north  line  and  1,  188.  49 
feet  east  of  the  west  line  in  Section  14,  Township  3  South,  Range  98  West 
of  6th  P.  M. ,  Rio  Blanco  County,  Colorado,  which  corresponds  to  geodetic 
coordinates  of  108°2r59"  west  longitude  and  39°47'35"  north  latitude. 


RBOSP  PROf3R£SS  REPORT 


AS-48 


ML 


- 


rV. 


I 


5* 


f 


f 


I. 


1 

t. 


V.  AS-49 

s 

4  ^ 


i:: 

\ 


% 


V 


TABLE  OF  CONTENTS 


Page 


PREFACE  .  ^ 

1.  introduction  .  1 

2.  ground  motion .  ^ 

2.  1  Instrumentation  .  . .  •  •  • 

2.2  Data  . .  **•  ^ 

3.  seismicity  and  ground  motion . . .  23 

3.  1  Seismicity  . 

3,  2  Correlation  of  Structural  Damage  to  Low  Level 

Ground  Motion  ...  ...  . . 

3,  3  Structural  Response  to  Ground  Motion  ...  . .  •  •  •  •  29 

3,4  Ground  Spall  Documentation  Program  ...  .••  ••• 

3  5  Sandia  Laboratories  Pore  Pressure  and  Ground 

*  30 

Motion  Measurements . 

30 

3.  6  Structural  Bracing . 

3.7  Bridges  . 

3.  8  Oil  and  Gas  Wells . 

3.9  Rock  Falla  . 

3.10  Wells  and  Springs . 

3.  11  ■  Microtremor  Surveys  . 

3.  12  El  Paso  Natural  Gas  Company  Program  . 

3.13  Mines  and  Quarries  ...  . . . 

3.  14  Damage  Claim  and  Adjustments  ^3 

3.15  Structural  Inventory  . . ' 

4.  CONCLUSIONS  . . 

45 

4.1  General  . . . 

4.2  Ground  Motion  .••  . . .  •** 

4.  3  Aftershocks . ^ 

4Q 

REFERENCES . 


RBOS 


P  PROGRESS  REPORT 


J 


AS-50 


I 


ILLUSTRATIONS 


Figure 


Page 


L 

2. 

3. 

4. 

5. 

6. 


7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 
19. 


Ground  Motion  Recording  Sites . 

PHC  of  Velocity  vs  Slant  Distance,  Hard- Rock 

Stations  . 

Predicted  and  Observed  PHC -of- Velocity  Values, 

Hard- Rock  Stations  . 

PHC  of  Velocity  vs  Slant  Distance,  Unknown  or 

Alluvium  Stations . 

PHC  of  Acceleration  vs  Slant  Distance,  Hard- Rock 

Stations  . . . .  ... 

Predicted  and  Observed  PHC-of- Acceleration 
Values  for  Hard-Rock  Stations  ...  ...  .  . . 

PHC  of  Acceleration  vs  Slant  Distance,  Unknown 

or  Alluvium  Stations  . 

Rio  Blanco  Aftershock  Hypocenter  Location  Stations 
Number  of  Aftershocks  per  5-Minute  Interval  . . . 
Number  of  Aftershocks  per  1-Hour  Interval  after 

Detonation .  . 

Location  of  Epicenters  . ... 

Depth  Distribution  of  Hypocenters  ...  . 

Rockfall  Locations . ... 

Rifle  Microtremor  Survey  . 

Grand  Valley  Microtremor  Survey  . 

De  Beque  Microtremor  Survey . 

Rangely  Microtremor  Survey . 

Meeker  Microtremor  Survey  . . ... 

Meeker  Microtremor  Survey  Results  vs  Damage 
Distribution . 


•  •  • 


•  •  • 

•  •  • 

•  •  • 

•  •  • 


•  «  • 


4 

12 

13 

14 

15 

16 

17 

24 

25 

26 

27 

28 
34 

36 

37 

38 

39 

40 

41 


I  R30SP  PH0GRE:SS  RliPoSTj 


AS-51 


( 

K 


i 


TABLES 


1. 

2. 

3. 

4. 

5. 


Ground  Motion  Data  . .  •••  ' . 

Regression  Data  .  ...  ...  ••.  ••• 

Comparison  of  Rulison  and  Rio  Blanco  Data  from 
the  Same  Recording  Locations  ...  ...  ... 

Property  Damage  Claims  ...  ...  ...  ••• 

Comparison  of  Effects  Predictions  and  Actual  Effects 


Page 

5 

11 

19 

44 

46 


[rBQSP  PRQGFIESS  REPORt'J 


AS-52 


! 


i". 

I  ■ 

r 


J 


1.  INTRODUCTION 


Several  seismic  effects  programs  were  carried  out  in  connection  with  Pro¬ 
ject  Rio  Blanco.  They  were: 

1.  A  seismic  array  to  monitor  the  ground  motion  resulting  from  the 
detonation  at  distances  varying  from  3,7  miles  to  93  miles  from 
the  emplacement  well  (EW). 

2.  Two  seismic  arrays  to  detect  the  aftershocks  of  the  detonation  and 
to  locate  the  hypocenters  of  these  aftershocks, 

3.  A  seismic  array  to  study  the  correlation  of  structural  damage  to 
low-level  ground  motion, 

4.  A  seismic  array  to  study  structural  response  to  the  ground  motion, 

» 

5.  A  seismic  array  to  monitor  the  ground  spall  caused  by  the  detona¬ 
tion.  1^) 

6.  An  instrumental  array  to  monitor  changes  in  pore  pressure  result¬ 
ing  from  the  seismic  wave.  (H) 

7.  A  structural  bracing  program  to  minimize  the  damage  caused  by 
the  ground  motion. 

8.  A  pre-  and  postdetonation  inspection  of  selected  bridges  in  the 
area  to  detect  any  change  due  to  ground  motion, (12) 

9.  A  pre-  and  poatdetonation  inspection  of  the  oil  and  gas  wells  within 
7  miles  of  the  EW,  (1) 

10,  A  postdetonation  inspection  to  locate  and  document  the  reported 
rock  falls  due  to  the  ground  motion,  (1) 

11,  A  pre-  and  postdetonation  inspection  of  springs  and  wells  out  to  a 
distance  of  10  miles  from  the  EW,(1»  1^) 

12,  Predetonation  microtremor  surveys  in  Rifle,  Meeker,  Rangely, 

De  Beque,  and  Grand  Valley  v/ith  postdetonation  correlation  to 
damage  claims,  (1) 

13,  El  Paso  Natural  Gas  Program.  The  primary  objective  of  this  pro¬ 
gram  was  to  improve  the  accuracy  of  effects  predictions  of  under¬ 
ground  nuclear  explosion  induced  ground  motion  on  manmade  and 
natural  structures  and  the  bioenvironment  near  the  detonation  point. 
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Project  Rio  Blanco  provided  an  opportunity  to  improve  capabilities 
for  evaluations  relative  to  Project  Wagon  Wheel.  Tlie  specific 
objectives  were  to  monitor  pore  pressures  for  potential  sand  liqui- 
faction,  monitor  the  effect  on  small  bridges,  determine  the  poten¬ 
tial  for  rock  falls,  measure  wavelengths  in  the  ground  and  monitor 
the  effect  on  vegetation,  and  monitor  the  effect  upon  mechanical 
structures  characteristic  of  conventional  gas  field  operations. (2) 

14.  A  pre-  and  postdetonation  inspection  was  made  of  all  the  mines 
and  quarries  within  53  miles  of  the  EW.  As  a  safety  precaution, 
all  mines  and  quarries  in  this  area  were  evacuated  at  detonation 
time.  ) 

«  i 

15.  Damage  claims  and  adjustments.  (1 ) 

16.  Structural  inventory  and  damage  assessment.  An  inventory  was 
made  of  all  structures  within  28  miles  of  the  EW  and  assessments 
made  of  the  seismic  hazard,  A  prediction  was  made  of  the  ex¬ 
pected  damage  In  dollars. 


RBO'SP 


AS-55 


r  • 


I 


2.  GROUND  MOTION 


2.  1  INSTRUMENTATION 

Thirty- seven  locations  were  instrumented  for  ground-motion  documenta¬ 
tion  with  Li-7  velocity  systems  by  the  U,  S,  Geological  Survey /Seismic 
Engineering  Branch  (USGS/SEB).  Eighteen  of  these  locations  were  for 
CER  Geonuclear  Corporation's  (CER)  basic  documentation  program;  five 
were  for  the  structural  response  studies  conducted  by  Kenneth  Medearis 
&  Associates  (KMA);  ten  were  installed  in  five  close-in  communities  for 
a  ground-motion  structural-damage  correlation  study  conducted  by  John 
A.  Blume  &  Associates,  Engineers  (JABAE);  one  location  at  the  Colony 
Tower  was  instrumented  for  the  Atomic  Energy  Commission  (AEC);  and 
one  location  at  the  base  of  Douglas  Pass  was  instrumented  for  the  USGS/ 
SEB  program. 

Two  locations  were  instrumented  for  ground-motion  documentation  with 
strong-motion  accelerographs;  the  Shell  Plant  location  with  a  Teledyne- 
Geotech  RFT  250  and  the  Cascade  Pump  Station  with  a  Kinemetrics  SMA-1. 
The  Hayden  Power  station  was  instrumented  by  KMA  with  an  accelerograph. 
Several  strong  motion  stations  were  operated  by  Atlantic  Richfield  Co,  , 
two  of  which  were  reported,  Arco  No.  1  and  Arco  No,  10,  Arco  No,  1 
station  used  an  Endevco  accelerometer  while  Arco  No,  10  used  a  Geotech 
RFT  250  strong-motion  accelerometer.  Figure  1  shows  the  location  of 
all  these  recording  sites. 


2.  2  DATA 

Where  the  data  were  recorded  as  particle  velocity  (L-7  stations),  the 
records  were  digitized  and  differentiated  to  obtain  accelerations;  for 
those  stations  where  the  data  were  recorded  as  accelerations,  the  records 
were  digitized  and  integrated  to  obtain  velocity.  These  data  are  shown  in 
Table  1. 

Regression  analyses  were  made  for  the  peak  horizontal  component  (PHC) 
of  velocity  and  the  PHC  of  acceleration.  The  data  we^e  separated  as  to 
surface  conditions  into  hard  rock,  alluvium,  or  unknowns  using  the  follow¬ 
ing  criteria:  if  the  site  was  visibly  on  hard  rock  (consolidated  material 
such  as  sandstone  or  shale)  or  microtremor  measurement  showed  a  reso¬ 
nant  frequency  of  15  Hz  or  greater,  the  site  v/as  classified  as  hard  rock. 

If  there  was  geologic  data  indicating  the  presence  of  more  than  5  ft  of 
alluvium  or  if  a  microtremor  measurement  at  the  site  showed  a  resonant 
frequency  of  less  than  15  Hz,  the  site  was  classified  as  alluvium.  When 
geologic  data,  microtremor  data,  or  visual  data  were  not  available,  the 
surface  condition  at  the  site  was  classified  as  unknown. 
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Figure  1. 


Ground  motion  recording  sites. 
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Table  1.  Ground  motion  data. 


t 


I 

I 

C*  j 


Slant 

o 

Peak 

Peak 

Station 

Distance 

m 

Velocity 

Accel 

Station 

Medium 

(km) 

P 

(cm/ sec) 

(s) 

Comments 

Arco  ^1 

Unknown 

4.  34 

V 

2.6-9 

Spall 

T 

— 

Strong  motion  data 

R 

— 

— 

Endevco  accelerometer 

SR-5 

Unknown 

5.  26 

V 

T 

R 

V 

46.3 

— 

Spall 

RB-1 

Hard  Rock 

6.  00 

? 

Spall 

Equity  Plant 

T 

25.  8? 

— 

Tape  recorder 

R 

47.  5? 

— 

affected  by  ground  motion 

RB-2 

Hard  Rock 

7.74 

V 

18.4? 

1.44? 

Tape  recorder 

White  River 

T 

13.  3? 

0.  967? 

affected  by  ground  motion 

Antenna 

L 

14.  2? 

1.73? 

RB-3 

Hard  Rock 

11.  3 

V 

5.  05 

0.  337 

Rock  School 

■ 

T 

4.  54 

0.  296 

L 

4.  34 

0.  282 

Arco  #10 

Unknown 

13.  0 

V 

— 

0.  35 

Strong  motion  data 

T 

— 

0.39 

Geotech 

R 

— 

0.  40 

• 

RB-4 

Alluvium 

17.3 

V 

8.  75 

0.  204 

Strong  motion  data 

Shell  Plant 

T 

8.4 

0.  136 

Geotech  R.  F.  T.  250 

R 

9.  5 

0.156 

RB-5 

Hard  Rock 

18.9 

V 

7.85 

0.  33 

Chadbourne 

T 

2.  58 

0.  132 

Plant 

R 

4.  05 

0.  162 

RB-8 

Unknown 

22.  8 

V 

4.  22 

• 

0.244 

Weiland 

T 

5.  28 

0.186 

Ranch 

R 

8.85 

0.434 

RB-6 

Unknown 

25.  3 

V 

5.4 

0.  144 

Strong  motion  data 

Cascade 

T 

8.9 

0.  193 

Kinnemetrics  SMA-1 

Pump  Station 

R 

6.8 

0.227 

V  is  vertical,  T  is  transverse,  R  is  radial,  L  is  longitudinal  to  structure. 
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Table  1.  Ground  motion  data.  (Cont.  ) 


C* 


Station 

Station 

Medium 

Slant 

Distance 

(km) 

o 

m 

P 

Peak 

Velocity 

(cm/sec) 

Peak 

Accel 

is). 

Comments 

RB-7 

Alluvium 

29.6 

V 

5.08 

0.  203 

i 

Colony  Tower 

T 

5.  34 

0.  154 

L 

4.  08 

0.  124 

RB-7a 

Hard  Rock 

29.8 

V 

1.  27 

0.  043 

Inside  mine 

Colony  Mine 

T 

0.90 

0.  029 

Not  a  surface  location 

L 

0.75 

0.  024 

RB‘9 

Unknown 

35.9 

V 

1.78 

0.  069 

Rfo  Blanco 

T 

1.  35 

0.  035 

Dam 

R 

1.  23 

0.  038 

RD-10 

Hard  Rock 

43.  4 

V 

0.  585 

0. 0199 

Douglas  Pass 

T 

0.  576 

0. 0254 

FAA  Warehouse 

R 

0.  532 

0.  0246 

RB-ll 

Hard  Rock 

45.7 

V 

0.294 

0. 0075 

Inside  mine 

Mobil  Mine 

T 

0.  235 

0. 0067 

Not  a  surface  location 

R 

0.  342 

0.  0125 

RB-U 

Hard  Rock 

46.4 

V 

0.  447 

0. 0156 

Douglas  Pass 

T 

0.414 

0.  0106 

Base 

R 

0.  518 

0. 0094 

NSG-1 

Unknown 

46.  1 

V 

0.  509 

0.  020 

Grand  Valley 

T 

0.  474 

0.  016 

North 

R 

0.  933 

0.  024 

RB-12 

Hard  Rock 

46.4 

V 

0.296 

0.  0105 

Grand  Valley 

T 

0.466 

0.  014 

High  School 

R 

0.  518 

0.  014 

NSG’-  2 

Unknown 

46.  5 

V 

0.  309 

0. 0097 

Grand  Valley 

T 

0,  290 

0. 0073 

South 

R 

0.436 

0. 0130 

RB-13 

Hard  Rock 

46.9 

V 

0,  466 

0.019 

I  BBOSP  PROGPiESo  REPORT  "j 

Meeker 

T 

0,  206 

0. 0097 

Coop,  Store 

R 

0.674 

0.  024 

V  is  vertical,  T  is  transverse,  R  is  radial,  Lis  longitudinal  to  structure. 
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Table  1,  Ground  motion  data,  (Cont.  ) 

C* 

Slant  o  Peak  Peak 


Station 

Station 

Medium 

Distance 

(km) 

m 

Velocity 
(cm/ sec) 

Accel 

(g)  Comments 

NSM-1 

Unknown 

47.  6 

V 

0.  558 

0,  0205 

Meeker 

T 

0.  510 

0. 0216 

Hospital 

R 

0.  700 

0. 0298 

NSM-2 

Unknown 

48.  1 

V 

0.  576 

0.  0267 

Meeker  High 

T 

0.  584 

0.0242 

School 

R 

0.  888 

0.  0326 

RB-26E 

Unknown 

48.  1 

V 

—  . 

Rangely  Jr, 

T 

0.  403 

0. 0146 

College 

L 

0.  362 

0.  0122 

RB-26W 

Unknown 

48.  1 

V 

«  i-  «• 

.  . 

Rangely  Jr, 

T 

0.  305 

0. 0146 

College 

L 

0.  290 

0. 0100 

RB-14 

Hard  Rock 

48.9 

V 

0.  306 

0. 0114 

J,  L,  Supply 

T 

0.  301 

0. 0142 

Rangely 

R 

0.  371 

0. 0191 

NSA-2 

Unknown 

49.5 

V 

0.  520 

0. 0231 

Rangely 

T 

0.  665 

0. 0253 

South 

R 

0.  984 

0. 0303 

NSA-1 

Unknown 

50.  2 

V 

0.  348 

0.0123 

Rangely 

T 

0.  282 

0.0127 

North 

R 

0.480 

0.0203 

RB-15 

Hard  Rock 

52.4 

V 

0.331 

0.0103 

DeBeque 

T 

0.  239 

0. 0093 

School 

R 

0.  203 

0. 0084 

NSD-2 

U  nknown 

52.  5 

V 

0.  316 

0. 0090 

Debeque 

T 

0.  243 

0. 0117. 

West 

R 

0.  250 

0. 0088 

NSD-1 

U  nknown 

52. 7  . 

V 

0.  335 

0.0098 

DeBeque 

T 

0.  256 

0.0126 

East 

R 

0.280 

0.0104  j  RBOSP  PROGRtSS  REPORT 

NSR-2 

Unknown 

57.2 

V 

0.441 

0. 0208 

Rifle 

T 

0.  706 

0. 0160 

West 

• 

R 

0.  744 

0.  0192 

=5^  V  is  vertical. 

T  is  transverse,  R  is 

radial,  L  is 

longitudinal  to  structure. 
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Table  1.  Ground  motion  data,  (Cont.  ) 

C* 


'  ' 

Slant 

o 

Peak 

Peak 

Station 

Distance 

m 

Velocity 

Accel 

station 

Medium 

(km) 

P 

(cm/ sec) 

(g) 

RB-17 

Unknown 

57.  3 

V 

0.  315 

0. 0161 

Millers  Car 

T 

0.  254 

0.  0178 

Wash,  Rifle 

• 

R 

0.391 

0. 0206 

RB-16 

Hard  Rock 

58.  1 

V 

0.  276 

0. 0078 

Methodist 

T 

0.  253 

0. 0086 

Church,  Rifle 

1 

R 

0.  382 

0. 0137 

NSR-1 

Unknown 

58.  4 

V 

0.  695 

0.  0344 

Rifle 

T 

0.  673 

0.  0190 

East 

R 

0.  745 

0. 0181 

RB-18 

Alluvium 

64.  1 

V 

0.  339 

0. 0143 

Harvey  Gap 

T 

0.452 

0.  0168 

Dam 

L 

0.624 

0.  0185 

RB-19 

Unknown 

71.  5 

V 

0.  181 

0. 0127 

Cameo 

T 

0.369 

0. 0156 

R 

0.  382 

0.  0183 

RB-20 

Hard  Rock 

77.  1 

V 

0.  1325 

0. 00513 

Little 

T 

0. 0705 

0. 00254 

Bonanza,  Utah 

R 

0. 0955 

0. 00338 

RB-21 

Unknown 

82.  5 

V 

0.  125 

0.  0069 

Grand  Junction 

T 

0.  316 

0. 0106 

CER  Office 

R 

0.  264 

0. 0079 

RB-23 

Hard  Rock 

93.  5 

V 

0.  0296 

0. 00203 

Glenwood 

• 

T 

0.  0601 

0. 00368 

Springs 

R 

0. 0592 

0. 00325 

RB-22 

Unknown 

106.  1 

V 

0.  068 

0. 00175 

Craig 

■ 

T 

0.  132 

0.  00483 

R 

0.  124 

0. 00496 

Hayden 

Unknown 

129. 00 

V 

-  — 

0.  0012 

Comments 


T 

R 
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«  V  is  vertical,  T  is  transverse,  R  is  radial,  L  is  longitudinal  to  structure. 


AS-61 


Table  1.  Ground  motion  data.  (Cont.) 


C* 


r  • 

Station 

Station 

Medium 

Slant  ' 
Distance 
(km) 

o 

m 

P 

Peak 
Velocity 
(cm/ sec) 

Peak 

Accel 

i£] 

Comments 

RB-24 

Aspen 

Unknown 

148.  9 

V 

T 

R 

0.  117 

0.  115 

0.  156 

0. 00343 

0. 00284 

0. 00321 

- 

Note:  Similar  sets  of  data  have  been  published  by  USGS  in  PNE-RB-53  and  the 
Environmental  Research  Corporation  in  PNE- RB- 54. All  of  these  sets  of  data 
are  based  on  the  same  set  of  recordings.  Variations  of  as  much  as  10%  can  be 
expected  due  to  differences  in  the  individual  reading  of  the  recordings. 


V  is  vertical,  T  is  transverse,  R  is  radial,  L  is  longitudinal  to  structure. 
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The  CER  predetonation  predictions^ ^  ^  were  based  on  the  PHC  of  velocity 
and  acceleration  from  hard-rock  station  locations.  The  predictions  were 
based  on  a  two-line  regression  curve.  One  line  was  used  for  distances  of 
from  6  km  to  22  km  and  the  second  line  was  used  for  distances  of  from 
22  km  to  300  km.  The  rationale  for  using  the  peak  horizontal  component 
of  the  motion  is  that  structures,  in  general,  have  far  greater  resistance 
to  damage  from  vertical  motion  than  from  horizontal  motion. 

In  the  earlv  reporting  of  the  Rio  Blanco  results,  a  two-line  regression 
curve  was  assumed^^^.  Later  analysis,  however,  indicates  that  all 
Uald  Irom  sictiiuiii  loser  than  lU  km  are  highly  questionable,  bet  ause  of 
spall  or  the  ground  motion  affecting  the  tape  recorder.  In  this  report,  a 
single-line  regression  curve  is  assumed  from  11  km  to  150  km  for  both 
hard- rock  and  unknown  station  condition  analyses. 

Table  2  sh  ows  the  stations  and  the  data  used  in  the  regression  analyses. 
Figure  2  shows  the  results  of  the  regression  analysis  of  the  PHC  of 
velocity  vs  slant  distance  for  hard- rock  stations. 

Figure  3  shows  the  comparison  of  the  mean  of  the  PHC  of  velocity  at  hard- 
rock  stations  and  the  predicted  values. 

Figure  4  shows  the  results  of  the  regression  analysis  for  the  PHC  of 
velocity  vs  slant  distance  for  alluvium  or  unknown  station  conditions. 

Figure  5  shows  the  results  of  the  regression  analysis  for  the  PHC  of 
acceleration  vs  slant  distance  for  hard- rock  stations. 

Figure  6  shows  a  comparison  of  the  predicted  PHC  of  acceleration  at  hard- 
rock  stations  and  the  observed  mean. 

Figure  7  shows  the  results  of  the  regression  analysis  for  the  PHC  of,  accel¬ 
eration  vs  slant  distance  for  alluvium  or  unknown  station  conditions. 

The  formulae  for  the  regression  curves  based  on  the  Rio  Blanco  data 
are; 


1.  PHC  of  Velocity  vs  Slant  Distance 

Hard- rock  Media  1 1  to  100  km 

V=1352  R-2.  103  lo'=1.446 

2.  PHC  of  Velocity  vs  Slant  Distance 

Alluvium  or  Unknown  Media  17  to  150  km 
V=7039  R"^'^^  l(r-1.76 
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Table  2,  Regression  data. 


Station 

Hard  Rock 

RB-3 

RB-5 

RB-10 

RB-11 

RB-12 

RB-13 

RB-14 

RB-15 

RB-16 

$ 

RB-20 

RB-23 

Unknown 

Arco  #10 

RB-4 

RB-8 

RB-6 

RB-7 

RB-9 

NSG~1 

NSG-2 

NSM-1 

NSM-2 

RB-26E 

RB-26W 

NSA-2 

NSA-1 

NSD-2 

NSD-1 

NSR-2 

RB-17 

NSR-1 

RB-18 

RB-19 

RB-21 

RB-22 

RB-24 


Slant 

PHC 

PHC 

listance 

Vel. 

Acc. 

(km) 

(cm/  sec) 

(g) 

11.  3 

4.  54 

,  0.296 

18.  9 

4.  05 

0.  162 

43.4 

0.  576 

'  0.  0254 

46.4 

0,  5.18 

0. 0106 

46.4 

0.  518 

0. 0140 

46.  9 

0.674 

0. 0240 

50.  2 

0.  371 

0. 0191 

52.  4 

0.  239 

0. 0093 

58.  1 

0.  382 

0.0137 

77.  1 

0. 0955 

0. 00338 

93.  1 

0. 0601 

0. 00368 

13.  0 

0.40 

17.  3 

9.50 

0.  156 

22.  8 

8.85 

0.434 

25.  3 

8.  90 

0.  227 

29.6 

5.  34 

0.  154 

35.9 

1.  35 

0.  038 

46.  1 

0.  933 

0.  024 

46.  5 

0.434 

0.  013 

47.6 

0.  700 

0. 0298 

48.  1 

0.  888 

0. 0326 

48.  1 

0.403 

0. 0146 

48.  1 

0.  305 

0. 0146 

49.  5 

0.  984 

0. 0303 

50.  2 

0.  480 

0. 0203 

52.  5 

0.  250 

0. 0088 

52.  7 

0.280 

0. 0126 

57.  2 

0.  744 

0. 0192 

57.  3 

0.391 

0. 0206 

58.4 
64.  1 
71.  5 
82.  5 
106. 1 
148.  9 


0.745 
0.624 
0.  382 
0.  316 
0.  132 
0.  156 


0.  0190 
0. 0185 
0. 0183 
0.  0106 
0. 00496 
0. 00321 
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PHC  VELOCITY  (cm/sec^ 


-  OBSERVED  MEAN 

ALLUVIUM  OR  UNKNOWN  STATION  CONDITIONS  - ±\cr 


Figure  4.  PHC  of  velocity  vs  slant  distance, 
unknown  or  alluvium  stations. 
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Figure  5, 


PHC  of  acceleration  vs 
hard-rock  stations. 
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slant  distance. 
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PHC  ACCELERATION  g't 


* 


3> 


HARO  ROCK  STATIONS 


Figure  6.  Predicted  and  observed  PHC-of-acceleration 
values  for  hard- rock  stations,. 
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Figure  7.  PHC  of  acceleration  vs  slant  distance, 
unknown  or  alluvium  stations. 
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PHC  of  Acceleration  vs  Slant  Distance 


Hard- rock  Media 


11  to  100  km 


A^88.  7  R-2.  238 


l(r-  1.37 


4,  PHC  of  Acceleration  vs  Slant  Distance 

Alluvium  or  Unknown  Media  17  to  150  km 
A^135.  9  R-2*  204  lo-^l.  b39 


In  the  above  equations,  R  is  in  kilometers,  V  is  in  cm/sec,  and  A  is  in  g's. 

Table  3  shows  a  comparison  of  the  station  response  factor  for  stations  that 
were  in  the  same  location  for  both  Project  Rulison  and  Project  Rio  Blanco. 

The  calculated  means  for  the  Rulison  stations  were  computed  using  the 
Rulison  regression  equations(2)  at  the  appropriate  slant  distance.  The 
variation  is  the  observed  value  divided  by  the  calculated  value.  The  same 
analysis  applies  to  the  Rio  Blanco  data.  The  comparison  is  somewhat 
inconclusive,  although,  in  general,  if  the  station  response  was  high  (i,  e.  , 
a  variation  greater  than  1)  for  Rulison  it  was  also  high  for  Rio  Blanco,  and 
if  the  station  response  for  Rulison  was  low  (i.  e.  ,  a  variation  less  than  1) 
it  was  also  low  for  Rio  Blanco,  Exceptions  are  the  Rifle  Church  and  Harvey 
Gap  Dam. 
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Table  3.  Comparison  of  Rulison  and  Rio  Blanco  data 
from  the  same  recording  locations. 


Observed 


Station 

No. 

Dist. 

(km) 

Compo¬ 

nent 

Velocity 
(cm/ sec) 

Calculated 

Mean 

R-25 

10.6 

Z 

7.  34 

R 

8.  14 

T 

3.  30 

FHC 

8.  14 

6.24 

RB-12 

46.4 

Z 

0.  296 

R 

0.  518 

T 

0.  466 

PHC 

0.  518 

0.461 

R-08 

14.  6 

Z 

0.  830 

R 

0.  613 

T  , 

0.  665 

PHC 

0.  665 

3.  91 

RB-11 

45.7 

Z 

0.294 

R 

0.  342 

T 

0.  235 

• 

PHC 

0.  342 

0.492 

R-26 

20.2 

Z 

2.  88 

R 

1.89 

T 

2.  19 

PHC 

2.  18 

2.77 

RB-16 

58.  1 

Z 

0.276 

R 

0.  382 

T 

0.253 

PHC 

0.  382 

0.261 

Observed/ 

Calculated 

Variation 

Observed 
Acceler¬ 
ation  (g) 

Calculated 

Acceleration 

Observed/ 

Calculated 

Variation 

1.3 

0.  531 

0.  358 

0.  136 

0.  358 

0.268 

1.3 

1.  1 

0. 0105 

0. 0140 

0. 0140 

0. 0140 

0. 0291 

2.  1 

0. 0495 

' 

0.  17 

0. 0275 

0. 0287 

0. 0287 

0.  166 

0.  17 

0.70 

0. 0075 

0. 0125 

0. 0067 

0. 0125 

0. 0301 

0.41 

0.80 

0.  090 
0.0565  - 

0.  076 

0.  076 

0.  106 

0.71 

1.45 

0. 0078 

0. 0137 

0. 0086 

0. 0137 

.  0.0177 

1.  30 

\ 


t 


Table  3.  Comparison  of  Rulison  and  Rio  Blanco  data  from 
the  same  recording  locations.  (Cont. ) 


Station 

Location 


Debeque 

School 


Colony 

Tower 

CO 

I 

CO 


Harvey 
Gap  Dam 


Craig 


Station 

No. 

Dist. 

(km) 

Compo  - 
nent 

Observed 
Velocity 
(cm/ sec) 

Calculated 

Mean 

R-12 

22.8 

Z 

2.  07 

R 

1.34 

T 

1.73 

PHC 

1.73 

2.  13 

RB-15 

52.4 

Z 

0.  331 

R 

0.203 

T 

0.  239 

PHC 

0.239 

0.343 

16 

29 

Z 

0.904 

Rulison 

R 

7.  36 

T 

3.  36 

PHC 

7.  36 

1.26 

RB-7 

29.6 

Z 

5.08 

R 

4.  08 

T 

5.  34 

PHC 

5.  34 

1.55 

R-17 

32.4 

T(East 

1.07 

1.00 

end) 

RB-18 

64.,  1 

T(East 

0.  129 

0.201 

end) 

R-38 

127 

Z 

'0.0396 

R 

0.0437 

T 

0.0524 

PHC 

0.0524 

0.051 

Observed/ 

Calculated 

Variation 

Observed 
Acceler¬ 
ation  (g) 

Calculated 

Acceleration 

Observed/ 

Calculated 

Variation 

0.79 

0. 0965 

0. 0532 

0. 0473 
0.0532  . 

0.  080 

0.67 

0.70 

0. 0103 

0. 0093 

0. 0084 
0,0093 

0.022 

0.42 

NA« 

NA* 

NA* 

4.7 

NA* 

3.4 

1.07 

0.0353 

0.037 

1.06 

0.  65 

NA* 

0.014 

1.03 

0.000785 

0.00175 

0.00137 

0. 00175 

0.00155 

1.1 

♦Not  available 


I 


AS-74 


Station 

Location 

Craig 
(Cont,  ) 


Aspen 


Table  3.  Comparison  of  Rulison  and  Rio  Blanco  data  from 
the  same  recording  locations.  (Cont.) 


Station 

No. 

Dist. 

(km) 

Compo¬ 

nent 

Observed 

Velocity 

(cm/sec) 

Calculated 

Mean 

RB-22 

106 

Z 

0.  068 

R 

0.  124 

T 

0.  134 

PHC 

0.  134 

0.053 

R-30 

100 

Z 

0.  150 

R 

0.  193 

• 

T 

0.  244 

PHC 

0.  244 

0.  086 

RB-24 

149 

Z 

0.  117 

R 

0.  156 

T 

0.  115 

PHC 

0.  156 

0.0217 

Observed/ 

Calculated 

Variation 

Observed 
Acceler¬ 
ation  (g) 

Calculated 

Acceleration 

Observed/ 

Calculated 

Variation 

2.  53 

0. 00175 

0. 00496 
0.00483 

0. 00496 

0. 0047 

1.05 

3.  3 

0. 00428 

0. 00656 

0. 00594 

0. 00656 

0.00269 

2.4 

7.  1 

0. 00343 

0. 00321 
0.00284 

0. 00321 

0. 00223 

1.45 
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3.  SEISMICITY  AND  GROUND  MOTION 


3.1  SEISMICITY 

An  intensive  program  was  initiated  to  determine  whether  or  not  Project 
Rio  Blanco  caused  any  change  in  the  natural  seismicity  of  the  area  within 
40  km  of  the  EW. 

Geologic  and  seismic  surveys  showed  that  there  were  no  surface  or  subsur¬ 
face  indications  of  faulting  within  12,000  ft  of  the  EW.  A  seismic  station 
was  installed  November  20,  1971,  some  20,  000  ft  NNW  of  the  EW.  It 
monitored  the  seismicity  of  the  area  within  40  km  continuously  up  to 
December  17,  1973.  In  the  18-month  predetonation  period,  only  1  micro¬ 
earthquake  was  recorded.  At  detonation  time,  this  station  became  part 
of  a  six-station  array  deployed  to  monitor  the  seismic  activity  during  and 
after  the  detonation.  This  array  was  discontinued  19  days  after  the  detona¬ 
tion  and  the  single  station  was  returned  to  its  original  configuration  and 
continued  monitoring  the  seismicity  of  the  area  until  December  17,  1973. 

In  the  7 -month  period,  commencing  19  days  postdetonation,  only  1  micro¬ 
earthquake  was  recorded  within  40  km  of  the  station. 

A  6-instrument  array  was  deployed  17  days  before  detonation  to  record  and 
locate  the  hypocenters  of  the  subsequent  aftershocks.  The  location  of  these 
stations  is  shown  in  Figure  8.  One  station  was  located  900  ft  north  of  the 
^  EW  and  the  other  5  in  a  circular  array  at  distances  ranging  from  16,  000  ft 

to  20,  000  ft. 

In  the  period  from  D-19  days  to  detonation,  no  micro- earthquakes  were, 
recorded  in  the  area.  From  detonation  to  H+5  minutes,  the  records  were 
too  noisy  to  permit  identification  of  events.  From  H+5  minutes  to  D+8 
days,  120  aftershocks  were  recorded  on  at  least  4  of  the  6  seismic  stations. 
The  time  distribution  of  the  aftershocks  is  shown  in  Figures  9  and  10.  The 
largest  aftershock  had  a  magnitude  of  2.4  and  occurred  82  minutes  after 
the  detonation.  The  locations  of  the  hypocenters  were  computed  using  the 
computer  code  EPIC  originally  written  by  Turcotte  &  Bolt  of  the  University 
of  California.  The  epicenter  locations  are  shown  in  Figure  11.  The  great¬ 
est  epicentral  distance  is  2,  900  ft  in  a  northerly  direction.  There  is  a 
tendency  for  the  epicentral  pattern  to  be  elongated  in  a  northwest- southeast 

direction. 

The  depth  distribution  of  the  hypocenters  is  shown  in  Figure  12.  It  is 
apparent  that  the  majority  of  the  aftershocks  occur  in  the  vicinity  of  the 
explosives  and  are  chimney  related.  In  no  case  did  an  aftershock  occur 
near  or  in  the  oil  shale  horizon. 
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Figure  8.  Rio  Blanco  aftershock  hypocenter  location  stations. 
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Figure  9.  Number  of  aftershocks  per  5-minute  interval. 
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Figure  10. 


Nvmiber  of  aftershocks  per  l-hour  interval  after  detonation. 
(No  activity  detected  after  0629  GMT,  5/25/73. ) 
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Figure  11.  Location  of  epicenters. 
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Figure  12^.  Depth  distribution  of  hypocenters.- 
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On  the  basis  of  data  reported  by  19  stations  of  the  worldwide  network  to  the 
National  Earthquake  Information  Center,  Project  Rio  Blanco  had  an  average 
magnitude  of  5.  4  m^. 

A  second  seismic  array  was  fielded  by  the  U.S.  Geological  Survey/ National 
Center  for  Earthquake  Research  (USGS/NCER)  to  monitor  cavity  collapse 
and  aftershocks.  The  USGS/NCER  emplaced  33  seismometers  within  5.  8 
km  of  the  EW,  Twenty-five  stations  operated  successfully.  To  date,  the 
results  from  this  program  have  not  been  published;  however,  it  is  under¬ 
stood  that  there  is  close  agreement  with  the  results  of  the  CER  program. 

As  a  result  of  the  seismicity  monitoring  program  it  is  evident  that: 

1,  Project  Rio  Blanco  had  no  long  range  effect  on  the  natural  seismicity 
of  the  area. 


2.  The  aftershocks  appear  to  be  chimney  related. 

3.  No  aftershocks  occurred  in  or  near  the  oil  shale  horizon. 


3.  2  CORRELATION  OF  STRUCTURAL  DAMAGE  TO  LOW 

LEVEL  GROUND  MOTION 


An  experimental  program  to  investigate  the  correlation  of  structural  damage 
to  low  level  ground  motion  was  carried  out  by  JABAE  with  funding  from  the 
National  Science  Foundation  (NSF),  The  ground  motion  data  are  shown  in 
Table  1  and  are  from  those  stations  with  a  NS  prefix.  The  results  of  the 
correlation  are  not  available  as  yet>  although  publication  is  expected 
shortly. 


3.  3  STRUCTURAL  RESPONSE  TO  GROUND  MOT  ION 

t 

An  experimental  program  to  monitor  the  response  of  five  selected  structures 
to  ground  motion  was  carried  out  by  KMA  with  funding  from  NSF.  The 
structures  monitored  were  the  White  River  Electric  microwave  tower,  the 
teacher's  residence  at  the  Rock  School,  the  Chadbourne  depropanizer 
tower,  the  Rangely  College  gymnasium,  and  the  Harvey  Gap  Dam. 

The  ground  motion  data  at  these  structures  are  included  in  Table  1.  The 
objective  of  the  program  was  to  compare  the  response  of  the  structure  as 
predicted  by  predetonation  structural  modeling  with  the  observed  response. 
The  results  of  this  comparison  have  not  been  published  at  this  time. 
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3.4 


GROUND  SPALL  DOCUMENTATION  PROGRAM 


This  was  a  program  designed  and  executed  by  the  Lawrence  Livermore 
Laboratory.  The  instrumentation  and  results  are  described  in  reports 
published  by  the  Lawrence  Livermore  Laboratory^^" ^ The  conclu¬ 
sions  reached  are  that  the  measured  spall  occurrences  are  in  excellent 
agreement  with  the  predetonation  predictions namely  that  the  spall 
zone  can  be  described  as  a  shallow  dish  with  a  depth  of  less  350  ft  and 
a  radius  of  less  than  24,  000  ft.  Therefore,  since  the  oil  shale  exists 
from  depths  of  about  400  ft  to  2,  300  ft,  no  spall  has  occurred  in  the  oil 
shale  deposits. 


3.5 


SANDIA  LABORATORIES  PORE  PRESSURE  AND  GROUND 
MOTION  MEASUREMENTS^^^) 


Sandia  Laboratories,  operating  with  funds  obtained  from  the  National  Science 
Foundation,  measured  changes  of  water  pore  pressure  in  saturated  soils  of 
the  Fawn  Creek  Valley.  Changes  in  water  pore  pressure  are  a  good  index  of 
the  small  amount  of  consolidation  which  precedes  soil  liquefaction.  The 
ground  motion  waves  developed  by  Project  Rio  Blanco  impose  a  seismic 
input  similar  to  an  earthquake,  so  Project  Rio  Blanco  provided  an  oppor¬ 
tunity  to  study  pore  pressure  changes  in  undisturbed  soils. 


•  V 


-v' . 

'  •:  > . 


Three  effects  worthy  of  note  were  observed.  The  first  was  that  significant 
pressure  increases  were  observed  before  spall  closure  at  the  two  closer 
stations.  The  second  and  more  expected  trend  was  that  measurements  in 
less  permeable  material  showed  smaller  pressure  oscillation  than  those  in 
more  permeable  media  even  though  the  sustained  pressure  increases  were 
about  the  same.  Finally,  inspection  of  the  records  reveals  there  is  the 
expected  close  correlation  between  the  transient  values  of  pore  pressure 
and  vertical  acceleration.  It  is  also  apparent  that  radial  accelerations 
are  affecting  pore  pressures  in  the  later  portion  of  the  traces.  In  general, 
the  sustained  pore  pressure  increases  were  much  lower  than  would  be 
required  to  achieve  liquefaction.  We  found  no  evidence  of  liquefaction  in 
undisturbed  soils  in  the  postdetonation  search. 


'•  I  , 

3.6  STRUCTURAL  BRACING  . 

The  Project  Rio  Blanco  responsibilities  of  CER's  structural  response  •  ,•  ’ 

consultant  (KMA)  included  an  inventory  of  all  structures  in  the  project 
area,  structural  response  predictions,  slope  stability  evaluations,  esti- 
mates  of  damage,  safety  measure  recommendations,  and  the  bracing  of 
certain  structures.  All  but  the  bracing  and  inventory  have  been  covered 
in  previous  KMA  Project  Rio  Blanco  reports(^2» 
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It  should  be  noted  that  the  structural  bracing  program  for  Pr,oject  Rio 
Blanco  was  fairly  extensive.  The  most  time-consuming  effort  was  con¬ 
ducted  at  the  historic  Rock  School,  which  was  completely  regrouted  to 
a  depth  of  about  6  inches  from  the  outside  and  heavily  braced.  Other 
bracing  was  installed  at  the  adjacent  concrete  block  school,  the  Shell 
Oil  Shale  Plant,  Equity  Oil  Shale  Plant,  Bookcliff  Bowling  Alley,  Chad- 
bourne  Gas  Compressor  Plant,  El  Paso  Gas  Plant,  Rangely  College, 

Harvey  Gap  Dam  spillway,  as  well  as  a  number  of  trailers,  chimneys, 
and  miscellaneous  items. 

The  observed  ground  motion  mean  appeared  less  than  ground  motion  pre¬ 
dictions,  which  were  2,  5  times  the  CER  median  predictions,  used  to 
evaluate  seismic  damage.  In  view  of  the  lower  level  of  ground  motion, 
it  is  somewhat  difficult  to  evaluate  the  effects  of  the  bracing.  It  was 
visually  obvious  that  the  bracing  at  the  Equity  Oil  Shale  Plant  was  instru¬ 
mental  in  precluding  damage.  Attachment  bolts,  soil  anchors,  and  wood 
bracing  showed  the  effects  of  absorbing  dynamic  forces.  There  v/ere  also 
minor  visual  indications  at  the  Rock  School,  However,  the  bracing  was 
designed  to  absorb  only  the  dynamic  force  in  excess  of  that  which  the 
structure  could  withstand.  For  example,  tank  bracing  was  primarily 
placed  to  carry  the  dynamic  force  not  absorbed  by  soil-structure  friction. 
Low-level  ground  motions  may  not  have  caused  such  bracing  to  come 
into  play,  but  visual  observations  usually  would  not  be  adequate  for  assess¬ 
ing  the  situation  except  in  extreme  cases.  In  retrospect,  the  only  bracing 
that  might  have  been  eliminated  was  that  placed  at  the  Bookcliff  Bowling 
Alley  and  Harvey  Gap  Dam  spillway.  Both  were  sources  of  concern 
primarily  because  of  the  investigations  made  for  Project  Rulison.  The 
bracing  installed  was  in  accordance  with  the  safety  factors  added  to  the 
‘  mean  ground  motion  predictions.  As  the  latter  are  refined,  bracing 
requirements  can,  of  course,  be  modified  accordingly. 

As  part  of  a  structure  bracing  program  for  Project  Rio  Blanco,  AEG 
employed  the  services  of  JABAE  to  conduct  structure  response  analyses 
and  to  design  structural  bracing  for  the  190-foot-tall  Colony  Development 
Oil  Shale  Retort  Tower,  situated  on  Parachute  Creek  near  Grand  Valley, 
Colorado(l^).  Other  Colony  Development  structures  such  as  office 
trailers,  sheds,  etc. ,  were  evaluated  by  the  JABAE  resident  engineer  at 
the  tower  site  and  conservative  safety  measures  were  instituted  prior  to 
the  detonation. 

The  structural  analysis  was  based  on  the  plus-one  sigma  level  of  ground 
motion  as  predicted  by  ERC^^"^)  for  that  specific  location.  A  two-dimensional 
computer  analysis  was  utilized  by  JABAE  to  develop  the  modified  bracing 
scheme.  The  general  bracing  scheme,  including  sizes  of  steel  members, 
was  developed  by  JABAE.  Fabrication  drawings  and  connection  details 
for  all  new  and  existing  bracing  members  were  developed  in  the  field  by 
the  JABAE  resident  engineer  who  was  also  responsible  for  inspection. 
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Installation  of  bracing  on  the  tower  was  performed  by  the  Stearns- Roger 
Company  of  Denver,  Colorado,  under  contract  to  CER  during  the  period 
April  30  to  May  14,  1973.  The  work  consisted  of  fabricating  and  placing 
approximately  21  tons  of  structural  steel  shapes  (134  pieces)  and  750 
pounds  of  welding  electrode  to  strengthen  the  existing  tower  bracing 
system. 

The  site  was  evacuated  at  detonation  time  as  a  safety  measure  and  all  plant 
personnel  moved  to  a  safe  location  some  distance  down  the  canyon  from  the 
tower  site.  Ground  motion  was  perceptible  at  this  location  at  detonation 
time. 

No  damage  or  indications  of  adverse  effects  from  seismic  motion  were 
observed  during  an  immediate  postdetonation  inspection  of  the  tower 
structure  norat  any  other  facilities  at  the  tower  site  (except  some  minor 
rock  falls  observed  in  the  mine).  The  tower  response  appears  to  have 
been  at  lower  levels  than  those  used  to  develop  the  bracing  design. 

Minor  rock  falls  were  also  observed  on  the  roadway  up  to  the  Colony 
Development  guard  gate  and  on  the  roadway  above  the  tower  structure  to 
the  mine  entrance,  but  these  were  readily  cleared.  Near  a  waterfall  about 
800  ft  north  of  the  tower,  a  rock  about  1  cubic  foot  in  size  was  dislodged 
from  the  overhanging  cliff^^^h 

3.7  BRIDGES 

Although  KMA  predicted  no  damage  to  bridges  in  the  Project  Rio  Blanco 
area,  six  were  selected  for  documentation  by  JABAE,  the  AEC’s  structural 
response  contractor.  The  bridges  selected  for  study  in  this  program 
include: 

Steel  truss  bridge  crossing  Colorado  River  at  Grand  Valley  (47  km  from 


EW) 


180-foot  steel  bridge  owned  by  Colorado  Highway  Department  located 
at  White  River- -Highway  64  crossing  (34.  74  km  from  EW) 

90-foot  steel  truss  bridge  owned  by  the  Rio  Blanco  County  located  at 
White  River--Ivo  Shults  crossing  (33.  84  km  from  EW) 

69-foot  timber  bridge  owned  by  the  Colorado  State  Highway  Department 
located  at  Highway  64  and  Smith  Gulch  (33.  95  km  from  EW) 

108-foot  concrete  bridge  owned  by  Rio  Blanco  County  located  at  White 
River--County  Road  (38.  54  km  from  EW) 

80-foot  steel  bridge  owned  by  Rio  Blanco  County  located  at  White  River- - 
Piceance  road  crossing  (35.64  km  from  EW) 
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KMA  inspected  these  bridges  carefully  before  and  after  the  detonation  and 
found  no  visible  evidence  of  detonation- related  damage. 
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3.  8 


OIL  AND  GAS  WELLS 


Forty-two  oil  and  gas  wells  at  distances  between  1. 5  and  7  miles  from  the 
EW  were  inspected  and  photographed  before  and  after  the  detonation.  No 
surface  damage  to  any  of  the  wells  or  ancillary  equipment  was  found.  Some 
soil  movement  was  noted  around  the  well  casings  in  some  of  the  wells  within 
4  miles  of  the  EW. 


ROCK  FALLS 


3.9 


Rock  falls  causing  road  obstructions  occurred,  but  those  of  any  consequence 
were  limited  to  locations  in  the  Piceance  Basin.  Rock  falls  were  reported 
on  nine  locations  along  the  Piceance  Creek  Road,  eight  along  Black  Sulphur 
Creek  road,  eight  along  Fawn  Creek  road,  four  on  Hunter  Creek  road,  and 
one  bn  East  Stuart  Creek  road. 

No  rocks  reached  the  county  road  on  Parachute  Creek;  however,  small 
rocks  were  scattered  at  four  locations  between  the  West  Fork  Bridge  and 
the  Colony  Oil  Shale  site.  On  Cathedral  Creek,  three  to  four  rocks  drop¬ 
ped  on  the  road  and  some  rocks  fell  on  the  Douglas  Creek  Road  just  south 

of  Little  Horse  Draw, 

No  known  rock  falls  occurred  in  De  Beque  or  Glenwood  canyons,  the  Douglas 
Pass,  or  on  any  of  the  other  State  or  Federal  highways,  except  for  a  few 
large  boulders  on  the  west-bound  lane  of  Highway  64,  about  12  miles  west 
of  Meeker.  This  rockfall  occurred  within  one  of  the  exclusion  areas  con¬ 
trolled  by  roadblocks.  The  locations  of  the  .reported  rock  falls  are  shown 

in  Figure  13. 


3.10  WELLS  AND  SPRINGS 

A  pre-  and  postdetonation  inspection  was  made  of  all  springs  and  water  w 
within  10  miles  of  the  EW.  Water  samples  analyses  show  no  significant 
difference  between  the  chemical  compositions  of  the  water  samples  taken 
before  and  after  the  detonation.  No  increase  in  turbidity  was  detected. 

♦ 

An  increase  in  the  yield  of  springs  as  far  as  10  miles  from  the  EW  was 
observed  after  the  detonation.  To  date,  the  yield,  while  less  than  that 
observed  immediately  after  the  detonation,  has  still  not  returned  to  the 
predetonation  value. 
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In  the  Fawn  Creek  and  Hunter  Creek  valleys,  new  springs  have  been  noted 
and  old  spring  locations,  dry  for  decades,  appear  to  have  been  revived  by 
the  shock  wave.  The  damage  to  wells  in  the  area  out  to  10  miles  from  the 
EW  appear  to  be  mechanical  in  nature.  A  few  springs  over  the  divide  in 
the  Colorado  River  drainage  basin  appear  to  have  dried  up. 


3.  11  MICROTREMOR  SURVEYS 

Microtremor  surveys  were  made  in  1970  and  1971  of  the  towns  of  Rifle, 
Grand  Valley,  Meeker,  Rangely,  and  De  Beque  (Figures  14  through  18). 

A  comparison  of  the  surveys  and  damage  patterns  in  Rifle  and  Grand 
Valley  from  Project  Rulison  indicates  that  for  those  locations  where 
the  ground  resonance  is  less  than  12.  5  Hz,,  the  damage  incidence  to  resi¬ 
dential  structures  is  double  that  of  locations  where  the  ground  resonance 
is  above  12.  5  Hz(2,  1*7)^ 

Project  Rio  Blanco  did  not  cause  enough  damage  in  Rifle,  Grand  Valley, 
De  Beque,  and  Rangely  to  establish  a  damage  pattern.  The  damage  loca¬ 
tions  in  these  towns  in  relation  to  the  microtremor  surveys  are  shown  in 
Figures  14  through  17.  In  the  town  of  Meeker,  there  were  36  damage 
claims  paid  as  of  December  1973.  The  location  of  these  claims  in 
relation  to  the  microtremor  survey  is  shown  in  Figures  18  and  19.  An 
analysis  shows  that  there  are  approximately  twice  as  many  damage 
locations  where  the  ground  resonance  is  less  than  12,  5  Hz  than  where 
the  ground  resonance  is  greater  than  12.5  Hz. 


3. 12  EL  PASO  NATURAL  GAS  COMPANY  PROGRAM 

The  primary  objective  of  this  work  is  to  further  improve  the  accuracy  of 
effects  predictions  of  underground-nuclear-explosion  induced  ground  motions 
upon  man-made  and  natural  structures  and  the  bioenvironment  near  the 
detonation  point.  A  number  of  structures  exist  near  the  Rio  Blanco  EW 
which  are  similar  to  structures  near  the  proposed  Wagon  Wheel  Project., 
Project  Rio  Blanco  provides  an  opportunity  to  improve  capabilities  for 
evaluations  relative  to  Project  Wagon  Wheel, 


The  following  specific  objectives  will  be  supportive  of  the  primary  objective: 

1.  Determine  the  increase  and  decay  of  pore-water  pressure  in 

saturated,  loosely  consolidated  sands  over  several  hours  dura¬ 
tion  to  determine  the  potential  of  these  sands  for  liquefaction 
due  to  high  amplitude,  short  duration  pulses  of  a  limited  number 
of  cycles  characteristic  of  nuclear  detonations. 

The  measurements  showed  the  changes  in  average  pore  pressure 
were  positive  and  died  away  in  times  of  tens  of  seconds. 
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Figure  14. 


Rifle  microtremor  survey. 
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Figure  17,  Rangely  microtremor  survey. 
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Figure  18.  Meeker  microtremor  survey. 
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Figure  19-  Meeker  microtremor  results  vs  damage  distribution. 
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2. 


Determine  the  potential  for  the  occurrence  of  permanent  and 
transient  deformations  with  regard  to  small,  single-span 
bridges  typical  of  those  associated  with  ranching  operations 
and  rural  transportation. 

Engineering  surveys  of  all  bridge  dimensions  and  positions  were 
obtained  pre-  and  postdetonation.  Preliminary  results  indicate 
no  permanent  deformations  in  excess  of  the  survey  accuracy  of 
0,01  foot  occurred  to  any  bridge  components  with  the  exception 
of  minor  cracks  in  the  abutment  fill  of  one  bridge. 

3.  Determine  the  potential  for  rock  falls  and  minor  sloughage  from 
slopes  of  comparable  material  and  geometry  as  those  found  near 
Project  Wagon  Wheel. 

Several  bluffs  and  slopes,  including  steep  ditch  banks,  were 

selected  for  pre-  and  postdetonation  engineering  surveys  and 

photographs.  No  massive  slides  or  slope  movements  occurred 

on  any  of  the  bluffs  and  slopes  selected  for  the  EPNG  study.  J  ' 

Most  falls  consisted  of  individual  rocks  or  groups  of  rocks. 

4.  Determine  if  wave  lengths  in  ground  motion  exist  such  that  large 
strains  can  be  induced  in  structures  of  a  few  hundred  feet  in 
length  by  the  ground  motion  from  the  detonation. 

The  analysis  indicates  the  presence  of  at  least  two  distinct  arrivals 
in  the  high  amplitude,  early  part  of  the  record.  The  earlier  arrival  * 

'  '  *  5  Jr ' 

appears  to  exhibit  frequencies  of  about  6  Hz,  velocities  of  about 
8,200  feet  per  second,  and  wave  lengths  of  about  1, 325  feet.  The  !;  > 
later  arrivals  show  frequencies  of  about  4  Hz  with  corresponding  ■' 

velocities  of  about  4,  830  feet  per  second  and  wave  lengths  of  about  ' 

1,175  feet. 

5.  Investigate  the  effects  upon  mechanical  structures  characteristic  ^  - 

of  conventional  gas  field  operations, 

Accelerations  and  velocities  were  successfully  measured  at  several  ^ 

points  in  and  on  equipment  at  EPNG’s  gas  compression  plant  at  about  . 

18.  8  km.  Engineering  surveys  were  made  of  dimensions  and  positions  ’•* 
throughout  the  facility.  No  deformations  in  excess  of  the  survey 
accuracy  of  0.  01  foot  were  found, 

6.  Investigate  the  effects  upon  vegetation. 

Six  20- square-foot  vegetation  area's  were  established  and  quanti-  - 

tatively  and  photographically  documented  predetonation.  Four  of 
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these  locations  were  within  1, 000  feet  of  the  EW,  with  the  other 
two  locations  being  adjacent  to  EPNG  study  bluffs  and  slopes  at 
8,  0  and  11,  26  km  from  the  EW.  From  early  observations,  the 
only  significant  effect  of  ground  motion  was  that  relating  to 
minimal  sloughage  of  steep  banks  and  loss  of  plants  existing 
very  near  those  banks.  On  the  close-in  plots  and  adjacent 
relatively  level  surfaces,  some  soil  disturbance  was  observable 
but  was  not  reflected  by  plant  wilting  or  cessation  of  growth. 

The  early  results  of  this  program  were  reported  in  PNE- RB-48^^).  No 

further  information  is  available  at  this  time. 


3.  1  3  MINES  AND  QUARRIES 

A  pre-  and  postdetonation  inspection  was  made  of  all  mines  and  quarries 
within  53  miles  of  the  EW.  Eight  active  mine  operations  in  this  area 
were  evacuated  at  detonation  time. 

Minor  rock  falls  were  reported  from  the  walls  in  the  Colony  Oil  Shale 
Mine.  None  of  the  other  inspected  mines  or  quarries  suffered  any  damage 
as  a  result  of  the  detonation. 

The  evacuation  procedures  instituted  for  the  eight  active  mine  operations 
in  the  area  resulted  in  production  delays  in  only  five  of  the  mines;  the 
NuGap  Coal  Mine  near  Silt,  and  the  three  gilsonite  companies  in  Bonanza, 
Utah,  stopped  work  for  2  to  4  hours,  while  the  Redwing  Coal  Mine,  near 
Axial,  was  closed  for  one  8-hour  shift.  The  combined  effect  was  a  delay 
in  mining  of  about  750  tons  of  coal  and  about  100  tons  of  gilsonite, 

3. 14  DAMAGE  CLAIMS  AND  ADJUSTMENTS 

As  of  October  17,  1974,  169  claims  have  been  reported  by  Continental  Oil 
Company.  Of  this  total,  114  have  been  settled  for  a  total  expenditure  of 
$52,648,  06.  Fifty-four  have  been  denied  or  are  reported  as  no  further 
action  to  be  taken.  One  is  pending.  All  claims  are  for  real  and  personal 
property  and  are  generally  of  a  superficial  nature,  as  predicted.  (Use 
of  the  term  "claim"  here  is  equivalent  to  the  use  of  the  v/ord  "complaint" 
in  similar  discussion  on  Project  Rulison). 

Table  4  shows  the  claim  distribution. 
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Table  4.  Property  damage  claims. 


Area 

Claims 

Settled 

Amount 

NFAC* 

Pending 

Piceance  Creek 

24 

19 

$20,  039.  18 

4 

1 

Roan  Creek 

5 

4 

$5, 532. 80 

1 

0 

Meeker  (Rural) 

10 

5 

$6,  508.  15 

5 

0 

Rifle  &  Silt 
(Rural) 

10 

7 

$2, 938. 70 

3 

0 

Grand  Junction 

23  ‘ 

16 

$1,842.  20 

7 

0 

Clifton 

4 

3 

$320.  00 

1 

0 

Grand  Valley 

6 

4 

$318.  00 

2 

0 

Meeker 

55 

38 

$11,  128.48 

17 

0 

Rangely 

6 

4 

$145.  00 

2 

0 

Rifle 

6 

4 

$612.  10 

2 

0 

Silt 

3 

1 

$150.  00 

2 

0 

•Other 

17 

9 

$2. 513. 70 

8 

_ 0 

TOTALS 

169 

114 

$52,648.  06 

54 

1 

*No  further  action 

contemplated. 

3.15  STRUCTURAL  INVENTORY 

A  structural  inventory  was  made  of  the  352  structures  located  v/ithin  30  miles 
of  the  EW.  The  total  estimated  fair  market  value  of  these  structures  is 
$22.  5  million  with  an  $18.  5  million  valuation  for  16  industrial  locations  and 
$4,  0  million  for  the  remaining  336  locations.  A  seismic  hazard  evaluation 
was  made  of  each  structure  and^  where  required,  a  recommendation  for 
bracing  was  made.  A  predetonation  estimate  was  made  for  seismic  damage 
amounting  to  $51,  300  t  $13,  200^^2).  The  actual  damage  as  quoted  in  the 
previous  paragraph  was  $52,648.  06,  some  of  which  was  non- structure  related 
(i.  e.  ,  irrigation  ditch  repair,  etc.  ).  This  is  a  remarkable  agreement  in  view 
of  the  state-of-the-art  for  seismic  damage  prediction. 


4.  CONCLUSIONS 


4.  1  GENERAL 

% 

The  predicted  impact  on  the  environment  resulting  from  Project  Rio 
Blanco  has  been  discussed  extensively  in  numerous  reports  including 
the  Environmental  Statement^ ^  and  Addendum,  and  public  meetings 
held  prior  to  the  execution  of  the  project.  A  comparison  of  predicted 
and  observed  effects  is  shown  in  Table  5, 


;4.  2  GROUND  MOTION 

The  ground  motion  predictions  made  by  CER  used  prediction  techniques 
based  on  the  results  of  the  Rulison  detonation.  Regression  analyses  of 
observed  velocity  data  for  hard  rock  stations  were  less  than  predetonation 
predictions. 

The  mean  of  the  observed  acceleration  data  were  somewhat  less  than 
predicted  but  not  as  much  less  as  the  velocity  mean  indicating  a  higher 
predominant  frequency  for  the  Rio  Blanco  ground  motion  than  for  Rulison# 


4. 3  AFTERSHOCKS 

On  the  basis  of  previous  detonation  experience,  it  was  expected  that  no 
more  than  a  few  dozen  aftershocks  would  occur,  none  having  an  equivalent 
Richter  magnitude  greater  than  3  (compared  with  equivalent  Richter  mag¬ 
nitude  5.  5- -body  wave-  -  for  the  explosion  itself).  Aftershocks  were  pre¬ 
dicted  to  cease  within  a  few  hours, 

* 

The’  aftershocks  generated  by  Project  Rio  Blanco  were  located  within  a 
few  thousand  feet  of  the  explosion  cavity.  There  were  relatively  few 
aftershocks.  Most  occurred  within  hours  and  all  within  8  days.  All 
aftershocks  detected  had  magnitudes  less  than  2,  5, 
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Table  5.  Comparison  of  effects  predictions  and  actual  effects. 


CATEGORY 
Residential  structures 


PREDICTED  EFFECTS 

Minor  architectural  damage 
($51,  300i  $13,000) 


ACTUAL  EFFECTS 

Minor  architectural 
damage  (about  $51,  000) 


Industrial  facilities 

Spall 


CO 

I 

VO 

VO 


OJ 
O 
C/) 

u 

Tl 

g 

ej 

tn  i 
ini 


C: 


Bridges 

Tunnels 

Power  and  telephone 
lines 

Communication  facili¬ 
ties 

Gas  &  oil  wells 
Mines 

Hydraulic  features 

Airports 


1 1 


No  damage  (except  possibly 
the  Equity  Plant) 


Minor  damage  at 
Equity  Camp 


The  results  of  the  ground  spall  documentation  program  indicate 
a  close  agreement  of  observed  occurrences  with  predicted 
occurrences  (PNE-RB-90) 


No  damage 
No  damage 
No  damage 


No  damage 
No  damage 
No  damage 


No  damage 


No  damage 


No  damage 
No  damage 


No  surface  damage 

Minor  rockfalls  from 
walls  of  Colony  Mine 


No  damage,  except  minor  Minor  damage  to  small 

'damage  to  small  irrigation  irrigation  ditches 

ditches  (sloughing  of  banks)  .  (sloughing  of  banks) 


No  damage 


No  damage 


i) 
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Table  5.  Comparison  of  effects  predictions  and  actual  effects.  (Cont. ) 


CATEGORY 


PREDICTED  EFFECTS 


ACTUAL  EFFECTS 


Railroads 


No  damage 


No  damage 


Roads 


Water  wells  and 
springs 


No  damage  (possible  tem¬ 
porary  blockage  due  to 
rockfalls) 

Slight  temporary  increase 
in  flow  and  turbidity 


CO 

I 

o 

o 


Temporary  blockage 
due  to  rockfalls 


Larger  increases  in 
flow  rates  than  an¬ 
ticipated  and  to  a 
greater  areal  extent. 
In  many  instances, 
flow  rates  have  not 
returned  to  predeto¬ 
nation  levels.  Appre¬ 
ciable  but  temporary 
increases  in  sediment 
levels  to  Fawn  Creek 


1 
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SEISMICITY  IN  THE  RANGELY,  COLORADO,  AREA;  1962-1970* 

Bv  Jamhs  F.  Gibbs,  John  H.  Hcaly,  C.  Barry  Ralligh,  and  John  Coaklly 

ABSTRACT 

Soisnilc  data  recorded  for  a  7-ycar  period  at  the  Uinta  Basin  Observatory  were 
searched  for  earthquakes  ori'jinating  near  an  oil  field  at  Rangely,  Colorado,  65  km 
ESE  of  the  observatory.  Changes  in  the  number  of  earthquakes  recorded  per  year 
appear  to  correlate  nith  changes  in  the  quantity  of  fluid  injected  per  year.  Between 
November  1962  and  January  1970,  976  earthquakes  were  detected  near  the  oil 
field  by  the  UBO  station;  320  earthquakes  were  larger  than  magnitude  1. 

Richter  local  magnitudes  are  estimated  from  both  5-wave  and  P-3vave  measure¬ 
ments;  a  method  based  on  the  duration  of  the  seismic  signal  is  used  to  estimate  the 
magnitude  of  the  larger  shocks.  Magnitude  of  the  two  largest  shocks  was  3.4  and 
3.3.  The  total  seismic  energy  released  was  10'^  ergs.  During  this  same  period,  the 
energy  used  for  water  injection,  measured  at  the  wellhead,  was  10^’  ergs. 

Introduction 

Interest  in  the  triggering  of  earthquakes  by  injection  of  fluid  into  the  crust  of  the 
Earth  has  increased  since  Evans  (1966)  found  a  correspondence  between  injection  of 
chemical  waste  material  into  a  deep  disposal  w'ell  near  Denver,  Colorado,  and  earth¬ 
quake  activity  in  the  Denver  area  (see  also  Healy  et  al.,  1968).  At  least  one  study  is  now 
underway  in  an  area  w'here  earthquakes  and  injection  are  found  together  to  determine  if 
man  is,  in  fact,  inadvertently  triggering  earthquakes.  This  study  is  being  conducted  at 
Rangely,  Colorado,  where  the  U.S.  Geological  Survey  in  cooperation  wiUi  Chevron  Oil 
Company  is  attempting  to  control  the  seismic  activity  by  changing  the  fluid  pressure  in  a 
small  fault  zone.  Prior  to  this  detailed  study,  there  were  no  seismographs  located  in  the 
Rangely  area.  The  nearest  instruments  were  those  of  the  Uinta  Basin  Seismological 
Observatory  located  65  km  to  the  WNW.  In  this  paper,  we  examine  the  seismic  data 

available  before  the  experiment  was  begun. 

Seismograph  records  for  the  period  November  1962  to  January  1970  recorded  at  the 
Uinta  Basin  Seismological  Observatory  (UBO),  were  searched  for  earthquakes  originating 
near  the  oil  field  and  the  number  of  earthquakes  checked  against  the  quantity  of  water 
iniected  to  see  if  a  correlation  could  be  established.  For  6  of  the  7  years  studied,  changes  in 
the  amount  of  water  injected  per  year  have  a  general  correspondence  to  changes  in  the 
number  of  earthquakes  per  year. 

Estimates  of  earthquake  magnitude  were  calculated  from  5-wave,  P-wave,  and  signal- 
duration  measurements  in  order  to  provide  a  magnitude  cutoff  below  which  earthquakes 
would  not  be  counted  and  to  compare  magnitudes  with  new  data  being  recorded  by  the 
U.S.  Geological  Survey.  The  magnitudes  of  the  larger  earthquakes  were  determined  by  a 
method  based  on  duration  of  the  seismic  signal.  Energy  calculations  show  that  the  total 
seismic  energy  released  is  four  orders  of  magnitude  less  than  the  energyused  to  inject  water 
into  the  reservoir. 

Uinta  Basin  Observatory 

The  Uinta  Basin  Observatory,  UBO,  located  approximately  16  km  south  of  Vernal, 
Utah  (Figure  I)  and  65  km  WNW  of  the  Rangely  oil  field,  was  the  closest  seismograph 


*  Publication  authorized  by  the  Director,  U.S.  Geological  Survey. 
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Station  operating  continuously  during  the'  period  studied.  Nine  short-period  vertical 
seismometers  are  arranged  in  triangular  patterns  with  a  3-component  station  located  at 
the  center  of  the  array  (Geotech  Corp.,  1965).  Initially,  the  high-gain  seismographs  were 
operated  at  a  magnification  of  300,000  (R.  A.  Hartenberger,  December  1970,  written 
communication).  On  November  14,  1962,  the  gain  was  increased  to  400,000  and  on  July 
3,  1963,  to  600,000  where  it  remained  for  the  rest  of  the  study  period.  The  data  are 
recorded  on  16  mm  photographic  film  and  magnifications  reported  by  UBO  are  measured 
on  a  viewing  screen  that  magnifies  the  film  strip  10  times. 


Fio.  1.  Location  map  showing  the  Uinta  Basin  Seismological  Observatory  (UBO)  and  the  Rangely 

oil  field. 

Rangely  Oil  Field 

The  Rangely  oil  field  is  in  Rio  Blanco  County,  Colorado,  and  surrounds  the  town  of 
Rangely  (Figure  1).  Oil  was  discovered  in  1933,  but  full-scale  development  of  the  field 
was  not  started  until  1943  and  continued  into  1949.  When  the  field  was  completed,  481 
wells  had  been  put  down  into  the  Weber  sandstone  to  an  average  depth  of  6,700  ft.  The 
geological  structure  of  the  reservoir  is  an  asymmetrical  anticline  about  12  miles  long  and 
5  miles  wide  with  the  longitudinal  axis  trending  NW-SE.  A  normal  fault  with  a  NE-SW 
strike  has  been  recognized  cutting  the  longitudinal  axis  about  one-third  the  distance  from 
the  SE  end  of  the  reservoir  (Bleakley,  1964). 

Water-flooding  of  the  Rangely  field  was  begun  late  in  1957.  Water  under  pressure  is 
injected  into  converted  oil  wells  to  force  oil  that  cannot  be  recovered  by  conventional 
pumping  methods  toward  a  lower  pressure  well  where  it  can  be  pumped  to  the  surface. 
The  water-injection  wells  are  being  fed  by  a  centrally  located  pumping  station  that 
maintains  a  pressure  of  approximately  1,200  psi.  Wells  near  the  periphery  of  the  field 
were  the  first  to  be  converted  to  water  injection.  By  September  1965,  97  wells  had  been 
converted  and  by  September  1969,  water  was  being  pumped  into  202  wells  (Figure  2). 
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Fig.  2.  Oil  wells  (dots)  and  water  injection  wells  (squares)  in  the  R^ngcly  oil  field  in  September  1965  Oop) 

and  in  September  1969  {bottom). 

uses  Network 

In  November  1967  the  U.S.  Geological  Survey  installed  four  temporary  portable 
seismographs  (Figure  3)  around  the  Rangely  oil  field  in  an  effort  to  locate  the  active 
earthquake  areas.  These  instruments  recorded  continuously  for  an  8-day  period,  revealing 
a  pattern  of  seismic  activity  that  was  related  to  the  areas  of  the  field  where  fluid  pressure 
was  high.  In  September  1969,  a  16-station  seismograph  network  was  installed  around  the 
oil  field.  The  data  are  transmitted  via  telephone  line  from  the  oil  field  to  Menlo  Park, 
California,  where  they  are  recorded  on  16-min  photographic  film. 


AS-106 


; 

k 

[■ 


!» 


I ' 

f  . 


RBOSP  ^ROGFT^SS  REPORT  I 


1 


k 


d 


;  *  %xjiL  I  miU  I < 


h'r^irfr-  •r.r-iteir'^' 


-  H  JAMES  F.  GIBBS,  JOHN  H.  HEALY,  C.  BARRY  RALEIGH,  AND  JOHN  COAKLEY 


+ 


A 


I09«00' 

+  40’’05' 


A 


I08M7.5' 

+  40*I0' 


RANGELY 

[27 

+  + 


I  KM 

■  EARTHQUAKE  EPICENTER  '  ' 

A  PORTABLE  SEISMIC  STATION 

A 


Fig.  3.  Location  of  portable  recording  instruments  and  earthquake  epicenters  in  the  vicinity  of  Rangely 

oil  field  for  an  8-day  recording  period  in  November  1967. 


Seismic  Interpretation 


Earthquakes  occurring  in  the  vicinity  of  Rangely  were  identified  from  UBO  data  by 
the  characteristic  P-wave  arrival  pattern  and  the  S—P  time.  The  outer  triangle  of  stations 
was  used  to  determine  the  azimuth  from  UBO  to  the  epicenter  (110  +  4°).  The  distance 
was  determined  by  reading  the  time  interval  between  the  /*-wave  and  ^-wave  arrivals, 
using  the  3-component  station  at  the  center  of  the  array.  If  the  S—P  time  was  between 
5  and  10  sec  (Figure  4)  and  the  azimuth  was  correct,  the  earthquake  was  considered  to  be 
from  the  Rangely  area.  This  corresponds  to  an  epicentral  distance  of  about  40  to  85  km. 
Those  earthquakes  with  an  S—P  time  near  the  extremes  of  5  and  10  sec  are  probably 
outside  the  boundary  of  the  oil  field.  The  onset  of  the  .S'-wave  arrival  was  usually  easy  to 
pick  on  the  NS  horizontal  seismometer  and  together  with  the  P-wavc  arrival  sequence  on 
the  outer  triangle  of  stations  made  Rangely  area  earthquakes  easy  to  identify. 

The  plot  of  magnitude  versus  frequency  (Figure  5)  shows  that  most  earthquakes 
greater  than  magnitude  1  were  read  from  the  UBO  records.  A  complete  listing  of  the 
earthquake  data  is  contained  in  a  U.S.  Geological  Survey  report  by  Gibbs  er  al.  (1972). 
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Magnitude  Determinations 

Richter  magnitude  (A/J  was  estimated  from  Rangcly  area  earthquakes  using  ampli¬ 
tude  and  period  measurements  of  the  ^-wave  arrival,  usually  the  largest  arrival  on  the 
seismogram.  The  amplitudes  were  measured  peak-to-peak  in  millimeters  on  a 

viewing  screen  that  magnifies  the  film  strip  20  times.  The  estimated  Wood-Anderson 
amplitude,  was  calculated  from 


ir-; 


(I) 


The  factor  1  /4  is  necessary  to  reduce  the  measured  amplitude  to  center-to-peak  and  to 
account  for  amplitudes  being  measured  on  a  20-power  viewing  screen  (UBO  reports 
magnifications  for  a  10-power  viewing  screen).  R(H^)  is  the  ratio  of  the  magnification  of 
the  standard  Wood-Anderson  to  the  standard  Johnson-Matheson  as  a  function  of 
frequency  (Figure  6). 

The  Richter  magnitude  is  obtained  from 


A/jl  =  Tog  Log  Ao 


(2) 


where  —Log  Aq  =  2.8  at  65  km  (Richter,  1958,  page  342). 

Estimate  of  from  P-wave  amplitude  iAj„p)  was  also  calculated  using  the  largest 
P-wave  amplitude  in  the  first  second  of  record.  Estimates  of  A/^  were  obtained  from 


Mlp  =  Log 


-|^P(lT)j-l-2.8  +  0.71 


(3) 


Equation  (3)  includes  the  empirical  constant  0.71,  which  corrects  for  the  average 
amplitude  difference  between  P  and  S  waves.  The  graph  of  Mi^  versus  Mi^p  shows  a 


AS- 108 


Htnt'iiTnnjiiiir’ 


JAMF:S  F.  GIBBS,  JOHN  H.  HEALY,  C.  BARRY  RALEIGH,  AND  JOHN  COAKLEY 


Fig.  5.  Frequency  versus  magnitude  curves  for  earthquakes  in  vicinity  of  Rangcly  recorded  between 
November  1962  and  January  1970  at  UBO.  The  hwer  curve  shows  the  number  of  events  in  a  A-magnitude 
interval.  The  upper  curve  shows  all  events  greater  than  or  equal  to  a  given  ^-magnitude  interval. 


reasonable  agreement  between  the  two  sets  of  magnitudes  (Figure  7).  A  comparable 
magnitude  could  be  determined  for  those  earthquakes  when  the  S  wave  could  not  be 
measured. 

Magnitudes  Reported  by  UBO 

We  calculated  Unified  Magnitude  as  described  by  the  staff  at  the  Uinta  Basin  Observa¬ 
tory  as 


RBOSP  PROGRESS 


AS-109 


i*«W-iii 


«#*M‘  *»  ■liii,l><><^#>l 


y 


ite*MiC/«nr  j 


JAMES  F.  GIBBS,  JOHN  H.  HEALY,  C.  BARRY  RALEIGH  AND  JOHN  COAKLEY 


S-MflGNITUDE 

Fig.  7.  Relation  of  F-Magnitude  to  ^-Magnitude  for  earthquakes  in  vicinity  of  Rangely.  A  constant 
(0.71)  has  been  added  to  P  Magnitude  and  a  line  with  a  slope  of  1  is  drawn  through  the  points. 


B  in  this  special  case  was  extrapolated  from  tabulated  values  supplied  by  Carl  F.  Romney, 
AFTAC  (written  communications),  for  65  km  and  has  the  value  0.726.  Mi^  versus  M  is 
plotted  in  Figure  8. 


Estimated  from  Coda  Length 

Data  from  the  U.S.  Geological  Survey  network  (Rangely,  Colorado)  of  16  stations 
operating  since  mid-September  1969  were  used  in  conjunction  with  data  recorded 
simultaneously  at  UBO  to  determine  a  magnitude  scheme  based  on  signal  duration 
(Tsumura,  1967).  The  coda  length  is  defined  as  the  time  in  seconds  that  the  amplitude  of  a 
seismic  trace  remains  above  an  arbitrary  level  (1  cm  on  the  viewer  screen  was  chosen). 
The  seismographs  of  the  U.S.  Geological  Survey  network  at  Rangely  are  driven  off  scale 
by  local  earthquakes  of  about  magnitude  1.  A  coda-length  magnitude  curve.  Figure  9, 
was  constructed  by  making  all  coda-length  measurements  on  seismograms  recorded  by 
U.S.  Geological  Survey  instruments.  The  magnitudes  were  calculated  using  amplitude 
measurements  from  U.S.  Geological  Survey  data  for  the  small  earthquakes  and  from 
UBO.  data  for  the  larger  earthquakes.  A  magnitude  is  determined  from  either  Rangely 
or  UBO  and  compared  with  a  signal-duration  measurement  obtained  from  the  Rangely 
network  in  order  to  plot  the  curve  (Figure  9).  By  straight-line  extrapolation,  an  estimate 
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Fig.  8.  Relation  of  Unified  Magnitude  to  5  magnitude  (MU  for  earthquakes  in  vicinity  of  Rangely. 
A  line  with  a  slope  of  1  drawn  through  the  point  yields  a  {/-Magnitude  of  1.1  for  a  5-Magnitude  of  0. 

of  larger  magnitude  earthquakes  can  be  read  directly  from  the  graph,  simply  by  measur¬ 
ing  signal  duration  from  the  Rangely  seismograph.  Magnitudes  between  2.5  and  5.0  can 
be  estimated  with  an  accuracy  of  about  a  ±0.5  magnitude  unit  using  this  method. 

The  largest  earthquake  recorded  from  the  vicinity  of  Rangely,  of  Richter  local  magni¬ 
tude  3.4  (estimated  from  coda  length),  occurred  on  August  5,  1964.  This  earthquake 
caused  minor  damage  in  the  town  of  Rangely. 

Energy 

The  energy  radiated  as  elastic  waves  from  each  earthquake  was  calculated  using  the 
Gutenberg- Richter  (Richter,  1958,  p.  366)  formula. 


Log  E  =  9.9  + 1.9  0.024  Mj} 


(5) 


The  largest  shock,  that  of  August  5,  1964,  has  an  energy  of  12  x  10^*  ergs.  The  total 
seismic  energy  released  per  year  is  plotted  together  with  the  number  of  earthquakes 
greater  than  magnitude  1  (Figure  10).  Total  seismic  energy  released  was  lO^"^  ergs, 
equivalent  to  the  energy  of  one  earthquake  of  magnitude  3.9.  The  energy  necessary  to 
inject  each  barrel  of  water  can  be  approximated  by  multiplying  the  volume  of  water  by 
the  pressure  of  injection,  a  constant  pressure  of  approximately  1,200  psi.  A  straight- 
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Fig.  9.  Coda  length-magnitude  curve  determined  from  Rangely  area  earthquakes.  Open  circle  indicates 

Rangely  data;  closed  circle  indicates  UBO  data. 


forward  calculation  of  this  energy  is  1.4  x  10*^  ergs  per  barrel  of  water  passing  through 
the  wellhead.  From  November  1962  to  January  1970,  the  volume  of  water  injected  is 
4.8  X  10®  barrels  and  the  energy  of  injection  is  approximately  6.7  x  10^^  ergs,  10“^  times 
the  energy  released  as  seismic  waves.  Although  the  number  of  earthquakes  occurring  in 
the  oil  field  area  during  this  period  is  quite  large,  approximately  1,000,  the  energy 
released  as  seismic  waves  is  quite  small. 
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Discussion 

The  number  of  earthquakes  recorded  per  month  does  not  appear  to  correlate  with  the 
net  fluid  injection  per  month  (Figure  1 1);  but  if  the  number  of  earthquakes  is  plotted  on  a 
yearly  basis  and  compared  with  the  fluid  injected  per  year,  there  is  an  apparent  correlation 
(Figure  12).  The  absolute  value  of  fluid  injection  docs  not  seem  to  affect  the  seismic 
activity;  rather,  changes  in  the  quantity  of  fluid  injected  are  related  to  changes  in  the 
number  of  earthquakes  recorded:  if  fluid  injection  increases,  there  is  a  corresponding 
increase  in  seismic  activity.  This  relation  holds  true  for  all  years  except  1969,  when  a 
modest  increase  in  fluid  injection  is  accompanied  by  a  dramatic  decrease  in  earthquake 
activity. 
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Fig.  10.  Total  number  of  earthquakes  per  year  (greater  than  magnitude  1.0)  and  total  energy  released  per 

year  for  earthquakes  occurring  in  the  Rangely  area. 

It  is  probable  that  a  correlation  could  be  developed  if  earthquake  hypocenters  could 
be  accurately  located  within  the  oil  field  and  checked  against  fluid  injected  into  wells  near 
the  active  areas.  In  parts  of  the  field  where  there  are  no  natural  faults,  injections  do  not 
seem  to  produce  earthquakes,  even  when  the  fluid  pressures  are  quite  high,  and  the 
inclusion  of  data  from  the  stable  parts  of  the  oil  field  tends  to  degrade  the  correlation. 
There  are  other  factors  that  tend  to  mask  the  relation  between  fluid  injection  and  earth¬ 
quakes.  For  example,  changes  in  the  pattern  of  injection  wells  can  redistribute  the  water 
available  for  injection  such  that  zones  that  are  near  their  critical  pressure  will  receive  less 
water  as  new  wells  are  initiated  in  zones  that  are  far  below  their  critical  pressure.  Even 
though  there  may  be  a  substantial  increase  in  fluid  injection  as  these  new  wells  are  intro¬ 
duced,  less  water  is  injected  into  the  zones  that  could  produce  earthquakes. 

A  comparison  of  the  number  of  injection  wells  as  of  September  1 965  with  the  number  of 
injection  wells  as  of  September  1969  (Figure  2)  shows  that  development  of  new  wells  in 
the  low-pressure  parts  of  the  field  may  account  for  the  apparent  decrease  of  seismic 
activity  during  the  latter  part  of  this  period.  The  development  of  a  complete  pressure 
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Fig.  11.  Total  and  net  fluid  injected  per  month  is  shown  along  with  the  number  of  earthquakes  per 
month,  Rangely  area.  The  net  fluid  injected  is  the  total  fluid  injected  minus  the  oil  and  water  produced. 


history  of  this  oil  field  would  be  a  major  undertaking,  and,  although  we  believe  it  would 
improve  apparent  correlations  between  fluid  pressures  and  earthquakes,  we  doubt  that 
this  improvement  would  warrant  the  effort  required.  The  data  we  now  have  strongly 
suggests  that  a  correlation  exists.  New  data  being  developed  in  a  restricted  part  of  the 
field  will  adequately  test  the  hypothesis  relating  fluid  pressure  to  seismic  activity. 
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MINIMUM  KECOIIDLXG  TIMES  FOR  DETERMINING  SHORT-TERM 
SEISMICITY  FROM  MICROEARTIIQUAKE  ACTIVITY 

By  Allan  R.  Sanford  and  Suuendila  Singh 

ABSTRACT 

The  minimum  microeorthquake  recording  time  needed  to  determine  short¬ 
term  seismicity  (10  years)  has  been  established  for  regions  having  the  same  de¬ 
gree  of  activity  as  Socorro  (400  ‘shocks  with  M  ^  0  each  year).  An  analysis 
of  more  than  five  years  of  microeorthquake  data  indicates  that  about  six  months 
of  recording  are  required  for  a  reliable  estimate  of  the  earthquake-frequency 
relation  from  which  seismicity  is  calculated.  If  the  slope  b  of  the  earthquake- 
frequency  low  is  assumed,  the  short-term  seismicity  can  be  determined  from 
four  weeks  of  recording.  The  danger  in  the  latter  method  is  that  a  small  error  in 
b  leads  to  relatively  large  errors  in  seismicity. 

long-term  seismicity  (100  years)  predicted  from  the  present  level  of  micro- 
earthquake  activity  in  Socorro  is  much  lower  than  the  known  seismicity  of  the 
region.  The  number  of  mlcroearthquakes  would  have  to  increase  by  a  factor  of 
20,  or  the  value  of  b  decrease  approximately  20  per  cent,  before  prediction 
would  agree  with  the  long-term  seismic  activity  of  the  region. 

Introduction 

Generally  determination  of  seismicity  begins  with  establishing  the  earthquake- 
frequency  relation  (Richter,  195S,  p.  359) 

logio  iV  =  A  —  bMf  (1) 

where  .1/  is  magnitude,  N  is  the  number  of  shocks  of  magnitude  M  or  greater, 
and  A  and  h  are  constants.  The  two  constants  define  the  best  straight-line  fit  be¬ 
tween  logio  N  and  .1/,  h  equaling  the  slope  and  A  being  the  logarithm  of  the  number 

of  shocks  of  magnitude  zero  or  greater.  ^  •  j  t 

For  most  areas,  equation  (1)  must  be  determined  from  a  limited  peiiod  of  re¬ 
cording.  To  determine  seismicity,  the  short-term  relation  is  extrapolated  to  longer 
time  periods,  e.g.,  10  or  100  years,  and  a  calculation  of  the  strongest  shock  to  occur 

within  that  period  is  made. 

A  recent  modification  of  this  procedure  is  to  assimie  a  value  of  I  (usually  about 
1.0)  and  to  detemiine  A  by  simply  counting  the  number  of  shocks  above  a  mim- 
mum  magnitude  within  a  specified  period  of  time  (Isacks  and  Oliver  1964,  p.  1965). 
llie  relation  obt.ained  in  this  manner  can  also  be  extrapolated  to  longer  time  periods 

f 

for  determining  seismicity. 

An  obvious  factor  that  may  lead  to  errors  in  seismicity,  whether  one  assumes  or 
measures  b,  is  time  variations  in  the  level  of  activity.  If  microcarthquakes  are  used 
to  establish  the  earthquake-frequency  relation,  they  must  be  recorded  for  a  length 
of  time  sufficient  to  average  out  exjiectcd  variations  in  the  number  and  strength 
of  the  shocks.  Minimum  recording  times  for  reliable  estimates  of  seismicity  can  be 
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established  in  regions  where  microcarthqiiakcs  have  been  monitorcKl  for  a  number  of 
years.  At  Socorro,  recording  of  nearby  microearthquakes  lias  been  continuous  for 
more  than  five  years.  Time  vanations  in  these  data  liavc  been  analyzed,  and  some 
reasonable  estimates  can  be  made  of  the  minimum  recording  times  required  to 
establish  the  short-term  (10  years  or  less)  seismicity  of  a  region  having  the  same 
degree  of  activity  as  Socorro. 
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Minimum  amplitude  of  ground  motion  —  18 m/i 
Fig.  1.  Number  of  microcarthquakes  versus  S-P  interval  for  one-year  periods. 


Seismicity  Basijd  on  Microe-arthquake  Count 

The  simplest  method  of  establishing  seismicity  is  to  assume  a  value  of  6  and 
determine  A  from  a  count  of  inicroearthquakes.  ^Minimum  recording  time  for  this 
method  was  arbitrarily  taken  as  that  period  for  which  there  is  a  95  per  cent  prob¬ 
ability  of  being  vithin  a  factor  of  3  of  the  expected  number  of  shocks.  A  threefold 
variation  in  the  number  of  microearthquakes  leads  to  a  difference  of  0.5  in  the  mag¬ 
nitude  of  the  strongest  shock  predicted  for  a  specified  interval  of  time  (assuming 
b  =  1.0). 
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MINIMUM  MK'KOEAKTHQUAKE  HECORUING  TIME  .  4 

Minimum  recording  time,  as  defined  above,  was  found  b}'  subdividing  the  total 
five  years  of  Socorro  microearthquake  activity  into  progressively  smaller  intervals  of 
time.  The  first  interval  considered  was  one  year.  Figure  1  shows  the  number  of 
shocks  as  a  function  of  S-P  for  the  years  19G2  through  19GC.*  The  minimum  level 
shock  considered  for  these  graphs  had  a  maximum  vertical  ground  motion  of  18.2 
m/i.  The  magnitude  corresponding  to  tliis  ground  amplitude  is  0.5,  a  figure  based 
on  an  average  epicentral  distance  of  10  km,  an  average  angle  of  incidence  of  20 
degrees,  and  a  factor  of  2  correction  for  converting  vertical  to  horizontal  ground 
motion.  A  relatively  large  minimum-level  shock  was  used  for  the  one-j  ear  interv’al 
to  include  data  taken  when  the  magnification  of  the  seismograph  was  quite  low. 


In  Figure  1,  the  average  number  of  shocks  a  yetir  is  183,  the  standard  deviation 
22.  As  the  variation  from  year  to  year  is  quite  small,  a  single  year  of  recording 
would  have  provided  an  accurate  estimate  of  the  seismicity  for  the  entire  five-year 
period. 

The  second  time  mterval  considered  was  three  months,  and  the  data  used  were 
from  the  years  1962,  1963,  and  19G4.  The  minimum-level  shock  in  the  analysis  was 
magnitude  0.2  (maximum  vertical  ground  motion  9.1  m/i).  For  the  three-rnonth 
interv'al,  the  average  number  of  shocks  an  interval  was  96,  the  standard  deviation 
28.  Thus,  if  Socorro  activity  had  been  monitored  for  only  three  months,  there  would 
have  been  about  a  95  per  cent  probability  of  being  within  a  factor  of  2  of  the  ex¬ 
pected  number  of  shocks  (96).  A  twofold  variation  in  microearthquake  activity  is 

•  The  marked  absence  of  events  with  S-P  less  than  0.0  second  and  greater  than  2.0  seconds 
plus  the  two  well-defined  peaks  in  the  curves  of  Figure  1  indicates  that  shocks  are  7wt  uuiformly 
distribated  in  space  about  the  recording  station. 
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equivalent  to  about  a  0.3  variation  in  -4  and  tlicrefore  i.s  fairly  insignificant  in  the 
determination  of  sei.smicity. 

The  third  time  interval  was  four  weck.s.  I'ifty-two  intervals  from  data  gathered 
in  19G2,  10G3,  19G4,  and  19GG  were  used,  and  the  minimum-level  .-hock  was  .set  at 
magnitude  0  (maximum  vertical  ground  motion  G.l  m^).  ligure  2  Is  a  graph  of  the 
cumulative  number  of  four- week  intervals  versus  the  number  of  events  an  interval. 
The  dashed  line  is  the  observed  distribution,  the  solid  line  is  the  the<jretical  dis¬ 
tribution  (Poisson),  as.suming  the  earthcpiakes  are  random,  mutually  independent 
events.  Tlie.se  cur\'cs  indicate  that  the  Poisson  distribution  may  be  of  little  value 
in  estimating  the  minimum  recording  time  required  to  establish  seismicity. 

The  ob.served  distribution  of  events  per  interv^al  in  I'igure  2  indicates  a  95  per 
cent  probability  that  the  number  of  events  for  any  arbitrary  selected  inten’al  would 


Number  of  Eorthquokes  an  Interval  (One  Week) 

Fio.  3.  Observed  and  theoretical  curv'es  showing  the  distribution  of  one-week  intervals 

according  to  the  number  of -events  an  interval. 


be  within  a  factor  of  3  of  the  expected  number  (37).  Thus,  four  ueeks  would  be 
the  minimum  recording  time  using  the  simple  count  method. 

An  indication  of  the  danger  involved  in  using  shorter  interx'als  of  recording  was 
obtained  by  subdividing  the  activity  in  the  years  1962  through  1964  into  one-week 
intervals.  In  this  anal5'sis,  the  minimum-level  shock  was  again  set  at  magnitude  0. 
Figure  3  is  a  graph  of  the  cumulative  number  of  one- week  interv'als  versus  the  num¬ 
ber  of  microearthquakes  an  interval.  As  before,  the  difference  between  the  ob¬ 
served  and  the  theoretical  Poisson  distribution  is  significant.  The  observ'ed  distribu¬ 
tion  indicates  a  95  per  cent  probability  that  the  number  of  shocks  in  an  arbitrarily 
selected  one-week  inter\'al  would  have  been  within  a  factor  of  10  of  the  expected 
number  (10).  A  tenfold  variation  in  the  number  of  shocks  is  equivalent  to  a  unit 
change  in  the  value  A  when  b  is  assumed  to  be  1.0.  Thus,  one  week  of  recording 
could  lead  to  estimates  of  seismicity  over  longer  periods  that  are  one  magnitude 
unit  too  high  or  too  low. 


Selsmicity  Based  on  the  Complete  E.arthquake-I'Heqi  enc\  Law 

More  reliable  estimates  of  seismicity  can  be  established  if  b  as  well  as  .4  is  deter¬ 
mined.  Tliis  atmroach  naturally  requires  more  recording  time 
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number  of  events  is  needed  to  accurately  define  b,  the  slope  of  the  logio  N  versus  M 
curv'e. 

Listed  in  Tabic  1  are  the  values  of  A  and  b  for  10  dilTcivnt  three-month  intervals 
of  micioearth(iuakc  activity  in  the  Socorro  area.  The  magnitude  interval  over 
which  the  values  of  -4  and  b  were  determined  (b}'  the  least  square  method)  was 
0.03  to  0.03.  The  last  two  columns  give  the  .‘Strongest  shocks  in  10-  and  lOO-ycar 
periods  as  determined  by  extrapolating  the  relation  obtained  for  three  montlrs. 

The  strongest  shocks  for  the  10-year  jieriod  range  in  magnitude  from  2.6  to  4.4. 
The  lowest  magnitude  was  obtained  from  a  three-month  interval  during  whieh  the 
total  microearthquake  count  was  low  (75)  and  the  value  of  6  high  (1.333) ;  the  high¬ 
est  was  obtained  for  a  period  during  wliich  the  count  was  near  average  (13S)  but 
the  value  of  b  low  (0.857). 

TABLE  1 

Valves  of  A,  6,  and  Strongest  Shocks  Determined  from  Three-Month 
Periods  of  Microeakthqvake  .Activity 


Three-Month  Inter^-als 

Total 

U  g  0.03 

A 

Stroni^est 
Shock,  10  yrs. . 

Strongest 

Shock. 

ICO  yrs. 

July  19G2-Sept.  19G2 . 

123 

2.118 

0.885 

4.2 

5.3 

Oct.  19G2-Dec.  1962 . 

185 

2.297 

0.942 

4.1 

5.2 

Jan.  1%3-March  1963 . 

129 

2.142 

0.938 

4.0 

5.1 

June  1963-Aug.  1963 . 

115 

2.088 

0.968 

3.8 

4.8 

Oct.  1963-Dec.  1963 . 

78 

1.917 

0.945 

3.7 

4.8 

Nov.  1963-Jan.  1964 . 

86 

l.%9 

1.258 

2.8 

3.6 

Jan.  1964-Mar.  1964 . 

138 

2.166 

0.857 

4.4 

5.6 

May  1964- July  1964 . 

139 

2.173 

1.004 

3.8 

4.8 

Sept.  1964-Nov.  1964 . 

141 

2.177 

0.943 

4.0 

5.1 

July  1965-Sept.  1965 . 

75 

1.918 

1.333 

2.6 

3.4 

Totals . 

1209 

10.073 

37.4 

.\verage . 

121 

1.007 

.3.74 

Standard  Deviation . 

34 

0.158 

0.59 

A  much  smaller  range  in  magnitudes  for  the  strongest  shocks  would  be  obtained 
if  h  were  assumed  equal  to  1.0,  the  average  value  for  the  10  three-month  interv'als. 
Although  this  procedure  ma}’  make  estimates  of  seismicity  more  consistent,  it  may 
not  necessarily  make  them  more  accurate.  If  the  assumed  value  of  b  is  15  per  cent 
in  error,  the  magnitude  of  the  strongest  shock  vill  be  in  error  about  0.5  (a.ssuming  a 
level  of  activity  about  the  same  as  Socorro). 

The  data  (listed  in  Table  1)  indicate  only  a  6S  per  cent  probability  of  obtaining  b 
Mithin  about  15  per  cent  of  the  expected  value  (1.0)  from  a  single  three-month 
period  of  recording.  An  alternate  interpretation  of  the  data  in  Table  1  is  that  a  re¬ 
liable  estimate  of  b  requires  a  certain  minimum  number  of  obsen'ations.  If  the 
three  interv'als  having  less  than  100  observations  are  removed  from  the  list,  (1) 
the  average  number  of  events  per  inteiwal  becomes  139  with  a  standard  deviation 
of  23,  and  (2)  the  average  value  of  b  becomes  0.934  with  a  standard  deviation  of 
0.049.  This  result  suggests  a  reliable  estimate  of  b  can  be  obtained  if  the  number  of 
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observations  exceeds  about  loO,  regardless  of  the  time  intcrv'al  over  which  the  ob¬ 
servations  are  made. 

The  computations  of  seismicity  described  above  indicate  that  small  variations  in 
h  can  be  as  important  as  large  variations  in  the  number  of  shocks  counted.  Thus, 
unless  observational  data  everywhere  indicate  that  h  is  some  type  of  universal  con¬ 
stant,  the  sim})le  count  method  must  be  used  with  caution  even  when  allowance  is 
made- for  time  variations  in  activity. 

J.,ong-Tekm  Selsmicitv 

It  is  probable  that  the  short-term  seismicity  of  the  Socorro  area  could  have  been 
accurately  c.stimatcd  from  a  single  six-month  j)eriod  of  recording.  The  .stronge.st 
shock  predicted  for  a  10 -year  period  would  be  about  3.7,  which  appear.^  reasonable 
on  the  ba.sis  of  .strong  shocks  that  have  occurred  during  the  past  41-  rears  (11 
shocks  3/  ^  2.6). 

The  place  where  the  microcarthquake  method  fails  is  in  predicting  long-term 
seismicity.  Socorro  has  had  relatively  severe  earthquakes  within  the  past  100 
years — shocks  that  all  e\'idence  indicates  must  have  originated  within  the  recording 
range  of  our  pre.sent  instrumentation  (Sanford,  1963,  p.  149).  The  stronge.st  .<hock 
within  the  last  100  years  occurred.  November  15,  1906,  was  felt  over  100,000  square 
miles,  and  had  a  maximum  intensity  of  VIII.  Magnitude  6  appears  to  be  a  reason¬ 
able  estimate  for  the  strength  of  this  shock  (Gutenberg  and  Richter,  1956). 

An  extrapolation  of  equation  (1)  using  the  average  values  given  in  Table  1  pre¬ 
dicts  the  strongest  shock  for  a  100-ycar  pei-iod  at  magnitude  4.7.  The  dilference 
between  prediction  and  observation  is  significant.  The  number  of  microcarth quakes 
recorded  in  the  Socorro  area  would  have  to  increase  by  a  factor  of  20,  or  the  value  of 
b  decrease  to  0.8,  before  prediction  would  agree  with  ob.serx’ation. 

The  difficulty  in  determining  long-term  sei.smicity  from  microeartlujuake  data 
is  not  surprising.  Short-term  variations  in  this  activity  make  it  hard  to  estimate 
seismicity  over  periods  of  time  as  short  as  10  years.  It  is  jirobable  that  long-term 
variations  in  activity  exi.st  and  can  only  be  accounted  for  by  basing  estimates  of 
long-term  {M?i.smicity  on  data  obtained  over  an  extended  period  of  time. 
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The  well-fracturing  operation  is  modeled  by  a  band  of  uniform  pressure  and 
TY-%  I  \a<U  of  unifoim  shear  stress  acting  in  a  cylindrical  cavity  in  an  infinite  body.  Two 
r-.v.-vj^ing  regions  of  induced  stress  are;  either  end  of  the  pressurized  interval  where  the  tan- 
f  --*..il  .‘trLss  is  zero  (the  vertical  stress  is  approximately  95  per  cent  of  ihe  pressure)  and  the 
n:  :v:  Cl  the  packed-off  interval  where  the  tangcuual  stress  equals  the  pressure  (the  veitJcal 
It  roro).  The  tectonic  stresses  are  the  overburden  load  and  two  un.known  principal 
<.  r  j.  n*.«l  stresses  that  cause  easily  determined  stress  concentrations  at  the  well  bore.  All 

•  A,  ’itircd  stresses  are  moJiUed  to  account  for  the  interstitial  po.^e-fiuid  pressure.  It  is  found 

'.Krce  situations  are  of  interest;  (1)  the  induced  vertical  stress  is  less  than  the  overburden 
.r'Tcurc;  (2)  the  induced  vertical  stress  and  the  instantaneous  shut-in  pressure  are  greater 
i.*  tiie  overburden  pressure;  (3)  the  induced  vertical  stress  is  greater  than  the  overburden 
Lut^the  instantaneous  shut-in  pressure  is  less  than  the  overburden  pressure.  la  (1) 
•r  f.Mrture  is  vertical  and  the  stresses  are  dcte.’-minable.  In  (2)  the  fracture  is  horizontal  and 
'  ••.•csh;s  are  indetenninate.  In  (3)  the  fracture  is  initisdly  horizontal  but  becomes  vertical 
*.  jTop.ag.ates  av/ay  from  the  well,  the  vertical  and  minimum  horizontal  compressions  are 
'v.arjablc,  and  the  other  principal  stress  is  bounded  by  a  set  of  inequalities.  Several  exam- 

•  \:e  presented  in  which  the  tectonic  stress  states  appear  to  be  rela.\cd — approximately 
• .  ...ih-2t  to  the  ‘hydrostatic’  pressure. 


Introduction 


% 


’  •-•-.■‘.r.itudes  and  orientations  of  the  prin- 
"  *'  -‘T?  in  tlic  earth's  crust  hold  the  key 
-'-Ivuental  understanding  of  the  defor- 
•'  I  t.'.c  c.arth’s  crust.  The  problem  con- 
•  -  t.'i’s  paper  is  the  determination  of 
of  the  components  of  the  stress 
*'•  i  paint  in  the  earth's"  crust.  This  is 
•‘•.'■.g  data  nomially  obtained  during  ?. 
--.i-fracturiiig  operation.  The  prob- 
•  ‘  r.rjpjitjjjjQjj  of  the  stresses  can  be 
vcvl.  It  will  be  seen  later  that  the 


•  :'>a  of  the  fracture  v/ill  be  perpen- 
lox't  compressive  principrd  .stress 
.  oric-ntation  can  bo  determined 

'  oi  an  oriented  inflatable  iiackcr, 

•’ •■vd  by  Fraser  and  Pettxtt  [19G2]. 
'"d  U  dlis  [1957]  wore  the  first  to 
c.-ndusively  the  influence  of  the 


■  ? tresses  on  the  orientation  of  a 
•••  induced  fracture,  but  they  did  net 
-  calculate  tlie  magnitudes 

from  tlio  well-fracturing  drta. 
-  ...  interpreted  as  the  natu- 

earth's  stresses  in  (he  \acinity  of 


the  hole  before  they  b?came  modified  by  the 
drilling  of  the  hole.)  Scheidegger  [1962],  using 
Hubbert  and  Willis'  theory  with  some  modifi¬ 
cation  calculated  the  tectonic  stresses  for  sev¬ 
eral  examples.  Scheidegger  [19C9]  expressed 
concern  over  the  fact  that  the  theory  of  Ilub- 
bert  and  Willis  did  not  account  for  the  creation 
of  vertical  stresses  ajid  hence  was  unable  to 
predict  the  fonnation  of  horizontrd  fractures 
e.xccpt  in  very  special  cases  under  a  set  of  re¬ 
strictive  assumptions.  In  hopes  of  alleviating 
this  difficulty,  Scheidegger  [1960]  suggested 
that  the  pressurized  interval  could  be  consid¬ 
ered  a  spherical  pressure  center.  Had  far-fold 
effects  been  of  interest,  tin's  assumption  could 
bo  used  by  virtue  of  St.  '\’cnant's  principle  [cf. 
Timoshenko  and  Goodier,  1951,  p.  33];  but  the 
rtresses  cxi.stnig  at  iiic  perimeter  of  tlio  borehole 
are  of  primary  interest,  for  it  is  here  that  the 
fracture  forms.  The  geometrical  differences  be¬ 
tween  a  uniformly  loaded  sp.Iierical  caadtv  r..nd 
a  cylindrical  hole  loaded  over  a  finite  vertical 
distance  are  so  great  that  slrc.-s  distributions 
cannot  be  expected  to  be  similar  near  the 
boundaries;  in  fact,  they  are  quite  different. 
The  tectonic  stresses  calculated  by  our  c.vact 
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analysis  agree  more  closely  ^yith  those  calcu¬ 
lated  using  Hubbert  and  Willis’  model  than 
with  those  obtained  by  Scheidegger. 


RALPH  0.  KEULE 

5ni,  S*.  =  maximum,  minimum  horiz/.:  -;. 
field  tectonic  stress.  . 


The  Stress  Field  Induced  by  the  Fracturing 
OrER-KTlON 


List  oj  symbols. 

X  =  radial  coordinate. 

$  *=  angular  coordinate. 

g  t=  vertical  coordinate. 

^  =  dummy  integration  variable. 

R  «=  unknown  function  of  r. 

5  *=  unknown  function  of  X. 

a  =  radius  of  borehole, 

cf'  =  Xa. 

p  sf  Xr. 

7o,  h  =  modified  Bessel  functions  of  the  first 
kind,  order  0,  1,  respectively. 

Ko,Ki—  modified  Bessel  functions  of  the  second 
kind,  order  0,  1,  respectively, 
c  *=  height  of  pressurized  interval. 

U  -  height  (length)  of  packer. 

V  =  Poisson’s  ratio. 

$  =  stress  function. 

V*  —  Laplace  operator. 

k  «  effective  boundarj*  porosity. 

S  =  value  of  shear  induced  by  packer. 
p  =  pressure  in  pressurized  zone. 

At  A',  B,  B'  -  integration  constants, 

p,  =  breakdown  pressure. 

pf  tt=  flowing  pressure, 

p,  c=  shut-in  pressure. 

P,  =  formation  fluid  pressure. 

Lo  =  total  overburden  pressure. 

Sie$,  Si„,  S,rr  =  tangential,  vertical,  and  radial 
components  of  the  induced  stress. 

=  effective  tangential,  vertical,  and 
radial  components  of  the  induced  stress. 
Srit,  S„rr  =  tangcTitial,  vertical,  and  radial 
components  of  the  natural  (tectonic) 
stress. 

<r«M,  =  effective  tangential,  vertical,  and 

radial  components  of  the  natural  stress, 
etc.  =  sum  of  the  tangential  components 
(<S,- » »  +  e ») . 

etc.  "  sum  of  the  clTcctivc  tangential 
components  (cr<  5#  +  <r„»j)* 

=•  natural  tectonic  stress  in  x  direction 
before  hole  was  drilled  (far-ficld  stress). 
=  natural  far-field  tectonic  stress  in  y 
direction. 


S$i, 

is.. 


Symbols  defined  in  text. 

Si(x)  =  sine  integral  0—>' X. 

Ci(x)  =  cosine  integral  x  — >  «> . 
h  =  spacing  used  in  numerical  hr.. 

=  integration  parameters  arising  ir. ;  , 
method  of  integration, 
functions  of  'F  that  aiLsc  i.a  } 
method  of  integration. 
a',  h'  -  integration  limits  in  Filon’s  r-.-.:- 

The  model  In  an  open-hole  wcll-fr . 
job,  the  prospective  producing  horizon  ; .. . 
off  and  fluid  is  introduced  into  this  z'-:.' 
fluid  pressure  is  increased  to  a  maxinvzr.  . 
P„  at  which  time  the  formation  fraev:: 
Figure  1).  (All  pressures  (Pej  Pf»  F,,  c.- 
to  be  taken  as  the  actual  bottom-hole  i 
If  pressures  at  the  surface  are  reconl  -i 
necessary'  to  correct  for  hydraulic  h'. 
viscous  losses.)  Fluid  continues  to  be  ,  - 
into  the  zone,  and  the  pressure  stabih:  - 
value  Pf,  the  flowing  pressure,  which  c-::.- 
is  intermediate  to  the  formation  fluid  ] 

(P.)  and  P..  (Occasionally  P.  =  Pf ' 
the  fracturing  operation  has  been  cor.' 
and  the  well  shut-in,  the  pressure  qui*;i- 
bilizes  to  a  constant  value.  This  is  ca. 
instantaneous  shut-in  pressure,  Pj.  If 
turing  fluid  viscosity  is  low,  it  is  fou... 

Pf  =  P,  when  appropriate  corrections  a.”-  ■ 
for  viscous  losses.  A  schematic  diagra.-^ 
well-fracturing  setup  is  shown  in  Fipne  - 
The  fracturing  fluid  induces  a  unih-nr.  "■ 
sure  over  the  packed-off  interval.  (1 1- 
variation  over  the  pressurized  interv.d  o  • 
gravitational  effects  is  smaller  than  the  r 
for  accuracy  of  the  method.)  This^'^  ‘ ' 
modeled  exactly  by  a  uniform  band  of 
in  a  cylindrical  hole  (Figure  1).  y'J 
also  tends  to  force  the  packers  away  f- yy 
pressurized  iiiterval,  but  any  such  ^ 
of  the  packers  is  prohibited  by  friction-- 
that  arise  at  the  contact  between  the  J  ■ 
and  the  wall  of  the  borehole.  These 
forces  are  modeled  by  a  uniform  ban* 
stress  over  the  interval  of  the  hole  ^  . 
the  packers;  the  magnitude. of  the 
anccs  the  force  applied  to  the  packer  • 
fracturing  fluid  pressure  (Figure  2). 


TECI 


V  1.  Schematic  of  pressure  ver? 

i;inionu  shear  stress  distt. 

.  ••  ,,.a,l  uiiifonn  pressure  cousin 

:  d.e  effect  of  the  fracturing  oi 
:•  hole.  These  distributions  are 
of  boundary  conaitions  iOi 

. of  the  stress  field  induced  i 
fracturing  operation.  Beca 
;.:oblcm,  the  total  solution 
■;  V  superposing  the  solutions 

•  : :  'iM-ms.  That  is,  the  effect  of 

•  •  :tc  and  the  shear  bands  ma; 

.  r-;:;»r.ately  and  then  combined 

-  •  fleet. 

'  ;.rtssurc.  The  problem  of  tl 
;  :-.-.*ure  acting  on  the  walls  o 

•  . -•.drical  hole  was  solved  by 
i"..-*  v;i!ucs  of  the  stresses  in  c} 

•  •.'•'1  .arc  given  by 

-  --  f"  lX’<irK,(Xa)K.(>.r) 

Tf  J, 

•a,  >.-.)K,CXr)  -  XrKo(Xr)iC.(X< 
-  2;t  -  r)]K.(Xa)K.(Xr)! 

••  t-rcfis-Xa 
'  a  • 

■  - aoX 

*  f  {Xq /vo(Xa)iv,(Xr) 

-5 

sm  Xa  sm  a 


\ » 


i)(Xa) 
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...  too  small  to  show  on  ih,. . 

3  r-^P.tecI  in  Figure  3  arc  for  a  r. 

T. 

cussion.  Inspection  of  Figure  3 
lain  the  complete  stress  state  ir ; 
oundarj'  of  the  hole  by  the  wdl-f:,. 
fion,  because  it  is  known  that  .5 
and  Si„  =  —P  over  tlie  pres;-;.*. 

.  It  is  of  particular  interest,  to  r. 
Juced  Si„  is  virtually  nonc.\i.-;t.  , 
ler  end  of  the  packed-ofT  interval  „• 
same  region  the  magnitude  of  e... 
o.  Thus  it  is  only  at  either  c;;*: 
i-off  interval  that  the  fracturi: 
iduces  a  stress  with  which  a  h  r 
re  could  be  fonned.  The  secor.d  ; 
t  is  that  over  most  of  the  press’;:.  • 
=  P.  This  is  exactly  the  v.i! 
for  Stt$  if  the  cylinder  is  loaded  x! 
length.  It  is  not  surprising,  tl’.er..  • 
e  found  that*  the  present  anai>v\ 
closely  with  Hubbert.  and  ; 

icnt  in  which  a  uniform  pressu:.- 
over  the  entire  hole. 

’ECTO^•Ic  Stress  at  the  BoREne: : 
concentration  of  the  regional 
xiste^^of  regional  stres.ses  in  t'.-  ■ 
hichl^oil  well  is  drilled  will  cau-- 
trations  in  the  vicinity  of  the  ! 
j  analysis  it  is  assumed  that  the  •.  • 
on  of  the  tectonic  stresses  i.s 
ue  thickness  of  the  packed-ofT  z< 
>tion  is  similar  to  the  prcviou'l; 
assumption  that  the  fluid  pre 
ty  significantly  over  the  packed-  ‘ 
Any  regional  stress  distributio; 
of  a  borehole  can  therefore  1).'  ’ 
by  the  superposition  of  two  stre.-- 
0  of  the  earth,  two  compression.'.  T 
lishcd  by  choosing  the  horizon  t;t!  - 
of  the  reference  coordinate  sy.-*- 
Old  to  the  direction  of  the  r"-- 
al  regional  stress.  We  arbitrar;!;.-  ‘ 
and  Sk!  =  5„„  where  and 
cipal  reglonal  stresses.  The  so!;;- 
ilcm  is  then  equivalent  to  an  it ' 
ith  a  hole  at  the  center,  tliat  ;• 

'  the  stresses  iS„,  and  in  tlie 
5,  respectively.  This  solution  oar.  ’ 

V  superpo.dng  the  solution  for  t;.  • : 

1  to  S*,  alone  and  the  solution  i  ' 


TECl'ONIC  STRESSES 

^  .»j  of  tlicsc  arc  just  the  problem  of  the 
"  a  hole  subjected  to  simple  tension, 

*  *’j  by  Kirsh  [ISSS]  [cf.  Timoshenko 

*  .  .It  f,  *1952,  pp.  7S-S0].  Superposing  the 

‘.ions,  and  setting  the  radial  variable 

*  >  t.hc  l^rebole  radius,  we  see  that 


problem.  They  are  the  central  region  which 
comprises  SO  per  cent  of  the  packed-off  interval 
where  the  total  stresses  are  given  by 

Srt  —  S,l  —  Si,  —  0 


-  +  2(5,.,  ~  cos  2d 

-  0 

-  0  • 


(31) 


r.vximum  values  for  5„s9  occur  at  cither 

■  -  '.T  or  :r(2n  -f  l)/2,  n  -  1,  2— depending 
.  ,-.,'ti;or  S«,  or  S„,  is  the  larger.  Tlie  mini- 

...  Tjiues  for  S„»e  occur  at  90°  from  the 
.:n  values.  From  this  point  in  the  paper, 
V  _  .a-sumed  that  5„  =  5.;  is  the  greatest 
■?.il  stress  (least  comprc.ssive)  and  that 
•  i.'  the  least  horizontal  stress  (greatest 

'■"ive).  (It  is  imperative  to  remember 
-j;ivc  stresses  are  tensions  and  that  com- 

•  - .  are  negative  stresses.  Furthermore, 
-  .•  indicated  by  P»  are  positive  quantities, 

•J  equation  a««  =  — P  states  that  (t.i  is 
'7-:«*.-^ive  (negative)  stress  equal  in  mag- 
-  t'»  the  pressure  P.)  Thus  the  least  com- 
■ (maximum)  is  given  by 

=  35„-  =  35„j-  5,.  (.32) 

•  *.-0  points  along  the  circumference  of  the 

•  -•'■.••‘re  5„,9  attains  this  maximum  value  are 
•*  p.-3int3  marked  by  the  intersection  of  the 

•  -'Uindary  and  the  y  coordinate  direction 
arcction).  These  points  are  important  bc- 
•*  die  well-fracturing  operation  superposes 

5,9,  over  the  existing  compressive  S„r>9 
tlie  fracturing  fluid  pressure  is  increased, 

■  d'.ccs  where  the  hoop  stress  first  become.? 

‘  .arc  c.xactly  the  points. in  question. 

completes  the  necessary  solutions.  The 
'■  state  at  any  point  of  tlie  borehole  wall  is 
'*  by  superposing  the  naturally  occurring 
(equation  32)  and  the  induced  stresses 
fmm  Figure  3. 

'‘•""Inre  criterion.  As  will  bo  seen  below, 
'  ''rcss  distributions  at  the  hole  boundary 
the  wcll-fracturing  operation)  are  such 

■  .’  is  not  really  necessary  to  use  a  complex 
.re  criterion  for  the  prediction  of  fracture 

“■'•'ion.  There  arc  two  di.'tinct  region.?  of  the 
‘•  i-ofi  interval  that  arc  of  interest  in  the 


Srr  -  -P 

5., =  -Lo  *  (33) 

5110.. ,  ^  35n,  —  5„y  +  P  —  3S„j  —  5„i  +  P 

and  the  verx*  ends  of  the  interval  where  the 
stresses  are 

Srr  =  ~P  Sr.  =  5,,  =  5,,  =  0 

s„  =  -lo  +  0.94P  (34) 

59,  =  35„  -  5„,  =  35,2  -  5„, 

where  the  value  of  5i„  (0.94P ,)  is  read  from 
Figure  3.  In  both  cases,  there  arc  two  large  com¬ 
pressive  stresses  and  one  stress  approaching 
tensile  values  as  the  hydraulic  fracturing  pres¬ 
sure  is  increased;  i.e.,  at  tlie  ends  of  the  zone 

5.,  is  approaching  tensile  values,  and  in  the 
center  Soo  is  becoming  a  tension.  Moreover,  in 
most  wells  the  fracturing  pressure  is  within  25 
per  cent  of  the  overburden  load,  and,  as  a  con¬ 
dition  of  incipient  fracture  is  approaclicd,  the 
interesting  stress  states  arc  composed  of  two 
nearly  equal  'confining'  pressures  and  a  third 
stress  that  is  either  zero  or  slightly  tensile. 
Such  stress  states  are  analogous  to  the  'triaxial 
extensions.’  tests  of  llandin  and  Fauburn  [1955] 
and  Heard  {Griggs  and  llandin,  1930,  p.  202]^ 
The  data  from  these’  tests  lead  me  to  believe 
that  a  value  of  20  kg/cm*  in  tension  would  1 
probably  be  adequate  to  fracture  the  rocks  dur¬ 
ing  a  well-fracturing  operation,  regardlf'ss  of  • 
the  depth  of  the  formation  being  treated,  and  ‘j 
this  value  will  be  used  as  the  fracture  criterion.  • 
Some  di.seretion  must  be  used  in  all  cases  before  ' 
this  value  is  insisted  upon.  However,  this  sim-  * 
pie  result  is  iictually  in  accord  with  both  the 
modified  Mohr-Coulomb  theory  of  failure  and  | 
Kchle's  [1961,  pp.  C9-70]  generalized  fracture  j 
criterion. 

The  fractures.  In  accordance  with  the  frac¬ 
ture  criterion  decided  up)on  in  the  preceding 
section,  all  that  is  necessary  for  the  initial  ion  of 
a  fracture  is  that  either  5,,  in  tlie  ccntml  por¬ 
tion  of  the  pressurized  interval  or  S,,  at  the 
ends  of  the  pressurized  interval  reach  a  value 
of  20  kg/cm’  in  tension.  In  the  latter  case  the 
fracture  that  is  fonned  will  be  horizontal;  in 
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.e  former  ease  the  fracture  \nll  be  vertical.  It 
fj  important  to  note  that,  if  it  forms,  the  verti¬ 
cal  fracture  will  begin  propagating  along  the 
plane  perpendicular  to  the  least  compressive 
horizontal  stress  The  problem  that  re¬ 

mains,  when,  is  to  decide  in  any  particular  case 
whether  the  horizontal  or  the  vertical  fracture 
■  forms. 

The  C.acuLATioN  of  the  Tectonic  Stresses 

• 

The  effective  stress.  .411  the  solutions  for  the 
stresses  used  in  the  paper  were  obtained  from 
the  classical  theoiy'  of  elasticity  which  assumes 
ft  simple  continuous  material.  In  reality,  how¬ 
ever,  wc  arc  dealing  with  porous  mediums.  Even 
though  a  continuum  theory  for  porous  mediums 
was  not  used,  a  first  approximation  to  such  a 
theory  is  obtained  when  the  results  arc  modified 
by  taking  into  consideration  the  concept  of 
‘effective’  stress.  The  effective  stress  is  defned 
as  ' 

total  "b  hPo5i;  (35) 

where  P,  is  the  interstitial  pore  fluid  pressure 
and  is  the  Kronecker  delta.  The  constant  k 
has  been  called  the  effective  boundary  porosity 
--  by  Tersagi  [1936].  Many  attempts  have  been 
^^nade  to  determine  values  for  the  effective 
boundaiy'  porosity,  and  arguments  as  to  its  sig¬ 
nificance  have  been  presented  in  the  literature 
for  more  than  twenty  years.  At  present,  I  feel 
safe  in  saying  that  the  value  of  the  effective 
boundary  porosity  varies  between  1.0  for  a 
sphere  pack  at  low  pressures  to  about  O.So  for 
some  sandstones  when  loaded  in  a  prescribed 
- '  manner  at  high  pressures  [Fatt,  195S].  G.  Mor¬ 
ris  (personal  conversations)  finds  that  the  value 
is  between  0.93  and  0.97  for  most  of  the  rocks 
tested  in  this  laboratory'  and  that  it,  in  effect, 
depends  on  the  compressibility  of  the  grains 
and  matrix  as  opposed  to  the  compressibility  of 
the  framework  a.s  a  whole.  Because  the  accuracy 
of  this  method  for  calculating  stresses  is  not  e.x- 
pected  to  be  better  than  ±10  per  cent,  and  also 
because  the  well-fracture  data  are  rarely  of 
better  accuracy,  for  simplicity  T  will  use  the 
value  of  1.0  for  the  effective  boundary  porosity. 
Tho  stresses  of  interest  to  us  (equations  33  and 
34)  become 


KEHLE 

=  3(S„  -b  P.) 


\ 


(Tfi 


C„  -  -(P  -  P.) 
<r„  “  ~{L.  —  P.) 


-  (5.,  +  P.)  +  (P  -  P.) 


;  ot -tho  largo 

/,th  the  edge  ot  »  f?* 

Vase  with  the  wcll-fracturuig 

;'.ho  cracUs  ^ubic^cd  to  an  ^ 


-  -r  i:.;  r  V-  ‘  ^  suDjccicu 

in  the  central  region  of  the  packed-off  int-.n  very  little  work  is  neces  . 

and  .  Ire.  In  the  tveU-fractunns  operation. 

IT,.  =  -(P  -  P.) 

<r„  =  -(L.  -  P.)  +  0.94(F  -  P.) 

<P|,  =  3(S„,  +  Po)  -*  (Sny  +  P*) 


at  either  end  of  the  “interval. 


(Xcous)'  Scf  »d\°'scparath«  ^ 

«  ia  a. 

■  respect  to  the  exic  stresses*,  thu 

Initiation  of  the  fracture.  Upon  rea. ..j  hold  the  crack  ope' 

the  pressure  P.  at  which  the  formation  more  than  the  far-field  stre. 

tures,  either  a„  at  either  end  of  the  inter.  ,-^t  of  the  frac 

<r„  in  the  center  of  the  interval  has  attain  .  ?  ,  ^osi  p.  490].  This  leads  t 

value  of  approximately  30  kg/cm’  in  tensi  j  lC  •  conclusions.  First,  if  the  well 

the  effective  induced  vertical  stress  cr,,,  =  ;  i  ‘  upon  completion  of^  t 

(P.  -  P.)}  is  less  than  the  effective  overbu:  )  ii)  the  pressure  wiU  quick 

load  {L.  -  P.},  (T..  could  not  have  attain-  . /.^niig  ^  ’  giigi^tly  greater  than  t 

value  of  20  kg/cm’  in  tension  and  therefor-.  j^cting  perpendicularly  to  t 

must  have  been  (To$  that  reached  the  enr  <  ^  j{  this  stabilization  press 

value.  The  fracture  that  forms  is  vertical.  T.  ‘  than  the  perpendicular 

the  relationsliip  calculation  [of. 

'‘?i  losi  P.  339]  .ovcals  U.at  streas  cone 

' '  \  nf  such  great  magnitudes  exi- 
;^:  ^^l^that  it  would  certain^  P* 
,..01  luc  _  the  volume  oi 

...  'Ti-iic  would  uicrea^e  mu  _  • 

■  ;J  nd  because  no  additional  flu.d 
S  to  the  system,  tho 
■,.ase.  The  final  stable 


35 


5.y  =  35.,  “  5., 

=  20  -  P.  -  P.  (• 

is  obtained  from  (36). 

If,  on  the  other  hand,  the  effective  fractcr. 
pressure  is  greater  than  the  o\erburdca  1 

either  a.,  or  a.,  could  have  reached  the  cn*.  uuai  .  maintain 

value,  because  our  analysis  to  tliis  point  i? :  pressure  sufficient  . 

definitive  enough  to  permit  us  to  decide  :•  ..jj-c  but  insufficient  to  ^ 

tween  the  two  jxissibilitics.  If  o-es  had  reac.  would  be  , 

the  critical  value,  a  vertical  fracture  would  !.  perpendicular  far-fieid  s  ^ 
been  formed  and  (3S)  is  again  obtained.  On '  is  that  the  pressure  ~ 

other  hand,  if  a.,  reached  the  critical  valu-  •  stress.  In  this  case,  e 
horizontal  fracture  would  have  formed  and  r.  ‘  collapse,  thereby  .  .^.-n 

would  have  been  less  than  20  kg/cra’  in  tenr.  •  fracturing  fluid.  * 

Thus  (36)  yields  1  .e  until  the  collapse  is  . 

■  -V  because  of  the  large  volume  of  bqu 

'•  ■  'd  h“^o  disperse  into  the  fonuatr 

;  .;  the  pressure  in  tbo  fluid  -s  b.ough 

■  .l-je  slightly  greater  than  the  fw-fi 

-  .la,  n-e  arrive  at  the  ■ 

■  -neons  shnt-in  pressure  is  approa 
'  -.1  to  the  regional  stress  acini  „  porf 

■  to  the  piano  of  tho  fracture. 
r:.e  second  important  point  coiiccin- 
■;,n  of  froctiiro  P'“l"Sn‘io“  ^ 
tlv  stressed  medium  when  propa. 


(36) 


3S„  -  S.,  <  20  -  P.  -  P. 

This  equation,  besides  serv'ing  as  a  bound  to  n 
horizontal  stresses,  has  an  important  inter:  ■' 
tation.  The  difference  between  the  princ’.:- 
horizontal  stresses  as  well  as  their  alisf  ■ 
magnitudes  relative  to  the  overburden  load  ■  ^ 
tcrmincs  whether  or  not  a  horizontal  friot.*' 
forms.  The  propensity  for  the  formation 
horizontal  fractures  increa.sc3  as  the  princ,- 
horizontal  stre.sscs  approach  equality. 

The  instantaneous  shit-in  pressure  P «• 
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;(S..  +  P.) 

-  (&.  +  >.)  +  (P  -  P.) 

ral  of  the  packed-off  ii.- 


’-P.) 

..  -  P.)  +  0.94(P  -  P.) 

^  +  P.)  -  (5,,  +  P,) 

.d  of  the  interval. 

I  of  the  fracture.  Upon  r- 
re  P,  at  which  the  formatio-  ; 
r  a„  at  either  end  of  the  into.-, 
renter  of  the  inter^'a!  h.as  atta;- 
iproximatcly  30  kg/cm’  in  ten.-, 
e  induced  vertical  stress  a*.,  = 

^  b  less  than  the  effective  ovc  r’r. . 
-  P«},  cr„  could  not  have  ati.;;- 
)  kg/cru*  in  tension  and  ihcM-: 
been  that  reached  llie  rr 
fracture  that  forms  b  vertical. ' 
oship 


3S,a  S, 


ml 


*=  20  -  P.  -  P. 

I  from  (36). 

.e  other  hand,  the  effective  fr;:'* 
'  greater  than  the  overburd..-, 
or  ^^^uId  have  reached  the 
auso^^  analysis  to  this  pol:'.'  . 
enough  to  permit  us  to  dcc'.  i 
two  possibilities.  If  veo  had  : 

I  value,  a  vertical  fracture  won! ! 
cd  and  (3S)  is  again  obtained. 
d,  if  a„  reached  the  crlticrd  v . 
fracture  would  have  fonr.cJ 
e  bcc-n  less  than  20  kg/cm’  in 
)  jiclds 

~  <  20  -  P.  -  Pc 

tion,  besides  scr\'ing  as  a  bound  ‘ 
stresses,  has  an  important  i:.’ 
no  difference  between  the  r*  ' 
stresses  as  well  as  their  r.-- 
•"3  relative  to  the  overburden  ! 
whether  or  not  a  hori?:onta'.  ” 
ho  propen-sity  for  the  fori’.v* 
fractures  increases  as  tlic  ; : 
stresses  approach  equality. 
tanla7icous  shut-in  pressure  /*. 


of  the  largo  stress  concentrations  associ- 
with  the  edge  of  a  crack  (especially,  as  is 
"/case  with  the  well-fracturing  operation, 

:  .....  the  crack  is  subjected  to  an  internal  pres- 

'  * very  little  w'ork  is  necessary'  to  extend  a 
[  In  the  well-fracturing  operation  most 

>  *.  energy  is  expended  in  overcoming  fric- 
■  -cd  (vbcous)  losses  and  in  separating  the 
,..i  of  the  fracture.  Furthemrore,  the  crack  is 

•  ..•.-■ximately  an  extremely  elongate  ellipse. 
's,f  f.nicture  wall  displacements  are  veiy  small 

•  «  .t  respect  to  the  extent  of  the  body  and  the 

•'.--d  regional  gradient  of  the  stresses;  thus 
pressure  necessary  to  hold  the  crack  open 
dightly  more  than  the  far-ficld  stress 
I  •;  pcr];>cndicularly  to  the  plane  of  the  frac- 
-  .’cf.  Sneddon,  1951,  p.  400].  This  leads,  to 
t  ;:r.portant  conclusions.  First,  if  the  well  is 
■  ui  immediately  upon  completion  of  the 
v'.rirg  treatment,  the  pressure  will  quickly 
. .  -..:e  to  a  value  just  slightly  greater  than  the 

•  i-.dd  stress  acting  perpendicularly  to  the 

re  plane.  If  this  stabilization  pressure 

•  :  'J!  bo  greater  than  the  perpendicular  re- 

:d  stress,  a  simple  calculation  [cf.  Muslchel- 
. ,  1053,  p.  339]  reveals  that  stress  concen- 

•.■..r.'*  of  such  great  magnitudes  e.xist  at  the 
-»  of  the  crack  that  it  would  certainly  prop- 
;  -  This  would  increase  the  voliune  of  the 
■  and,  because  no  additional  fluid  is  be- 
;  • .;  plied  to  the  sj'stcm,  the  pressure  would 
■•••••.  The  final  stable  pres.surc  would  be,  in 

•  tbit  pressure  sufficient  to  maintain  the 

•  .rv  but  insufficient  to  cause  further  prop- 

•  1.  and  it  would  be  approximately  equal 
*  perpendicular  far-ficld  stress.  The  other 

•  bty  is  that  the  pressure  is  less  than  the 
■li  stress.  In  this  case,  the  fracture  would 

-  collapse,  thereby  decreasing  the  volume 

•  !e  to  the  fracturing  fluid.  This  will  cori- 
^  xr.til  the  coUapsc  is  total  (which  is  not 

1  •■■cause  of  the  large  volume  of  liquid  that 
:  have  to  disperse  into  the  formation)  or 
pressure  in  the  fluid  is  brought  up  to 
“  .'lightly  greater  than  the-far-field  stress, 
“■e  .arrive  at  the  conclusion  that  the  in- 
■‘••'■us  shuMn  pressure  is  approximately 

•  *''•  the  regional  stress  acting  pcrpendicu- 
'  ■  ’.*.•?  plane  of  the  fracture. 

■•''•.■r.d  important  point  concerns  the  di- 
'■•f  fracture  propagation  in  a  nonuni- 
f'rc-ssed  medium  when  i)ropagation  is 


,*  . 

caused  by  a  uuifonn  pressure  acting  on  the 
walls  of  the  crack.  Because  of  the  large  stress 
concentrations  at  the  edge  of  the  crack,  little 
work  is  needed  to  cause  the  material  at  the  end 
of  the  crack  to  break,  regardless  of  the  orienta¬ 
tion  of  the  fracture.  The  work  necessary  to 
separate  the  walls  of  the  crack,  however,  is  di¬ 
rectly  proportional  to  the  magnitude  of  the 
compressive  stress  acting  perpendicularly  to  the 
plane  of  fracture.  Thus  the  easiest  fracture  to 
propagate  is  one  that  is  perpendicular  to  the 
least  compressive  regional  stress.  Invoking  the 
principle  of  least  work,  we  can  say  that  a  frac¬ 
ture  that  propagates  because  of  the  application 
of  an  internal  pressure  "will  propagate  in  that 
direction  that  minimizes  the  work  necessary  to 
extend  the  fracture,  and  that  direction  will  be 
along  the  plane  that  is  perpendicular  to  the 
least  compressive  far-field  stress.  This  is  exactly 
the  principle  used  by  Anderson  [1951]  when  he 
discus.sed  the  problem  of  dike  orientation,  and 
Ode  [1957],  using  this  principle,  successfully  ex¬ 
plained  the  formation  of  the  dikes  in  the  Span¬ 
ish  Peak.s,  Colorado.  The  most  important  corol¬ 
lary  of  this  principle  is  that  the  shut-in  pressure 
is  approximately  equal  to  the  least  compressive 
regional  stress.  In  the  present  notation,  we  have 


-P.  =  -  5, 


h3 


(40) 


A  second  important  conclusion  is  th.at  the  fluid 
pressure  near  the  end  of  the  crack  during  the 
treatment  is  also  appro.ximately  equal  to  the 
least  compressive  regional  stress.  Unfortunately, 
the  pressure  near  the  end  of  the  crack  may 
bo  considerably  less  than  the  pressure  in  the 
packed-off  interval  of  the  borehole  because  of 
large  viscous  resistance  to  flow  along  the  frac¬ 
ture.  However,  techniques  are  available  by 
which  one  could  calculate  those  losses  and  thus 
deduce  the  pressure  at  the  end  of  the  crack 
[cf.  Howard  and  Fast,  1957;  Christianovitch 
et  al.,  1959]. 

The  stress  calculation.  Under  Initiation  of 
the  fracture  it  was  discovered  that,  if  the  mag¬ 
nitude  of  the  effective  induced  ci,,  was  less  than 
the  effective  overburden  load,  the  fracture  that 
formed  could  not  have  been  horizontal.  In  this 
case,  the  fracture  is  vertical;  it  begins  to  propa¬ 
gate  in  a  radial  direction  along  the  plane  perpen¬ 
dicular  to  the  lea-st  compressive  regional  stress 
and,  in  agreement  witli  the  principle  of  least 
work,  it  will  continue  to  propagate  in  this  dircc- 

|~^OSP  PROGRESS  RlhPOn^  j 


AS-128 


1 


I 


'■  1# 


f 


j 


s 


^  RALPH  0.  KEHLE 

Substitutin';  (40)  into  (3S),  we  have 

=  P,  -f-  _  2P,  -  20  (41) 


TECTOXI 


Thus  nil  the  principal  regional  stresses  have 
been  detennined  because  =  ~P,  and  = 

-1.. 

The  second  possibility  is  that  the  induced 
elective  =:  {0.04  (P,  -  P.)}  is  greater 
than  the  overburden  load  (plus  tlie  20-kg/cin* 
tensile  cohesive  strength  if  no  bedding  plane  is 
available  for  separation).  If  the  instantaneous 
shut-in  pressure  is  also  greater  than  the  over¬ 
burden  load,  there  is  little  reason  to  believe 
that  a  vertical  fracture  would  form  because  the 
induced  a.,  is  sufTicicnt  to  initiate  a  horizontal 
fracture,  and  the  shut-in  pressure  is  sufficient  to 
support  such  a  fracture.  Furthermore,  if  we  as¬ 
sume  that  a  vertical  fracture  did  form,  (41) 
indicates^  that  the  maximum  horizontal  com¬ 
pression  is  approximately  twice  the  overburden 
load.-At  best,  such  large  horizontal  stresses  are 
to  be  c.xpcctcd  only  in  areas  of  cxircmely  ac- 
ti\e  tcctonTsm.  Thus  it  is  assumed  in  sucii  cases 
that  a  horizontal  fracture  is  formed,  and,  un¬ 
fortunately,  no  statements  may  be  made  about 
the  values  of  the  horizontal  stresses  save  what 
can  be  inferred  from  (39). 

Another  possibility  is  to  find  that  the  induced 
<ri,,  is  sufficient  to  create  a  horizontal  fracture 
but  that  the  shut-in  pressure  is  less  than  the 
overburden  load,  wliich  indicates  that  the  final 
fracture  orientation  is  vertical  rather  than  hori¬ 
zontal.  If  it  is  assumed  that  a  vertical  fracture 
was  formed  initially,  we  again  use  (41)  and  ar¬ 
rive  at  a  large  horizontal  stress.  Smaller  values 
for  the  stresses  are  arrived  at  if  it  is  assumed 
that  a  horizontal  fracture  is  initiated  at  the 
^\ell  bore  and  that,  when  it  propagates  away 
from  the  hole,  its  orientation  is  more  and  more 
influenced  by  the  regional  stres.ses  rather  than 
by  the  local  stress  concentrations  caused  by  the 
e.xistence  of  the  hole  and  the  well-fracturing 
procedure.  Actually,  at  a  distance  of  onlv  three 
well-boro  radii  away,  these  concentration  effects 
are  negligible,  and,  because  of  the  con-sidcra- 
rions  of  th?  previous  section,  the  fracture  would 
be  expected  to  find  the  plane  across  which  the 
compression  was  smallest.  As  a  simple  example, 
consider  a  treatment  at  a  deptli  of  1525  meters,' 
so  that  the  overburden  load  is  350  kg/cm*.  As¬ 
sume  that  the  horizontal  stresses  are~  both  2S0 
kg/cm*  and  that  the  fluid  pressure  is  155  kg/cm*. 


Thus  the  natural  (max)  is  255  ke/ni  , ,  n  n  r  t  ^ 

the  fracturing  pressure  would  have  P,  =  P.  =  L..  In  tins 

proximatelv  435  ke/cm*  to  initiate  a  Penetrated  a  bedding  plane;  th. 

fracture.  On  the  other  hand,  only  370  ke  V  horizontal  strcssc. 

necessary  to  initiate  a  horizontal  fractu/*^™’”^^^’ 
tainly  the  horizontal  fracture  will  fona  Applic.vtion  to  Field  D.\ta 

few  feet  from  the  hole,  however,  the  strnt;  m  j  •  » , 

is  considcrablv  different  from  that  at  ^  desirable  tj'po  of 

boundary.  Tlic  fracture,  should  it  rem.nff 

zontal,  would  do  the  maximum  amount'.’  measurements  is  a 

in  propagating  rather  than  the  mininn- e-^'Pensive  procedure  and 
is  not  at  all  reasonable,  and  it  is  ren-'^'eri/^  is  rarely  employed  in  normal  field 

the  fracture  pronnirates  away  from  enough 

it  would  follow  the  lines  of  lea.^t  reri=tani'^^°"  available  when  measurements  are 
becoming  a  vertical  fracture  within  ffic  "r^ly 

feet  from  the  hole.  The  fracture  then  ^  of  the  fracturing  fluid 

into  the  hole,  thereby  pennittin-  the  automatically  recorded, 

have  direct  entry  into  the  vertical  regarding  the  vis- 

The  assumption  that  the  horizontal  f  fracturing  fluid,  the 

forms  and  then  becomes  vertical  additives,  the  cou- 

•  ion  of  propping  agent,  the  instantaneous 
pressure,  and  the  ‘pumping  rate’  are 


cr.. 


>  <tc9.  From  equation  39  and  tliC' 


equation  37,  the  definition  5„.  >  S„ 


cause  S^,  =  it  is  po.csible  to  7cdi!^°  of  this  type  of 

■presented  in  Figure  I.  The  pressure  at 


~P.  >  =  5.,  >  —  3P,  -f  L,  formation  fractures  is  usually  indi- 

pn  the  log.  Even  if  it  is  not,  it  is  gen- 
-I-  0.94P ,  -f  O.OGP.  quite  easy  to  determine  the  appropriate 
which  serves  to  bound  the  value  of  5,^.  ^  graph.  To  obtain  bottom-hole  pres- 
As  a  final  comment,  it  will  be  noticed  ’jfoi^ation,  it  is  only  necessary  to  correct 
lower  bound  imposed  by  (42)  in  the  pressures  for  the  hydrostatic  head 

horizontal-vertical  fracture  assumption'  fracturing  pressure,  for 

lead  to  a  contradiction.  If  '  (viscous)  losses  in  the  tubing.  These 

fions  are  .easily  obtained  by  employing 
P.  <  LJA  +  0.2P,  +  0.01  oP.  ‘rapids  that  are  presented  in  various  hand- 
this  would  say  that  5„,  >  S^,  which  coiu^®^’  Frac.  Guide,  19.59].  The  fric- 

the  definition  Sr^  >  In"'this  case  corrections  wliich  depend 

vious  that  such  a  fracture  system  did  nd"'  °°  ''^Jocity  of  the  fluid  in  the  tub- 
and  that  instead  a  vertical  fracture  was  Pimping  rate  neces- 

This  completes  the  analysis  for  the  tkro^  obtain  the  fracturing  pressure  is  the 
classes  of  well-fracturing  iata.  In  formation  fractured, 

tectonically  quiet  areas,  if  the  induced  than  the  rate  em- 

fufficient  to  cause  a  horizontal  fracture  f  the  remainder  of  the  treatment 

is  sufficient  to  support  the  overburden  t:i  ^  fracture  is  being  propagated  through 
ture  is  horizontal  and  the  horizontal  .strcs^^.^'°"‘  the  pur7iping  rate 

indeterminate.  If  the  induced  a<.,  is  c-i  t/ie  rate 

to  cau.se  a  horizontal  fracture,  but  breakdown  rather  than  the  in¬ 

sufficient  to  support  the  overburden  ilf^'  almost  all  of  the  information  on 
ture  is  horizontal  and  then  vertical  number  of  well  treatments  is  of  no 

and  5^,  is  bounded  by  (42).  If  tlie  indill'"  /  , 

is  less  than  the  efTcefive  overburden,  boffora-hole  pressure  measurc- 

ture  is  vertical  regardless  of  the  value'of/ during  wcll-fracturing  treatments 


P.,  and  S,.'is  given  by  ('(1)  Tiu> li'er.nturc  by  Godbej 

}  Km.  lui  Horner 
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TABLE  1.  Principal  Strcssca  Calculated  from  Well-Fracturing  Data 


T1 


Well 

Oklahoma  D* 

Oklahom.a  F* 

Venczuclaf 

West  Texasf 

North  Oklahoma  5 

Depth  meters 

1005 

8S5 

565 

485 

2135 

P„  kg/cm* 

70 

50  ■ 

70 

50 

190 

P„  kg/cm^ 

225 

170 

115 

85 

360 

kg/cm* 

290 

220 

125 

7511 
or  120  § 

405 

S„  kg/cm* 

226 

-195 

-125 

-105 

-490 

<rn.  (5.  -f  P,), 
kg/cm’ 

-150 

-1*45 

-  55 

—  55 

-300 

“  —0.94 
(F.  -  P.), 
kg/cm* 

-205 

-ICO 

-  50 

-  25 
or  —  65 

-200 

iSi,  kg/cm* 

Indeterminate 

-170 

-115 

-  85 

-360 

jSi,  kg/cm* 

Indeterminate 

<-170 

>-320 

-170 

-150 
or  —105 

-505 

«Si,  kg/cm* 

-365 

-290 

-170 

-150 

-505 

Si,  •*  kg/cm’ 

-410 

-365 

-ISO 

-175 

-145 

•  -625 

I- 


•  Godbey  and  Hodges  (195S].  Fractured  through  perforations,  thus  the  model  does  not  necessarily  ho!: 
in  these  cases,  . 

t  Van  Dam  and  Horner  [1957], 
j  Fraser  and  Pcitiit  [1961], 

§  North  Dover  Field,  Oklahoma.  Unpublished. 

11  Breakdown  pressure  listed  as  75  kg/cra*  but  does  not  appear  to  have  been  corrected  for  stau' 
head  of  -lo  kg/emP. 

^  Calculated  using  Scheidegger’s  [19G2]  modification  of  Hubberi  and  Willis’  [1957]  model  of  the  ^r.- 
acturing  operation. 

••  Calculated  using  Scheidegger’s  [19C0]  model  of  the  well-fracturing  operation. 


tho  strain  is  pennanent  and  thus  is  not  asso¬ 
ciated  a  buildup  of  the  clastic  energy' 

necessary  for  an  efiective  earthquake.  In  the 
present  method  the  stresses  are  mca.sured  di- 
.  rectly,  and  it  is  independent  of  any  previous 
permanent  defonnation;  thus  it  allows  immedi¬ 
ate  calculations  of  the  stored  elastic  energy’.  As 
such,  it  might  prove  to  be  a  valuable  tool  in 
the  prediction  of  earthquakes. 

Finally,  it  should  be  emphasized  that  only 
when  we  are  acquainted  with  a  detailed  history 
of  the  spatial  and  temporal  variations  of  the 
crustal  stress  fields  will  we  be  in  a  jK)sition  to 
understand  the  nature  and  causes  of  orogenesis. 
The  present  method  makes  it  possible  to  obtain 
part  of  this  necessary  infonnation. 
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EARTHQUAKE  DATA  FILE 


1974:  108,926  events  v;orl dwido  in  total  file. 

Abbreviations  used  - 

Autbori cy,  (follows  origin  time): 

A  Parameters  of  explosion  supplied  by  U.S.  Atomic  Energy  Coirymission 
(AEG). 

B  Parameters  of  hypocenter  supplied  by  University  of  California  at 
Berkeley. 

E  Some  or  all  parameters  of  explosion  (controlled  or  accidental) 
supplied  by  any  group  or  individual  other  than  AEG. 

G  Parameters  of  hypocenter  supplied  by  the  U.S.  Geological  Survey 
(U.S.G.S.)  for  any  area  other  than  Island  of  Hawaii. 

•H  Parameters  of  hypocenter  supplied  by  the  U.S.G.S.  Hawaiian  Volcano 
Observatory. 

J  Parameters  of  hypocenter  supplied  by  St.  Louis  University. 

L  Parameters  of  hypocenter  supplied  by  Lamonc-Doherty  Geological 

Observatory. 

M  Hypocenter  based  on  macroseismic  information. 

P  Parameters  of  hypocenter  supplied  by  California  Institute  of 
Technology. 

R  Param.eters  of  hypocenter  supplied  by  University  of  Nevada. 

S  An  NEIS  solution  based  on  use  of  dense  local  networks,  a  local 
crustal  model,  or  other  methods  not  routinely  applied  by  NEiS. 

V  Parameters  of  hypocenter  supplied  by  Virginia  Polytechnic  Institute 
and  State  University. 

VJ  Parameters  of  hypocenter  supplied  by  University  of  Washington. 

Z  Moninstrumental  time  and  location. 

Second-order  hypocenter  determination  by  CG5/N0S/ERL/GS  using 

incomplete  or  less  reliable  data. 

Depth  control:  D,  restrained  to  agree  with  reported  . 

cr,  r ^Tr¥i"n^  by  geophysicist;  H,  restrained  at  nornal  deptn  (33  km] 

when  data  are  not  sensitive  to  depth  for  a  shallow  locus. 
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Magnitudes :  BRK,  Berkeley,  California;  BUG,  Bucarest,  Rumania; 
C'CS/N'OS,  Coast  and  Geodetic  Suv'vey/Natiorial  Ocean  Survey;  ERL, 
Environmental  Research  Laboratories;  GOL,  Golden,  Colorado;  GS, 
Geological  Survey;  MOS,  Moscov/,  USSR;  MOX,  Moxa,  Germany;  OXF, 

Oxford,  Mississippi;  PAL,  Palisades,  Mew  York;  PAS,  Pasadena, 
California;  QUE,  Quetta,  Pakistan;  SLM,  St.  Louis,  Missouri;  TEH, 
Tehran,  Iran;  TUL,  Tulsa,  Oklahoma;  UPP,  Uppsala,  Svaeden. 

Int  (intensity)  map:  USE,  United  States  Eoj’tluiv.akos ;  EQd,  Earthqvxike 
Notes;  PDE,  Preliminary  Determination  of  Epicenters^ 


Phenomena: 

D  •.  Diastrophism 


F,  surface  faulting;  U,  upl i ft/subs idence; 

D,  faulting  and  upVi ft/subsidence. 

-  T,  tsunami  generated;  Q,  possible  tsunami. 

-  S,  seiche;  Q,  possible  seiche. 

-  V,  earthquake  associated  with  volcanism. 

-  C,  coal  bump;  E,  explosion;  I,  collapse;  R,  rockburst. 


T  Tsunami 
S  Seiche 
V  Volcanism 
N  Montectonic 

0  Waves  generated  -  A,  acoustic  wave;  G,  gravity  vave;  B,  both 

A  &  G;  T,  T-wave. 


Region  number  (RN):  Geographic  region  number  as  described  in  "Seismic 
and  Geographical  Regionalization,"  E.A.  FI  inn,  E.R.  Engdahl ,  and  A.R. 
Hill,  Bulletin  of  the  Seismological  Society  of  America^  Vol .  64, 

Number  3,  Part  II. 

Cultural  Effects  (CE):  C,  casualties;  D,  damage;  F,  felt. 

Qua! ity/Number  of  Stations  (Q/$) :  Qual i ty  i ndi caters  are  usual ly 
on  an  A,  B,  C,  or  D  basis  (A-very  accurate;  B~gobd;  C~fair;  D-poor) 
when  used  in  combination  for  deep  shocks  (e.g.  BBA)  these  represent 
the  accuracy  of  (in  order)  epicenter,  origin  time,  and  depth. 

Marsden/Degree  Square  (MAR  DG):  Numbering  system  dividing  the  world 
into  10^  squares  (MAR)  and  P  subsquares  (DG).  Complete  Earthquake 
Data  File  available  sorted  by  this  method. 
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All  abbreviations  used  to  indicate  an  earthquake  source  in  the  data 
file  are  listed  below.  (Sources  for  each  event  are  coded  on  the 
tape  in  columns  6-8.) 

AEC  U.S.  Atomic  Energy  Commission 

ALG  Algiers,  Algeria 

ALI  Alicante,  Spain 

ALM  Almeria,  Spain 

APA  Apatity,  RSFSR,  USSR 

API  Apia,  Samoa  Is. 

BCI  Bureau  Central  International  Sesmologie,  Strasbourg,  France 
BLA  Blacksburg,  VA 

BNS  Bensberg,  German  Federated  Republic 
BOG  Bogota,  Colombia 
BRA.  Bratislava,  Czechoslovakia 
BRK  Berkeley  (Haviland),  CA 

BSS  Bulletin  of  the  Seismological  Society  of  America 
BUC  Bucharest,  Romania 

CAN  Canberra,  Australia  Capital  Territory 
CAR  Caracas,  Venezuela 

CFR  Charles  F .  Richter  (see  Richter,  1958,  in  References) 

CGS  Coast  and  Geodetic  Survey 

CHC  Chapel  Hill,  NC 

DJA  Djakarta,  Java,  Indonesia 

EQH  Earthquake  History  of  the  United  States  (See  References) 

G-R  Gutenberg-Richter  (see  Gutenberg  and  Richter,  1954,  in  References) 

GS  U.S.  Geological  Survey,  Denver,  CO 

HEL  Helsinki,  Finland 

HVO  Hawaiian  Volcano  Obsy. ,  Kilauea,  HI 

ISS  International  Seismological  Summary,  Kew,  England 

JMA  Japan  Meteorological  Agency,  Tokyo,  Japan 

JOH  Johannesburg,  South  Africa 
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bb 
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LIS  Lisbon,  Portugal 

LJU  Ljubljana,  Yugoslavia 

LWI  Lv/iro,  Zaire 

MAL  Malaga,  Spain 

MAN  Manila,  Philippine  Is. 

MER  Merida,  Mexico 
MOS  Moscow,  RSFSR,  USSR 

NES  Northeastern  Seismological  Association,  Weston,  MA 

NOS  National  Ocean  Survey 

NOU  Noumea,  New  Caledonia 

OAX  Oaxaca,  Mexico 

OTT  Ottawa,  Canada 

PAS  Pasadena,  CA 

PDE  Preliminary  Determination  of  Epicenters 

PEK  Peking,  China 

PMG  Port  Moresby,  Papua 

PRA  Praha  (Prague),  Czechoslovakia 

QUE  Quetta,  Pakistan 

REY  Reykjavik,  Iceland 

RIV  Riverview,  New  South  Wales,  Australia 

ROM  Rome,  Italy 

SEA  Seattle,  WA 

SHI  Shiraz,  Iran 

SHL  Shillong,  India 

SSS  San  Salvador,  El  Salvador 

SIR  Strasbourg,  France 

STU  Stuttgart,  W.  Germany 

SYK  Lynn  R.  Sykes,  Lamont-Doherty  Geol .  Obsy.,  Palisades,  NY 
TAC  Tacubaya,  Mexico 
TEH  Teheran,  Iran 
TOK  Tocklai,  India 
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TRI  Trieste,  Italy 

TRN  Trinidad,  Trinidad,  W.  I. 

UCC  Uccle,  Belgium 

UPP  Uppsala,  Sweden 

USE  United  States  Earthquakes 

WAR  Warsaw,  Poland 

WEL  Wellington,  Mew  Zealand 

ZUR  Zurich,  Switzerland 
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